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Abstract 
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SCHOOL OF PHYSICS AND ASTRONOMY 

Doctor of Philosophy 

PHENOMENOLOGY OF THE MINIMAL B - L EXTENSION 
OF THE STANDARD MODEL AT THE LHC 

by Lorenzo Basso 

A well-motivated framework to naturally introduce neutrino masses is the B — L 
model, a U{1) extension of the standard model related to the baryon minus lepton 
gauged number. Besides three right-handed neutrinos, that are included to cancel the 
anomalies (thereby naturally providing neutrino masses), this model also encompasses a 
complex scalar for the spontaneous symmetry breaking of the extended gauge sector and 
to give mass to the Z' boson. We present the phenomenology, the discovery potential 
at the LHC, and the most up-to-date experimental and theoretical limits of the new 
particles in this model. In the gauge sector, a Z' boson is present. We study its 
properties (i.e., production cross sections, branching ratios, total width), showing that 
it is dominantly coupled to leptons. We also present a detailed discovery power study 
at the LHC and at Tevatron for the Z' boson. In the fermion sector, after implementing 
the see-saw mechanism, we end up with three heavy neutrinos. We show that they can 
be long-lived particles (therefore providing displaced vertices in the detector), and that 
they can induce spectacular multi-lepton decays of the Z' boson. We also study the full 
signature pp ^ Z' ^ i^h^h^ and present a parton level and a detector level analysis for 
the tri-lepton decay mode of the Z' boson via heavy neutrinos. In the gauge sector, 
the two Higgs fields mix. We derive the unitarity bound and the constraints from the 
renormalisation group equations study. In the allowed region of the parameter space, we 
delineate the phenomenology of the Higgs bosons and we show characteristic signatures 
of the latter, at the LHC, involving the Z' boson and the heavy neutrinos. 
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Chapter 1 

Introduction 



A new era in particle physics has started. The Large Hadron ColHder (LHC) at CERN, 
Geneva, is finally taking data and a huge effort from both theoretical and experimental 
communities is required to meet the challenge. 

The standard model (SM) ^ of electroweak (EW) and strong interactions ^ is once 
again going to be severely tested, as well as many of its extensions that have been 
proposed to cure its flaws. The observed pattern of neutrino masses [8; 9], the existence 
of dark matter [10] and the observed matter-antimatter asymmetry [10; 11] are the most 
severe evidences the SM fails to explain. Deep criticism can be moved against the SM 
also from the theoretical point of view. The so-called 'hierarchy problem' (see, e.g., [12] 
and references therein) is an example. The inadequacy of the SM is reinforced by the 
fact that one of its fundamental component is still missing: the Higgs boson. 

It is widely accepted that the SM ought to be extended, but no one knows if the 
proper way has already been explored in the literature. A joint collaboration is therefore 
needed between the experimental and the theoretical communities. 

The project of this Thesis has been guided by these principles. The aim is to fill 
some of the gaps in the overall preparation towards real data, as well as to interact 
proficiently with experimentalists. An extension of the SM has been systematically 
studied, from the definition of its parameter space to the collider signatures, leading to 
some novel and exciting possibilities for the LHC to shed light on. Within the ambition 
of a complete study, in all its aspects, the experimental help has been fundamental to 
efficiently concentrate on aspects of actual interest in a way useful for both communities. 

The main motivations ^ for the extension of the SM that we will describe concern the 
lack of a natural explanation for the observed pattern of neutrino masses, the unknown 
origin of a global and not anomalous accidental ^7(1) symmetry in the SM (related to 
the baryon minus lepton (B — L) quantum numbers), and the absence of any observation 
of a fundamental scalar degree of freedom (the Higgs boson). 

'^For reviews, excellent textbooks exist, see Refs. [6; 7]. 

^We do not consider here the inclusion and description of gravity at the quantum level, often kept 
aside of the traditional SM framework. 

Other fundamental reasons to extend the SM being the inclusion of the gravity, the lack of a dark 
matter candidate, the dark energy, the hierarchy problem, the matter-antimatter asymmetry. 
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Besides the global U {1)b-l symmetry of the SM, that can be thought as an acciden- 
tal symmetry (not imposed, as the custodial symmetry [13]), the previous issues can be 
tackled separately. The neutrino masses and mixing angles are traditionally explained 
by means of the so-called 'see-saw mechanism' [14; 15; 16; 17; 18], simply including 3 
right-handed (RH) neutrinos and imposing that their effective mass term is much bigger 
than the Dirac one, that couples the traditional left-handed (LH) neutrinos and the RH 
neutrinos to the Higgs field. The right scale for the masses of the 3 SM-like neutrinos is 
then reproduced by an effective mass term for a RH neutrino of 0(10^^ ^ 10^^) GeV. 

By now, the Higgs boson is the only undetected particle of the SM, and its properties 
are still unknown. Therefore, it is not unreasonable to think of modifications of the scalar 
sector that are still compatible with experimental constraints. The most economical way 
to modify the scalar sector of the SM is to include one (or more) scalar singlets, either 
real [19; 20; 21; 22; 23] or complex [24], whose phenomenology at hadronic and leptonic 
colliders has been studied in great detail, as well as their impact on precision observables. 
In fact, the latter are able to constrain the viable parameter space of the extended Higgs 
sectors (see, e.g., [2; 22; 25] and references therein). 

Augmenting the scalar sector only does not provide an explanation for the observed 
pattern of the neutrino masses and mixing angles. Following a bottom-up approach, 
this Thesis discusses a well motivated framework that remedies at once such flaws of the 
SM: the B — L model [26; 27; 28; 29]. This is a triply-minimal extension of the SM. It is 
minimal in the gauge sector, in which a single U{1) factor is added, related to the B — L 
number, by simply promoting to local the already existing U{1)b-l global symmetry of 
the SM. It is minimal in the fermion sector, in which a SM singlet fermion per generation 
is added, to cure the new C/(1)b~l related anomalies. These fermions can naturally be 
interpreted as the RH neutrinos. It is minimal in the scalar sector, in which a complex 
neutral scalar singlet is added to spontaneously break the new U{1) symmetry, and at 
the same time to give to the new gauge boson a mass (to evade present experimental 
bounds). The two latter points, once the U{1)b-l symmetry is spontaneously broken, 
naturally provide a dynamical implementation of the see-saw mechanism, explaining the 
neutrino masses. As we will see, the remnant degree of freedom of the new complex 
scalar severely impinges in the phenomenology of the scalar sector. 

As mentioned, we follow a bottom-up approach, being mainly interested in studying 
the phenomenology of the model at colliders. In the light of this, we are not here 
interested in seeking an embedding of this model into grand unified theories (GUT). 
Direct consequence is that the gauge coupling related to the extra U{1) factor, controlling 
the interaction of the Z' boson, is a free parameter. Even though we are not interested 
in studying the details of it, we are aware that extra U{1) factors appear naturally in 
GUT scenarios, conceptually motivating our minimal extension. 

The general model we introduce is a one-dimensional class of U{1) extensions of the 
SM, in which each element is characterised by the properties of the new gauge boson 
associated to the extra U{1) factor. In all generality, the latter couples to fermions 
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proportionally to a linear combination of the standard liypercharge (Y) and the B — L 
number. The coefficients of this linear combination are 2 ^ gauge couplings, that in 
our bottom-up approach are free parameters. A variety of standard benchmark models 
is then recovered. For example, the sequential standard model (SSM), the C/(l)/j, the 
U{\)^, and the 'pure' B — L model, all anomaly-free with the SM particle content 
augmented by the RH neutrinos. The latter one is the benchmark model on which we 
will focus our numerical study. 

The 'pure' B — L model is identified by the fact that the extra gauge boson, or Z'^_j^, 
couples to fermions proportionally to their B — L number only. On the one side, this 
directly implies a vanishing Z — Z' mixing (at the tree-level), that is consistent with 
the existing tight constraints on such mixing, compatible with a negligible value [33]. 
Moreover, the B — L charge does not distinguish the chirality, i.e., the LH and the RH 
degrees of freedom of the same fermion have the same B — L quantum numbers. Hence, 
Dirac fermions have just vectorial couplings to the Z'^_^ boson: 



9z' = 9z' 



9l = '-^^ = iB-L)g[, 
g^, = = 0. 



(1.1) 



On the contrary, the mass-eigenstates for the neutrinos, after the see-saw mechanism, 
are Majorana particles. Therefore, they have pure axial couplings to the Z'^_^ boson 
[34]. 

As a consequence, we decided not to study the asymmetries of the decay products 
stemming from the Z'^_^ boson, given their trivial distribution at the peak. However, 
asymmetries can be important in the interference region, especially just before the Z' 
boson peak, where the Z — Z' interference will effectively provide an asymmetric distri- 
bution somewhat milder than the case in which there is no Z' boson. This is a powerful 
method of discovery and identification of a Z' boson and it deserves future investigations. 

It is important to emphasise that this is a TeV scale extension of the SM. This means 
that the U{1)b-l breaking vacuum expectation value (VEV) is of O(TeV). Hence, also 
the new particles will have masses at the TeV scale. In particular, we will consider 
both the Z' boson and the new Higgs particle with masses in the 0.1 10 TeV range. 
Regarding the neutrinos, with a suitable choice of the parameters, the 3 light neutrinos 
will have masses in the sub-eV range, while the 3 heavy ones will have masses of O(IOO) 
GeV. 

An attractive feature of the B — L model is that a successful realisation of the mech- 
anism of baryogenesis through leptogenesis (to explain the observed matter-antimatter 
asymmetry) might be possible in this model, as shown in Refs. [35; 36]. It is very 



n + l 

^In extensions of the SM with n (7(1) factors, the Abehan gauge group mix. For consistency, k 

fe=i 

gauge couplings appear. In our model, where n = 1, we will have 2 more gauge coupling besides the SM 
one, g [30; 31; 32]. More details are given in section 2.1. 



4 



1. INTRODUCTION 



suggestive that this realisation takes place at the TeV scale. However, its correct imple- 
mentation results in constraints on the neutrino masses, both on the light ones and on 
the heavy ones, and on their CP-violating phases. 

As already stated, in this Thesis we are mainly concerned with the impact of the 
model at colliders, on which the above phases do not play any role. Also, to maximise and 
highlight the interesting patterns that can be observed, and for illustrative purposes, we 
decided to not include here the above constraints on the neutrino masses. To reconcile, 
the results that we will show can be thought to be valid for at least one generation. 

Bearing this in mind, we will consider all the heavy neutrinos as degenerate and 
with masses that are free parameters, varying them in the 50 500 GeV range. Hence, 
in most of the cases, smaller than Z' boson mass. We will show that the Z' boson in 
this model, with TeV scale heavy neutrinos, can decay into pairs of the latter. Another 
difference with respect to the traditional Z' literature is, therefore, that the Z'^_^ boson 
branching ratios (BRs) are not fixed. Neither its intrinsic width is, being the gauge 
coupling as well a free parameter. 

The presence of new coupled matter, the heavy neutrinos, has important phenomeno- 
logical consequences. The possibility of the Z' boson (and of the Higgs bosons, as we 
will show) to decay into pairs of them, will provide new and exciting signatures. It is 
worth to briefly mention that the peculiar decays of the Higgs bosons into pairs of heavy 
neutrinos, or into pairs of Z' bosons, is a distinctive signature of this model, offering 
the chance to distinguish it from the plethora of the otherwise identical, concerning the 
scalar sector, singlet extensions of the SM in the literature. 

In this Thesis we decided to study the details of the gauge and fermion sectors. Their 
mutual interactions are fully included in the Z' decay into pairs of heavy neutrinos. 
Altogether, this decay provides new and spectacular multi-lepton signatures of the Z' 
boson. One of the main results of this Thesis is the study of one of them, the tri-lepton 
decay mode (i.e., when the Z' decays into exactly 3 charged leptons and other particles, 
such as jets and/or missing energy), together with the related backgrounds. We will 
present a parton level strategy for reducing the latter in order to isolate the signal, that 
will be validated at the detector level. 

It is very interesting that heavy neutrinos can be long-lived particles. For the exper- 
imental community, this model represents a test laboratory to study trigger efficiencies 
and detector resolutions for such a clear signal of physics beyond the SM (see Ref. [37] 
for a review). Moreover, from the simultaneous measurement of both the heavy neutrino 
decay length (from the displacement of its secondary vertex in the detector) and mass 
(we will show that this task is achievable in the tri-lepton signature) one can estimate 
the absolute mass of the parent light neutrino, for which at present, only limits exist. 
Altogether, this realises a spectacular and very peculiar link between very high energy 
and very low energy physics. 

It shall be noticed that a first insight on the phenomenology of the pure B — L model 
at the LHC was presented in Ref. [38]. 
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The general structure of each chapter is to present the constraints for the new pa- 
rameters in each sector and to study the phenomenofogy related to the new particles in 
the B — L model, with respect to the SM. They concern the parameters in the gauge 
sector (i.e., the Z' mass and the gauge coupling g'-^), in the fermion sector (i.e., the light 
and the heavy neutrino masses), and in the scalar sector (i.e., the scalar masses and 
mixing angle a). 

Both experimental searches and theoretical arguments can give informations about 
the regions of the parameters that are allowed. 

From the experimental side, past and existing facilities have looked for new particles, 
directly and indirectly, still with no success. This poses constraints on the parameters 
to evade limits. Tevatron, still accumulating data, is updating its constraints from 
direct searches, but for many cases the most stringent bounds come from the searches 
performed at LEP. 

The LEP searches are still setting some of the most stringent bounds in all the sectors. 
Key factors for their success are manifold: the strength of the interaction between the 
Higgs boson and the gauge bosons, in the SM; the clean experimental environment of 
such a leptonic collider; the enormous amount of data collected at the SM Z boson peak. 
They allow, in turn, to set stringent limits also in the B — L model, on the scalar masses 
(as a function of the scalar mixing angle), on the Z' mass and g'^ coupling, and on the 
heavy neutrino masses. 

The studies performed with the data collected by LEP (at the SM Z boson peak 
and beyond) are the so-called precision tests of the SM, since certain observables are 
measured with per mil accuracy. These observables are usually referred to as preci- 
sion observables. The level of agreement with the SM predictions is remarkable. This 
agreement results in tight constraints for most of the SM extensions, since they gener- 
ally predict new contributions to the precision observables (usually to be evaluated at 
loop- level). In studying these limits, we did not pursue a complete beyond the tree-level 
study; we rather refer to the literature for that. In fact, we decided to concentrate on 
the study of the collider signatures and to leave the precision test analysis for future 
investigations. 

Tevatron is currently collecting data and updating its analyses regarding the SM 
Higgs boson and Z' boson masses (of a discrete choice of models) from direct searches. 
However, for the latter, the Z'^_^ is not part of its traditional literature. Hence, the 
most up-to-date limits on the Z'q_^ mass at Tevatron, presented in section 3.1.1.1, are 
a novel contribution (firstly in Ref. [39]). 

From the theoretical side, the model has to be well-defined. Several meanings apply 
to this concept, the main ones regarding the scalar potential to be bounded from below 
and with a stable vacuum solution, and the perturbative expansion to hold. The most 
powerful techniques require to go beyond the tree-level order, by studying the evolution 
of the renormalisation group equations (RGEs) of the parameters of the model. The 
imposition of specific relations results in constraints on the yet free parameters. At the 
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same time, strong constraints, especially on the scalar masses, come from requiring the 
model to be unitary. 

We present the analysis of the theoretical constraints on the new parameters of the 
model, namely the unitarity bound and those coming from the study of the RGEs. The 
unitarity assumption [3] relies on demanding that the cross sections for the scattering 
of all the particles in the model are unitarised, i.e., with an occurring probability less 
than one. The tightest bounds come from the scattering of the longitudinally-polarised 
gauge bosons, and they constrain mainly the scalar masses. 

The RGE constraints [40] come from the analysis of the model at one loop, assuming 
the evolved parameters do not hit any Landau pole or spoil the theory in the scalar 
sector, e.g., the vacuum of the model. The latter request is known as the vacuum 
stability condition, while the former is the triviality condition. 

The one-loop RGEs we collected from the literature have been completed with direct 
calculations of the top quark one and of those regarding the scalar sector. The complete 
set is collected in Appendix B. As introduced in section 2.1, all the parameters of the 
model are running parameters, and the requirement of some to vanish at the EW scale 
just sets the boundary conditions of their evolution. In our case, the mixing parameter 
in the gauge sector, i.e., g, is set to vanish at the EW scale, defining the 'pure' B — L 
model, argument of the numerical analysis of this Thesis. However, its evolution must 
be taken into account. 

The numerical study of the collider phenomenology in this Thesis has been done with 
the CALCHEP tool [41] (where not specified otherwise). The model has been implemented 
therein with the lanhep module [42]. The availability of the model implementation 
into CALCHEP in both the unitary and t'Hooft-Feynman gauges allowed us to perform 
powerful cross-checks to test the consistency of the model itself. 

The RGE have been solved and studied in mathematica 5.0 [43]. 

1.1 Outline 

The Thesis is organises as follows. Chapter 2 introduces the model and contains a 
description of its consequences, i.e., spontaneous symmetry breaking and the particle 
spectrum. 

Chapter 3 presents the study of the gauge sector. It starts with the constraints on 
the free gauge coupling and on the Z' mass and it studies the Z' boson production cross 
sections and 2-body decays BRs. Finally, the discovery reach at Tevatron and LHC of 
the Z' boson in Drell-Yan processes is shown. 

Chapter 4 describes the study of the fermion sector. After reviewing the limits on 
neutrino masses, a detailed description of the decay properties of the heavy neutrinos 
is given. Their main production mechanism is in pair via the Z' boson, and, bringing 
pieces together, we show that they can induce new spectacular decay patterns for the 
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latter. Among them, we perform detailed parton and detector level analyses of one 
specific pattern, the tri-lepton decay mode. This is the main result of this Thesis. 

Chapter 5 analyses the theoretical bounds, i.e., unitarity, triviality, vacuum stability, 
and some experimental constraints on the scalar sector, presenting also production cross 
sections and BRs for the Higgs bosons, as well as cross sections for some signatures 
involving decays of the Higgs bosons into pairs of Z' bosons and heavy neutrinos. 

In appendix A is explained the implementation of the model in the numerical tools 
used and some related technicalities. We also list the complete set of Feynman rules 
for the pure B — L model. In appendix B are collected and described the RGEs of the 
model. 



Chapter 2 

The Model 



In this chapter we introduce the model that we wih study in this Thesis. We will 
introduce the notation and highlight the differences with respect to the SM. For the 
motivations, we refer to the introduction. 

The model under study is a triply-minimal extension of the SM. It is minimal in 
the gauge sector, where the SM gauge group is augmented by a single U{1) factor, 
related to the baryon minus lepton {B — L) gauged number. It is minimal in the scalar 
sector, where a complex scalar singlet is required for the spontaneous breaking of the 
new U (1) gauge factor, that we assume to take place at the TeV scale. Finally, it is 
minimal in the fermion sector, in which a SM singlet RH fermion per generation is 
introduced. Its B — L charge, in particular, is chosen to cure the new gauge and mixed 
i7(l)— gravitational anomalies, related to the U{1)b-l gauge symmetry factor [29]. It 
is important to stress that three and only three RH fermions are introduced (one per 
generation, as required by the anomally cancellation conditions), that can naturally be 
identified with the RH neutrinos. 

The extra Higgs singlet needs to be charged under B — L to trigger the spontaneous 
symmetry breaking of the related U{1) factor. Its charge is then fixed by the gauge 
invariance of the Yukawa term that couples it to the RH neutrinos. Altogether, this 
provides a natural and dynamical framework for the implementation of the type-I see- 
saw mechanism, contrary to its usual realisation through an effective mass term for the 
RH neutrinos. Hence, in this model, neutrinos are naturally massive particles and their 
masses are generated through the Higgs mechanism. 

2.1 A minimal U{1) extension 

As mentioned, we present here the general model under study, a minimal U{1) extension 
of the SM, related to the B — L number. The classical gauge invariant Lagrangian, 
obeying the SU{3)c x SU (2)^ x [/ (l)y x U (1) b-l gauge symmetry, can be decomposed 
as: 

if = if, + ^YM +^f+^Y, (2.1) 
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where the terms on the right hand side identify the scalar part, the gauge (or Yang-Mills) 
part, the fermion part and the Yukawa part, respectively. 

2.1.1 Scalar sector 

The spontaneous symmetry breaking of the new extra U{1) factor in the gauge group 
is achieved by introducing in the particle content a further scalar field, that is required 
to acquire a VEV at the TeV scale. To keep the realisation of the Higgs mechanism 
minimal, only one extra degree of freedom (beside the one that will become a physical 
particle) is required, to give mass to the new neutral gauge field associated to the extra 
U{1) factor. Hence, the minimal choice for the new scalar is a complex singlet, that we 
will denote as x- 

The scalar Lagrangian is 

= (D'^H)^ D^H + (i?^x)^ D^x - V{H, x) , (2.2) 
with the scalar potential given by 

ViH,x) = m'^H^H + \x\^ + ( H^H |xP 

= m'^H^H + fi^\x\'^+XiiH^Hf + X2\x\'^+hH^H\xW (2.3) 

where H and x &re the complex scalar Higgs doublet and singlet fields. Regarding the 
new B — L charges for these scalar fields, we take and +2 for H and X; respectively. 
Notice that the charge of the x field has been chosen to ensure the gauge invariance 
of the fermion sector for the minimal model under discussion. A non-minimal matter 
content would lead to a different requirement for the charge of x (see, e.g., Ref. [44]). 

2.1.2 Yang-Mills sector 

Moving to the ^ym^ the non-Abelian field strengths therein are the same as in the SM 
whereas the Abelian ones can be written as follows: 




where 



F^,,. = d^B,-d,B^, (2.5) 
F' = d^Bl-d^B'. (2.6) 
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We define 2?^ and B'^ as the U{1)y and U{1)b-l gauge fields, respectively. In this field 
basis, the covariant derivative is ^: 

D^ = d^ + igsT'^G^^ + i^T^iy/ + igiYB^ + i{gY + g[YB-L)B'^ . (2.7) 

In our bottom-up approach, we will not require gauge unification at some specific, 
yet arbitrary, energy scale ^, fixing the conditions for the extra gauge couplings at that 
scale. Therefore, in this model, the gauge couplings g and g'l are free parameters. To 
better understand their meaning, let us focus on eq. (2.7). This form of the covariant 
derivative can be re-written defining an effective coupling and an effective charge 

gEY^ = gY + g'^YB^L. (2.8) 

As any other parameter in the Lagrangian, Ij and g'^ are running parameters [30; 31; 
32], therefore their values ought to be defined at some scale. A discrete set of popular 
Z' models (see, e.g., Refs. [45; 46]) can be recovered by a suitable definition of both g 
and g'l- 

We will focus in our numerical analysis on the 'pure' B — L model, that is defined 
by the condition ^{Qew) = 0, i.e., we nullify it at the EW scale. This implies no 
mixing at the tree level between the B — L Z' and the SM Z gauge bosons. Other 
benchmark models of our general parameterisation are for example the Sequential SM 
(SSM), defined hyY^ = Y (that in our notation corresponds to the conditions g'l = Q 
and g = gi aX, the EW scale) and the U{l)ii model, for which RH fermion charges vanish 
(that is recovered here by the condition Ij = —2g'^ at the EW scale). 

It is important to note that none of the models described so far is orthogonal to 
the U{1)y of the SM, therefore the RGE running of the fundamental parameters, 'g and 
g']^, will modify the relations above. The only orthogonal U{1) extension of the SM is 
the 'S'O(10)-inspired' U{1)^ model, that in our notation reads g = — fg'i. Although 
the 5 and g'^ couplings run with a different behaviour, the EW relation 'g/ g'l = —4/5 is 
preserved (at one loop) at any scale. 

2.1.3 Fermion sector 

The fermion Lagrangian (where k is the generation index) is given by 

3 

k=l 

+ilkLltiD'^lkL + iekRlfiD'^ekR + ii^7M^''^fe/?) > (2-9) 

^In all generality, Abelian field strengths tend to mix. The kinetic term can be diagonalised with a 
GL{2,R) transformation, leading to the form of the covariant derivative in eq. (2.7), where the mixing 
between the two U{1) factors becomes once again evident [31; 32]. For the curios reader, we refer to 
Ref. [29] for a more comprehensive explanation. 

■^The gauge unification could in principle happen in a multi-step process, passing through interme- 
diate gauge groups until the final single group. 
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where the fields' charges are the usual SM and B — L ones (in particular, B — L = 1/3 
for quarks and —1 for leptons with no distinction between generations, hence ensuring 
universality). The B — L charge assignments of the fields as well as the introduction 
of new fermion RH heavy neutrinos {i'rs, charged —1 under B — L) are designed to 
eliminate the triangular B — L gauge anomalies of the theory. Regarding the new scalar 
Higgs field [x] the charge +2 under B — L is chosen to ensure the gauge invariance of 
the model (for a detailed discussion, see Ref. [29]). 

Therefore, the B — L gauge extension of the SM gauge group broken at the TeV scale 
does necessarily require at least one new scalar field and three new fermion fields which 
are charged with respect to the B — L group. 

2.1.4 Yukawa interactions 

Finally, the Yukawa interactions are 

■^Y = -y'-kOjidkRH - VjkqjLUkRH - y^ljiekRH 

-y^JJ^i^kRH-yfiJ^ji^kRX + i^-c., (2.10) 

where H = ia^H* and k take the values 1 to 3, where the last term is the Majorana 
contribution and the others are the usual Dirac ones. Without loss of generality, we work 
on the basis in which the RH neutrino Yukawa coupling matrices, y*^, are diagonal, real, 
and positive. These are the only allowed gauge invariant terms ^. In particular, the last 
term in eq. (2.10) couples the neutrinos to the new scalar singlet field, Xi and it allows 
for the dynamical generation of neutrino masses, as x acquires a VEV through the Higgs 
mechanism. Due to this, in the B — L model, the neutrinos couple to the scalar sector, 
in particular the right-handed (RH) ones couple strongly to the singlet scalar. 

Neutrino mass eigenstates, obtained after applying the see-saw mechanism, will be 
called vi (with / standing for light) and Uh (with h standing for heavy), where the 
first ones are the SM-like ones (see section 2.2.3). With a reasonable choice of Yukawa 
couplings, the heavy neutrinos can have masses m^^ ~ O(IOO) GeV. 

2.2 Symmetry breaking 

We generalise here the SM discussion of spontaneous EW symmetry breaking (EWSB) 
to the more complicated classical potential of eq. (2.3). To determine the condition for 
V{H, x) to be bounded from below, it is sufficient to study its behaviour for large field 
values, controlled by the matrix in the first line of eq. (2.3). Requiring such a matrix to 
be positive definite, we obtain the conditions: 



4AiA2-Ai > 0, (2.11) 
Ai,A2 > 0. (2.12) 
Notice that an effective mass term sucfi as MiynYvR is now forbidden by the gauge invariance. 
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If the conditions of eqs. (2.11) and (2.12) are satisfied, we can proceed to the min- 
imisation of 1^ as a function of a constant VEV for the two Higgs fields. In the Feynman 
gauge, we can parametrise the scalar fields as 



H 




—i{w^ — inP') 
V + [h + iz) 



1 



X 



V2 



{x + {h! + iz')) 



(2.13) 



=F inP' , z and z' are the would-be Goldstone bosons of , Z and 



where 

Z\ respectively. Making use of gauge invariance, for the minimisation of the scalar 
potential, it is not restrictive to assume 



1 







V2 



(2.14) 



with V and x real and non-negative. The physically most interesting solutions to the 
minimisation of eq. (2.3) are obtained for v and x both non-vanishing: 



A1A2 — 



^3 



AiAo-^ 



(2.15) 
(2.16) 



Physical solutions must have both > and > conditions satisfied. Eq (2.11) 
implies the denominators are always positive, so we just need to require the numerators 
to be positive too. In terms of the parameters in the Lagrangian, this in turn means 



Aom 



Ai//' 



< 



< 



A3 2 



^3 2 

—m 



It is interesting to note that there are no solutions for both fj?,m'^ > 0, irrespectively of 
A3. Solutions with (^^ > 0, rrP < 0) or (^^ < 0, rrP > 0) are allowed only for A3 < 0. 
Solutions for (^^ < 0, m? < 0) are allowed for both signs of A3. 



2.2.1 Gauge eigenstates 

To determine the gauge boson spectrum, we have to expand the scalar kinetic terms as 
for the SM. We expect that there exists a massless gauge boson, the photon, whilst the 
other gauge bosons become massive. The extension we are studying is in the Abelian 
sector of the SM gauge group, so that the charged gauge bosons will have masses 
given by their SM expressions, being related to the SU{2)l factor only. Using the 
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unitary-gauge parameterisation, the kinetic terms in eq. (2.2) become: 



-d^hd^h + -{h + vf (0 1) gW^^Oa + 9iB^ + gB''^ 



z o L 



2 / 



+ {gW^''-giB>^-gB'>^f 



(2.17) 



and 



^-^^n- - \d>^h'd^h' + ]^{h' + xf{g[2B'^^f, 



(2.18) 



where we have taken Y^~^ = 2 in order to guarantee the gauge invariance of the Yukawa 
terms (see eq. (2.10)). In eq. (2.17) we can recognise immediately the SM charged gauge 
bosons W^, with My/ = gv/2 as in the SM. The other gauge boson masses are not so 
simple to identify, because of mixing. In fact, in analogy with the SM, the fields of 
definite mass are linear combinations of i?'^, W^^ and B'^. The explicit expressions are: 



B^' \ 






cos 1? 




h 




sin "!?, 












— sin t?^ cos sin i9^„ sin 
cos "dw cos t?' — cos i9w sin i9' 
sin cos t?' 





MM 






/ 





(2.19) 



with -f < 1?' < f , such that: 



and 



tan2i?' 



25vV+5f 



g2 + lQ[^fg'2_g2_g2 



(2.20) 



Ma 
Mz,z' 



where 



0, 



2 V g^+gi / 



sin2'!9'Y^^9i 



sin 2t?' 



25vV+5f 



(2.21) 
(2.22) 

(2.23) 



[g^ + mfg\^ -g^- glf + {2gf{g^ + gl) 

The LEP experiments [33] constrain ^ 10~^. Present constraints on the VEV x (see 
section 3.1.1) allow a generous range of g. 

In the pure B — L model, that we remind is defined by the condition 5^ = 0, one can 
easily see, from eq. (2.23), that also the mixing angle vanishes, implying no mixing, 
at the tree level, between the SM Z and the new Z' gauge bosons, identically satisfying 
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the LEP bound. The Z and Z'q_^ masses then simphfy, respectively, to 

Mz = (2.24) 
Mz'^_^ = Hx. (2.25) 

2.2.2 Scalar eigenstates 

To compute the scalar masses, we must expand the potential in eq. (2.3) around the 
minima in eqs. (2.15) and (2.16). We denote by hi and /12 the scalar fields of definite 
masses, nih^ and rrih^, respectively, and we conventionally choose m^^ < rnj^^. After 
standard manipulations, the explicit expressions for the scalar mass eigenvalues and 
eigenvectors are 

^li = Ai?;2 + Aax^ - ^(Ai^;^ - \^x^f + {\^xvf , (2.26) 
ml^ = W + A2x2 + V(Ai^;2 - A2x2)2 + {\^xvY , (2.27) 




cos a — sm a 
sin Q cos a 




(2.^ 



where -f < a < f fulfils ^ 



sin2Q = = , (2.29) 

V(Ai^;2 _ A2x2)2 + (A3x^;)2 ^ ^ 

o Ai^^ - Asa;^ , . 

cos 2a = — , (2.30) 

V(Ai^;2 - A2x2)2 + (A3x^;)2 

Eqs. (2.3), (2.28), and (2.29) describe a rather general scalar sector as in the mini- 
mally extended SM with scalar singlets, in which the light(heavy) Higgs boson couples 
to SM matters proportionally to cosa(sina). Proper differences come from having new 
coupled matter, as we will explain in chapter 5. Notice that the light (heavy) Higgs 
boson couples to the new matter content (heavy neutrinos and Z' boson) with the com- 
plementary angle, sina(cosa), respectively, as the latter is directly coupled only to the 
scalar singlet x- 

It is interesting to study the decoupling limit (when a — t- 0). In this limit, hi is 
purely the SM-like boson, while /i2 does not couple to SM particles, making its discovery 
impossible at the LHC. We will further comment on this in section 5.2. The other 
possible decoupling limit, for a — )• ^ (where /i2 is the SM-like Higgs boson), is excluded 
by present experimental constraints [2]. 

For our numerical study of the extended Higgs sector, it is useful to invert eqs. (2.26), 
(2.27) and (2.29), to extract the parameters in the Lagrangian in terms of the physical 



^In all generality, the whole interval < a < 27r is halved because an orthogonal transformation 
is invariant under a — ^ a + vr. We could re-halve the interval by noting that it is invariant also under 
a — >■ —a if we permit the eigenvalues inversion, but this is forbidden by our convention m^^ < . 
Thus a and ~a are independent solutions. 
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quantities mh^, rrih^ and sin 2a: 

2 2 

mf mr 
Ai = -cos2a) + ^(1 + cos2a), (2.31) 

2 2 

mf mf 
As = ^(1 - cos 2a) + ^(1 + cos 2a), (2.32) 

J . (2.33) 

2.2.3 Fermion eigenstates 

The last line in eq. (2.10) contains the Dirac and the Majorana mass terms for the 
neutrinos, respectively. In contrast to the usual effective implementation of the see-saw 
mechanism, the last term of eq. (2.10) is a proper Yukawa interaction with the new 
Higgs singlet. 

To extract the neutrino masses we have to diagonalise the neutrino mass matrix from 
eq. (2.10): 

\ ruD M 



(2.34) 



where 

mo 



^v, M = V2y^'x, (2.35) 
v2 

where x is the VEV of the x held. This matrix can be diagonahsed by a rotation about 
an angle a^, such that: 

tan2aj. = -— ^. (2.36) 

For simplicity we neglect the inter-generational mixing so that neutrinos of each 
generation can be diagonahsed independently. We also require that the neutrinos be 
mass degenerate. Thus, i'l,r can be written as the following linear combination of 
Majorana mass eigenstates ui^h '■ 



J^L \ / cos ai/ — sm a,y \ I 



(2.37) 



whose masses are respectively 



m 



2 



= (2-38) 
m^, ^ M. (2.39) 

We can now rewrite eq. (2.36) in term of the physical masses: 



777, 

tan2ai. ^ -2./— . (2.40) 
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For m^^ ~ 0.1 eV and m,y^ ~ 100 GeV, tana^ ~ 10^^. Hence, I'lii'h) are mostly 
the LH(RH) neutrinos. One consequence is that the heavy neutrinos can be long-hved 
particles, as we will show in section 4.2.1. Also, being mostly RH, the heavy neutrinos 
are strongly coupled to the heavy(light) Higgs boson for small(big) scalar mixing angle 
a, while the SM-like ones (being mostly LH) are not. 



2.3 Summary 

We summarise here the matter content of the pure B — L model, argument of the 
numerical analysis of this Thesis. We list the new independent parameters with respect 
to the SM, with the used symbols. 



2.3.1 Gauge sector 

We remind the reader that the pure B — L model is defined by ^ = 0. Therefore, the 
only independent parameters in the gauge sector are the g'^ coupling and the Z' boson 
mass Mz' ^ The B -L breaking VEV X is then defined as 

2.3.2 Fermion sector 

All the neutrino masses are independent parameters of the model. For simplicity, we ne- 
glect the mixing within the generations and, furthermore, we take all the light neutrinos 
and all the heavy neutrinos to be degenerate in mass: 



nij^i = m^2 = m^3 = m^^ , (2.42) 

(2.43) 



The mixing angle a^, is then fixed, as in eqs. (2.36) and (2.40). 



2.3.3 Scalar sector 

The Higgs boson masses, m/j^ and m/j2, and mixing angle, a, are the only new parameters 
in this sector. Notice that, by convention, m/i^ < m/jj and — f < a < f • 

Table 2.1 summarises the charge assignments to fermion and scalar fields, while 
table 2.2 shows the fields' labels and masses. 



^Hereafter, when not confusing, we will call simply Z' boson the Z'g_j^ (the gauge boson of the pure 
B ~ L model), and Mz' its mass. 
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Table 2.1: y and B — L quantum number assignation to chiral fermion and scalar fields. 



Name 


Quarks 


Leptons 


Neutrinos 


Hig 


gses 


^|, 


q 


i 






hi 


h2 


Mass 


nig 


mi 











Table 2.2: Mass eigenstates. 



Chapter 3 

Gauge sector 



In this chapter we discuss the gauge sector of the pure B — L model, starting by de- 
lineating the viable parameter space for its new independent parameters, as listed at 
the end of chapter 2. Both experimental searches and theoretical arguments can give 
informations about the regions of the parameters that are allowed. 

Prom the experimental side, the leading constraints come from Tevatron and LEP, 
in two different and complementary kinematic regions, and they are presented in sec- 
tion 3.1.1. Tevatron, still accumulating data, is updating its constraints from direct 
searches, and its constraints are tighter for low Z' boson masses. Por values of the mass 
beyond its kinematic reach, the most stringent bounds come from the searches performed 
at LEP. 

Prom the theoretical side, the study of the triviality bound (from a RGE analysis) of 
the parameters of the model gives an upper bound on the gauge coupling. Notice that 
the equations pertaining to the gauge sector decouple from the rest. Therefore, they can 
be studied independently. We will present their analysis in section 3.1.2. 

After the viable parameter space is individuated, in section 3.2 is studied the phe- 
nomenology at the LHC of the new particle in the gauge sector of the pure B — L model, 
the Z' boson. We present the results of our investigation by delineating its properties 
(i.e., production cross sections, intrinsic width, BRs). The possibility of the Z' boson to 
decay into pairs of heavy neutrinos is certainly the most interesting of its features. Also, 
a parton level discovery potential and exclusion power study at the LHC is presented in 
section 3.2.2. Pinally, we draw the conclusions for this chapter in section 3.3. 

The main results in this chapter are the presentation of the experimental limits on 
the Z' mass from Tevatron (published in [39]), the discussion of the theoretical bounds 
(published in [40]) and the analysis of the Z' boson properties (published in [47]) and 
discovery potential at the LHC (published in [39]). 

3.1 Constraints 

In this section is delineated the viable parameter space for the new independent param- 
eters in the gauge sector, i.e., g'^ and Mz'- 
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First, we present the experimental constraints using the data from Tevatron, and we 
summarise those from LEP. Then, we analyse the theoretical constraints, coming from 
the study of the RGEs. 

3.1.1 Experimental constraints 

A further neutral gauge boson is present in the B — L extension of the SM, the so-called 
Z' , whose mass, in eq. (2.19), is generically not predicted. In the pure version of the 
model, argument of this Thesis, eq. (2.19) simplifies to Afz^ ^ = 2(/^x. As specified in 
the introduction, we consider the gauge coupling free parameter of the model. 

That is, a lower bound on the Z' mass will, in general, be a function of the gauge 
coupling. 

Generally speaking, two ways exist to constrain the mass. The most simple one is 
to look for the decay products of the gauge boson in direct production. If no deviation 
from the background is observed, lower limits can be derived. Obviously, the direct 
production is limited by the kinematic reach of the experimental facility. In the case of 
a hadronic machine, the Parton Distribution Functions (PDFs) will also lower its reach. 
The most stringent bounds from direct detection are coming from Tevatron. They will 
be discussed in section 3.1.1.1. 

A leptonic machine, being LEP its latest realisation, offers other possibilities to look 
for neutral gauge bosons. In all generality, the latter will mix with the SM Z boson, 
affecting its mass. The LEP-I experiments have measured the mass and other properties 
of the SM Z boson rather accurately, confirming the SM predicted values at the per mil 
level [48; 49]. Therefore, the mixing of an extra neutral gauge boson to the SM Z boson 
is tightly constrained [33]. In the pure B — L model, this bound is evaded by definition, 
as the vanishing of the mixing gauge coupling 'g (that defines the model) also implies no 
mixing at the tree level in the gauge sector (see also section 2.2.1). 

Effective contact interactions, arising from an expansion in s/M^, (when M^, > s), 
are a complementary way for setting bounds when the Z' mass is larger than the largest 
collider energy (about 209 GeV for LEP-II). Once again, the non-observation of any 
deviation from the SM expectations sets lower limits on the ratio Mz' I g'\ , rather than 
on the Z' mass alone. They will be summarised in section 3.1.1.2. 

3.1.1.1 Tevatron bounds on the Z'^_^ boson 

The Tevatron collider has been collecting data for several years. So far, no deviations in 
SM Drell-Yan production have been observed. Regarding the Z'^^^^ boson of this Thesis, 
the decay into pairs of electrons and into pairs of muons are considered. In general, 
the experimental communities release the data in the form of 95% C.L. excluded cross 
sections. Even though an explicit lower limit on the Z'g_^ mass is not shown, it is then 
easy to extract it from the data. We extract it by comparing the 95% C.L. excluded 
cross sections of Ref. [4; 5] with our theoretical prediction for the pp — t- Z'g_j^ — ?• 
e'^e~ cross sections. As mentioned in the introduction, the gauge coupling g[ is 
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a free parameter, so in all generality it is possible to identify the lower limit on the Z' 
mass per each g'^ fixed value. 

The latest available analyses are the D0 analysis of Ref. [4] using 5.4 fb~^ and the 
CDF analysis of Ref. [5] using 4.6 fb~^ of data, respectively, for electrons and muons in 
the final state. For the numerical evaluation of the Z' boson signal, we used CTEQ6L 
[50] as default PDFs, evaluated at the scale = M|^. The leading order (LO) cross 
sections are then multiplied by a mass independent fc— factor of 1.3 [45], as in Refs. [4; 5], 
to get in agreement with the Next-to-Next-to-Leading-Order (NNLO) QCD corrections. 
The Tevatron limits for the Z'q_^ boson are shown in table 3.1 (for selected masses and 
couplings) [39]. Notice that these are the most conservative limits, as they are evaluated 
for decoupled heavy neutrinos, i.e., with masses bigger than M^//2. 



pp 


—?■ e+e 


pp 




9[ 


Mz' (GeV) 


9[ 


Mz' (GeV) 


0.0197 


300 


0.0179 


300 


0.0193 


400 


0.0189 


400 


0.0281 


500 


0.0456 


500 


0.0351 


600 


0.0380 


600 


0.0587 


700 


0.0544 


700 


0.0880 


800 


0.0830 


800 


0.1350 


900 


0.1360 


900 


0.2411 


1000 


0.2220 


1000 


0.3880 


1100 


0.3380 


1100 



Table 3.1: Lower bounds on the Z' boson mass for selected g'l values in the B — L model, 
at 95% C.L., by comparing the collected data of Ref . [4; 5] with our theoretical prediction 
for pp ^ Z'q_j^ — )• e~^e~ {fj,^ n~) at Tevatron. 



3.1.1.2 LEP bounds on the Z'^_^ boson 

When the Z' mass is larger than the largest collider energy (about 209 GeV for LEP-II), 
an expansion in s/M"^, can be performed. This results in effective four-fermion contact 
interactions that have been bounded by LEP-II. Recall that the gauge coupling g'^ is 
a free parameter. One could suppress the Z' boson impact at lower energies also by 
reducing the coupling. Hence, an inverse law dependence from g'^^ is expected, and, in 
fact, the bound reads [45] 

— ^>6TeV, (3.1) 

9i 

or, in a more recent re-analysis, at 99% C.L., as [51] 

Mz^ 



9i 



> 7 TeV . (3.2) 



In this Thesis, we adopt the bound of eq. (3.2), that is the most conservative one. 

The LEP experiments are also able to provide a lower bound for the B — L breaking 
VEV X. In fact, making use of eq. (2.25), the LEP bound on the Z'g_^ boson mass of 
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eq. (3.2) can be rewritten as a lower bound for the B — L breaking VEV 

x>3.5TeV. (3.3) 

3.1.2 Theoretical constraints 

The RGE evolution gives us indications for the validity of the model concerning the gauge 
couplings. In particular, their evolution must stay perturbative up to some particular 
scale. In the B — L model, the conditions that the free parameters in the gauge sector 
must fulfil are ^ 

g[{Q')<l yQ'<Q and giQEw) = 0, (3.4) 

where the second condition in eq. (3.4) defines the pure B — L model. 

Varying the scale Q, the maximum scale up to which we want the model to be well- 
defined, we get an upper bound on g'i{QEw) as a function of Q, as shown in figure 3.1. 
Typical results are summarised in table 3.2. 

. r. 0.9 

0.8 
0.7 
0.6 
0.5 
0.4 
0.3 

4 6 8 10 12 14 16 18 

LogJQ/GeV) 

Figure 3.1: Maximum allowed values by eq. (3.4) for g'liQEw) the B — L model as a 
function of the scale Q. 




^Notice that the triviahty condition, as known in the hterature, reads gi{Q) < k, where fc = 1 or 
y/4iv. As shown in Ref. [52], this 'ad-hoc' prescription could be systematically improved by combining 
the unitarity and RGE evolution techniques. 
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Lo5io(Q/GeV) 


3 


5 


7 


10 


15 


19 


g'liQEw) 


0.860 


0.693 


0.593 


0.497 


0.397 


0.342 



Table 3.2: Maximum allowed values by eq. (3.4) for g'liQEw) the B — L model for 
selected values of the scale Q. 

3.2 Phenomenology of the gauge sector 

In this section we present the analysis of the new state in the gauge sector, the Z'^_j^ 
boson (or simply Z' boson, when the meaning is clear). 

An important feature of this Z' boson is the chiral structure of its couplings to 
fermions: since the B — L charge does not distinguish between left-handed and right- 
handed fermions (as clear in table 2.1), the B — L neutral current is purely vector-like, 
with a vanishing axial part ^. In fact. 



9z' 



9z' - 9z' 
2 



0. 



(3.5) 



As a consequence, we decided to not study the asymmetries of the decay products 
stemming from the Z'^_j^ boson, given their trivial distribution at the peak. 

The Z'^_^ boson is not always considered as a traditional benchmark for generic 
collider reach studies [46; 53; 54; 55; 56; 57; 58; 59] or data analyses [4; 5]. Main 
differences with the most popular Z' models in the literature (aside the vanishing axial 
coupling) are the fact that the gauge coupling g[ is a free parameter (as we do not 
consider here constraints from possible embeddings into GUTs) and the presence of new 
coupled matter (the f^.'s). Since their mass is as well a free parameter, they will affect 
the Z'^_^ boson BRs, which will be a function of the neutrino masses. Finally, in the 
pure B — L model, no mixing between the Z'^_j^ boson and the SM Z boson is present. 

In this section we give a detailed description of the Z' properties. In section 3.2.1, 
the production cross sections and decay BRs are shown. A detailed investigation of the 
discovery potential for the LHC is presented in section 3.2.2, for the foreseen center-of- 
mass (CM) energies of 7 and 14 TeV and integrated luminosities up to 1 fb~^ and to 
100 fb~^, respectively. For the former CM energy, a comparison with Tevatron is also 
shown, for an integrated luminosity up to 10 fb~^. 

3.2.1 Production cross sections and decay properties 

The most efficient hadro-production process involving a Z' boson is the Drell-Yan mode 



Z', 



(3.6) 



^This is strictly true for Dirac fermions, like quarks and charged leptons. Regarding the neutrinos 
(both light and heavy), they are Majorana states and, hence, the coupling to the Z'g_j^ boson is just 
axial [34]. 
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where q is either a valence quark or a sea quark in the proton. At the parton level, 
the Z' production cross section for process (3.6) depends on two main parameters: Mz' 
and g[. In figure 3.2a we present the Z' boson hadro-production cross section at the 
LHC as a function of both Mz' and §[, in the ranges 0.5 TeV < Mz' < 5 TeV and 
0.1 < g[ < 0.5, respectively, while figure 3.2b presents the contour levels in the (M^/, 
g[) plane for cross sections of 4 pb, 0.3 pb, 50 fb and 5 fb. 




Ct) 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

M^, (TeV) 

(a) 



0.5 



0.4 



0.3 



0.2 



0.1 




<3=4pb 
<5=0.3pb 
c=50fb 
a=5fb 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
M^. (TeV) 

(b) 



Figure 3.2: Z' hadro-production cross sections at the LHC over the (Mz', g'l) plane: 
(3.2a) as a function of Mz' for various g[ values and (3.2b) in the form of contour lines 
for four fixed values of production rates. In the left frame, the vertical ticks indicate 
the region excluded experimentally (on the left from the ticks), in accordance with sec- 
tion 3.1.1. In the right frame, the black shaded area is the portion excluded by eq. (3.2) 
(LEP bounds). The red shaded area is the region excluded by Tevatron, in accordance 
with table 3.1. 



For a fixed value of the coupling, we can compare the production cross sections 
for different hadronic machines and CM energies. For g[ = 0.1, the cross sections at 
Tevatron and the LHC (for y/s = 7, 10 and 14 TeV) are shown in figure 3.3. Note 
that although at Tevatron the production cross section is smaller than at the LHC, the 
integrated luminosity we are considering here for the LHC at 7 TeV (i.e., 1 fb~^) is 
smaller than for Tevatron (i.e., 10 fb^^). 

3.2.1.1 Decay properties 

As discussed earlier, the extra U{1)b~l gauge group provides an additional neutral 
gauge boson, Z', with no mixing with the SM Z boson. Therefore, the Z' boson decays 
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-5 



10 ^ 

1 1.5 2 2.5 3 3.5 4 



Figure 3.3: Cross sections forpp{p) — )■ Z'^ j^ at Tevatron and at the LHC (for y/s = 7, 10 
and 14 TeV) for g[ = 0.1. 



only to fermions at the tree level and its width is given by the following expression: 



Mz' 
Utt 



Ml 



(3.7) 



where mj is the mass and Cj the number of colours of the fermion type / and v-^ = 
[B — L) ■ q'y is the coupling (see table 2.1). 

Clearly, the Z' boson BRs depend strongly on the heavy neutrino mass and figure 3.4 
shows how they change with fixed (although arbitrary) values of m^^, for the following 
three cases: a heavy neutrino (i) much lighter than, (ii) lighter than, and (iii) comparable 
in mass to the Z' boson, in the range 0.5 TeV < Mz' < 5 TeV, after summing over the 
generations. 

A feature of the current B — L model illustrated in the previous figure is that the Z' 
boson predominantly couples to leptons. In fact, after summing over the generations, 
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Figure 3.4: Z' boson BRs as a function of Adz' for several heavy neutrino masses: 
ruy^ = 50, 200 GeV, and 1 TeV, from left to right, respectively. A summation over 
all lepton/ neutrino flavours is implied throughout, whereas, in the case of quarks we 
distinguish between light flavours (q = d,u,s,c,b) and the top quark. 



k = 1...3, we roughly get for leptons and quarks: 

fe k 

In particular, BR{Z' — )• varies between 12.5% and 15.5% {i = e,/i), while 

BR{Z' — qq) varies between 4% and 5%. 

Not surprisingly, then, for a relatively light (with respect to the Z' gauge boson) 
heavy neutrino, the Z' BR into pairs of such particles is relatively high: ~ 18% (at 
most, again, after summing over the generations). As we will see in section 4.2.1, the 
pair production from Z' boson decays is an effective way to produce the heavy neutrinos 
at the LHC, also because of the relatively high BR. The importance of this channel is 
due to the fact that it allows for the B — L model discrimination: as already mentioned, 
asymmetries of the standard Z' decay modes are trivial at the peak, as true for a variety 
of other models. To remove this degeneracy, one should look for further consequences 
of our model. Among these, the decay into heavy neutrino pairs is certainly the most 
spectacular, allowing for unusual Z' decays, such as multi-lepton and/or multi-jet decay 
patterns (see, for instance, Refs. [34; 47; 60]). After introducing the heavy neutrino 
properties, the whole section 4.2.2 will be devoted to the study of this signature, with 
particular emphasis on one specific decay pattern of the heavy neutrino pair, the so- 
called tri-lepton decay mode. 

One should finally note that possible Z' decays into one light and one heavy neutrino 
are highly suppressed by the corresponding (heavy-light) neutrino mixing and thus they 
can safely be neglected. 

In figures 3.5a and 3.5b we present the total decay width of the Z' boson as a 
function of Mz' and g'^^, respectively (with the other parameters held fixed to three 
different values), assuming that the partial decay width into heavy neutrinos vanishes. 
Also, figure 3.5c presents the relative variation of the total width as a function of the Vh 
mass for three different values of Mz' and for g'^ =0.1. 
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Figure 3.5: Z' total width as a function of: (3.5a) Mz' (for fixed values of g'l), (3.5b) 
9i (for fixed values of Mz') and (3.5c) rrij^^ (for fixed values of Mz' and = 0.1). In 
plots (3.5a)[(3.5b)] the portion of the curves to the left[right] of the vertical ticks (when 
appearing) are experimentally excluded, in accordance with eq. (3.2) and table 3.1. Notice 
finally that in plot 3.5b the curve for Mz' = 0.5 TeV is shown just for sake of comparison, 
as it is all excluded. 

From the first two plots we see that the total width of a Z' gauge boson varies from 
a few to hundreds of GeV over a mass range of 0.5 TeV < Mz' < 5 TeV, depending on 
the value of g'l, while from the third plot one can gather the importance of taking into 
consideration the heavy neutrinos, since their relative contribution to the total Z' width 
can be as large as 25% (whenever this channel is open). 
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3.2.2 Discovery power at the LHC 

In this section we determine the discovery potential of the LHC considering several CM 
energies, 7 and 14 TeV, using the expected integrated luminosities. For ^/s = 7 TeV only, 
we also compare our results for the LHC to the expected ultimate reach at Tevatron. 

For sake of comparison, a complementary discovery power study at the LHC for 
the Z'^_j^ boson also appears, for CM energy of 7 TeV, in Ref. [61] and for 14 TeV in 
Ref. [62] (for the Z' — )• e'^e~ channel only). 



Figure 3.6: Feynman diagram for the Drell-Yan di-lepton production (i = e, fi). 

The process we are interested in is the di-lepton production, shown in figure 3.6. 
We define our signal as pp — )■ 7, Z, Z'^_^ — )• {£ = e, //), i.e., all possible sources 

together with their mutual interferences, and the background as — )• 7, Z — )• £~^£~ 
{I = e, /i), i.e., SM Drell-Yan production (including interference). No other sources of 
background, such as WW, ZZ, WZ or it have been taken into account. These can be 
suppressed or are insignificant [4; 59]. For both the signal and the background, we have 
assumed standard acceptance cuts (for both electrons and muons) at the LHC: 



and we apply the following requirements on the di-lepton invariant mass, Mu, depending 
on whether we are considering electrons or muons. We distinguish two different scenarios: 
an 'early' one (for ^/s = 7 TeV) and an 'improved' one (for t/s = 14 TeV), and, in 
computing the signal significances, we will select a window as large as either one width 
of the Z'q_j^ boson or twice the di-lepton mass resolution ^, whichever the largest. The 
half windows in the invariant mass distributions respectively read, for the 'early scenario': 




p^T > 10 GeV 



r]^\ < 2.5 {£ 



e, //) 



(3.8) 



electrons: \M^ 



ee 




(3.9) 



muons: |M^^— M^/j < 




(3.10) 



and for the 'improved scenario': 



electrons: IM^ 



ee 




(3.11) 



muons: \Mun — Mz'\ < 




(3.12) 



^We take the CMS di-electron and di-muon mass resolutions [59; 63] as typical for the LHC envi- 
ronment. ATLAS resolutions [64] do not differ substantially. 
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Our choice reflects the fact that what we will observe is in fact the convolution between 
the Gaussian detector resolution and the Breit-Wigner shape of the peak, and such 
convolution will be dominated by the largest of the two. Our approach is to take the 
convolution width exactly equal to the resolution width or to the peak width, whichever 
is largest ^, and to count all the events within this window. Finally, only 68% of signal 
events are considered: intrinsically, when the peak width is dominating, effectively (by 
rescaling the signal), otherwise. 

In figure 3.7 we compare the LHC resolutions for electrons for the two aforementioned 
scenarios (eqs. (3.9) and (3.11)) with F^'/^- It is clear that, whichever the Z'^^j^ mass, 
for a value of the coupling g'^ smaller than roughly 0.4, the peak will be dominated by 
the early experimental resolution, i.e., the half window will contain an amount of signal 
as big as the one produced with |M« — Mzi\ = Vz'/'i- The region of interest in the 
parameter space we are going to study almost always fulfils the condition g'^ < 0.4, as we 
will see from the plots in the following section. The muon resolution is much worse and 
in such a plot it would be an order of magnitude higher than the other curves. Hence, 
for this final state, the peak is always dominated by the experimental resolution, for the 
values of the gauge coupling we are considering. 



Figure 3.7: Comparison of the CMS electron resolution according to eqs. (3.9) and (3.11), 
with the two different constant terms as described in the text, compared to Tz'/2 for 
g[ = 0.23 and g[ = 0.44. 

In the next section we will compare the LHC and Tevatron discovery reach. In the 
derivation of the experimental constraints, we refer to the latest publications, being the 
D0 analysis of Ref. [4] for the electron case and the CDF analysis of Ref. [5] for the 
muon final state, discussed in section 3.1.1.1. Hence, we have considered the typical 

^In details, for resolutions below r/2, we take the convolution equal to the resolution width. For 
resolutions above SF, we take the convolution equal to the peak width. When the resolution G [F/2, 3F], 
the convolution is taken as a linear interpolation between the two regimes. 




500 



1000 1500 2000 2500 3000 3500 4000 



Mz.(GeV) 



30 



3. GAUGE SECTOR 



acceptance cuts (for electrons and muons) for the respective detector: 



> 25 GeV, 
pf^ > 18 GeV, 



rj^l < 1.1 
^1 < 1, 



(3.13) 
(3.14) 



and the following requirements on the di-lepton invariant mass, Mu, depending on 
whether we are considering electrons or muons ^: 



The selection of an invariant mass window centred at the Z' boson mass is comparable 
to the standard experimental analysis, as in Ref. [4] (the electron channel at D0), where 
signal and background are integrated from Adz' — WTz' (where Tz' is the Z' boson 
width obtained by rescaling the SM Z boson width by the ratio of the Z' to the Z boson 
mass) to infinity. Since the background (in proximity of the narrow resonance) can be 
reasonably thought as flat, while the signal is not, the procedure we propose enhances 
the signal more than the background and it is expected to be more sensitive than the 
aforementioned one. Ref. [5] applies a different strategy and figure 3.8a shows that our 
procedure is comparable to it, although less involved. A Bayesan approach is being used 
at the LHC [66], similar to the CDF case. Hence, we present our results for a comparison 
'a posteriori'. 

In our analysis we use a definition of signal significance a, as follows. In the region 
where the number of both signal (s) and background (b) events is 'large' (here taken to 
be bigger than 20), we use a definition of significance based on Gaussian statistics: 



Otherwise, in case of smaller statistics, we used the Bityukov algorithm [67], which 
basically uses the Poisson 'true' distribution instead of the 'approximate' Gaussian one. 

Finally, for the numerical evaluation of the cross sections, we use the same setup as 
in section 3.1.1.1. In particular, the same mass independent fc— factor of 1.3 is also used 
for the LHC (as in Ref. [66] 

A typical detector resolution has effectively been taken into account by our proce- 
dure, that consists in counting all the events that occur within the window (in invariant 
mass) previously described, and by rescaling to 68% the signal events when the peak 
is dominated by the experimental resolution. Nonetheless, our simulation does not ac- 
count for Initial State Radiation (ISR) effects. ISR can have two main sources: QED-like 

^We take the D0 di-electron [65] and the CDF di-muon [5] mass resolutions as a typical Tevatron 
environment, in accordance with the most up-to-date limits. 

^Notice that in Ref. [66] the fc-factor used was mass-dependent. Here we use the average value. 



electron: |M, 




ee 




a = s/Vh. 



(3.17) 



31 



ISR (i.e., photon emission), that has the effect of shifting the peak and of creating a 
tail towards smaller energy, and QCD-like ISR (i.e., gluon emission), that has similar 
effects and might also induce trigger issues in the intent of removing backgrounds (e.g., 
by cutting on final state jets). Although we are aware of such effects, we believe that 
their analysis goes beyond the scope of this work and it will be argument of future in- 
vestigations. Altogether, we are confident that, while particular aspects of our analysis 
may be sensitive to such effects, the general picture will not depend upon these sub- 
stantially. Also, the only background considered here was the irreducible SM Drell-Yan. 
Reducible backgrounds, photon-to-electron conversion, efficiencies in reconstructing elec- 
trons/muons, jets faking leptons etc., whose overall effect is to deplete the signal, were 
neglected (being tt the most important source, at the level of < 10%). However, for 
this analysis they are not quantitatively important [4; 59; 68]. The net effect of the 
factors above is usually regarded as an overall reduction of the total acceptance, being 
the lepton identification the most important source, about 80 90% per each lepton. 
We comment on this in the conclusion of the chapter. 

3.2.2.1 LHC at ^/s = 7 TeV and comparison w^ith Tevatron 

The first years of the LHC work will be at a CM energy of 7 TeV, where the total 
integrated luminosity is likely to be of the order of 1 fb~^. Figure 3.8 shows the discovery 
potential under these conditions, as well as the most recent limit from LEP (see eq. 3.2) 
and from Tevatron (as in section 3.1.1.1). 

In the same figure we also include for comparison the Tevatron discovery potential 
at the integrated luminosities used for the latest published analyses [4; 5] (5.4 fb~^ 
and 4.6 fb~^ for electrons and muons, respectively) as well as the expected reaches at 
£ = 10fb-\ 

Notice that the Tevatron excluded area are based on the actual data, while the 
dot-dashed 2a curves are theoretical curves. Thus, if from the one side theory cannot 
reproduce experiments, from the other side we are comparing two methods of extracting 
the results. As mentioned previously, figure 3.8 shows that the procedures used in exper- 
imental analyses for the electron channel [4; 69] are not quite optimisied for maximising 
the signal significance. The alternative analysis described in this work has the potential 
to improve sensitivities and can be easily developed even further. 

It is then clear that Tevatron will still be competitive with the LHC (for -y/i = 7 TeV 
CM energy), especially in the lower mass region where the LHC requires 1 fb^^ to be 
sensitive to the same couplings as Tevatron. The LHC will be able to probe, at 5cr level, 
the Z'^_^ boson for values of the coupling down to 3.7 — 5.2 • 10~^ (for electrons and 
muons, respectively), while Tevatron can be sensitive down to 4.2 • 10~^ with electrons. 
The kinematic reach of the two machines is different. The LHC for 1 fb^^ can discover 
the Z'q_^ boson up to masses of 1.20 — 1.25 TeV, while at Tevatron a 2>a evidence will 
be possible up to a value of the mass of 1 TeV in the electron channel, for a suitable 
choice of the coupling. As clear from figure 3.8b, the muon channel at Tevatron requires 
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more than 10 fb~^ to start probing (at 3a) points in the Mz' — g'\ plane allowed by the 
CDF constraints, as in table 3.1. This total integrate luminosity appears to be more 
than what can be collected, due to the announced shutdown by the end of the year 2011 
[70; 71]. 
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Figure 3.8: Significance contour levels plotted against g[ and Mz' at the LHC for -y/s = 7 
TeVfor 0.1 - 1 fb'^ and at Tevatron (y/s = 1.96 TeV) for (3.8a, electrons) 5.4- 10/6^^ 
and (3.8b, muons) 4.6 — 10/6"^ of integrated luminosity. The shaded areas correspond 
to the region of parameter space excluded experimentally, in accordance with eq. (3.2) 
(LEP bounds, in black) and table 3.1 (Tevatron bounds, in red). 

Figure 3.9 shows the integrated luminosity required for 3cr evidence and So" discovery 
as a function of the Z'^ ^^ boson mass for selected values of the coupling for both electron 
and muon final states at the LHC, and for the electron channel only at Tevatron. The 
muon channel at Tevatron requires more than 10 fb~^ to start probing the Z'^^j^ boson 
at 3o", and, hence, we do not present it. 

We now fix some values for the coupling ((7^ = 0.158, 0.1, 0.08 for the LHC analysis, 
g'^ = 0.1, 0.08 for Tevatron) and we see what luminosity is required for discovery at each 
machine. For g'^ = 0.1 the LHC requires 0.35 — 0.70 fb~^ to be sensitive at 5cj, with 
electrons and muons, respectively, while Tevatron requires 10 fb~^ with electrons. For 
the same value of the coupling, Tevatron can discover the Z'q_j^ boson up to Mz' = 840 
GeV, with electrons in the final state, with 12 fb^^ of data. The LHC can extend 
the Tevatron reach up to Mz' = 1.0(0.9) TeV for g[ = 0.1. Regarding g[ = 0.08, 
a discovery can be made chiefly with electrons, requiring 0.4(12) fb~^, for masses up 
to 900(780) GeV at the LHC (Tevatron). For muons, the LHC requires 0.85 fb~^ for 
masses up to 800 GeV. Both machines will be sensitive at 3cr with much less integrated 
luminosities, requiring roug hly 0.12 - 0.15(0.25 - 0.35) fb"i to probe the Z^„^ at the 
LHC for electrons(muons) in the final state, for (7^ = 0.1 — 0.08. At Tevatron, 5.5 fb~^ 
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Figure 3.9: Integrated luminosity required for observation at 3(T and Scr vs. Mz' for 
selected values of g'l at the LHC for ^/s = 7 TeV for (3.9a) electrons and (3.9c) muons 
and at Tevatron (^/s = 1.96 TeV) for (3.9b) electrons (the muon channel requires more 
than 10 fb^^ and is, hence, not shown). Only allowed combinations of masses and 
couplings are shown. 



are required for probing at 2>a both values of the coupUng, for electrons only. Finally, 
bigger values of the coupling, such as g'^ = 0.158, can be probed just at the LHC, that 
is sensitive at 3cj to masses up 1.45(1.35) TeV using electrons (muons) and at 5cr to 
masses up to 1.2(1.15) GeV, requiring 0.2 — 0.6(0.3 — 0.9) fb~^ at least, at 2>a — 5(t for 
electrons (muons) . 
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The 5(7 discovery potential for the LHC at ^/s = 7 TeV and for Tevatron are sum- 
marised in the table 3.3, for selected values of couplings and integrated luminosities. 



LHC 


pp — >■ e+e 




pp fi+fi 




L (fb-i) 


g[ = 0.08 


g[ = 0.1 


g[ = 0.158 


g[ = 0.08 


g[ = 0.1 


g[ = 0.158 


0.2 


820(-) 


925(-) 


IIOO(-) 


-(-) 


-(-) 


-(-) 


0.3 


900(-) 


1000(800) 


1200(-) 


-(-) 


850(-) 


IIOO(-) 


0.5 


1000(775) 


1100(875) 


1300(-) 


825(-) 


950(-) 


1200(-) 


1 


1130(900) 


1250(1000) 


1450(1200) 


950(800) 


1080(900) 


1360(1130) 


Tevatron 


pp — )• e^e 




pp fi+fi- 




L (fb-i) 


g[ = 0.08 


g[ = 0.1 


g[ = 0.158 


g[ = 0.08 


9'i = 0.1 


g[ = 0.158 


8 


825(-) 


895(-) 


-(-) 


-(-) 


-(-) 


-(-) 


10 


850(-) 


915(825) 


-(-) 


-(-) 


-(-) 


-(-) 


12 


870(775) 


930(830) 


-(-) 


-(-) 


-(-) 


-(-) 



Table 3.3: Maximum Z'^_j^ boson masses (in GeV) for a 3<t(5(t) discovery for selected 
g'l and integrated luminosities in the B — L model, both at the LHC (for ^/s = 7 TeV) 
and at Tevatron (for ^/s = 1.96 TeV). No numbers are quoted for already excluded 
configurations. 



3.2.2.1.1 Exclusion power 

If no evidence for a signal is found at the LHC at the considered energy and luminosity 
configuration, 95% C.L. exclusion limits can be derived: in this subsection we present 
exclusion plots for the early stage of the LHC CM energy (-y/s = 7 TeV). We also show 
the expected exclusions at Tevatron for £ = 10 fb~^. 

We start by looking at the 95% C.L. limits presented in figure 3.10 for Tevatron and 
for this stage of the LHC (for 10 fb~^ and 1 fb~^ of integrated luminosity, respectively). 

One can see that the different resolutions imply that the limits derived using electrons 
are always more stringent than those derived using muons in excluding the Z'^^j^ boson. 
As for the discovery reach, Tevatron is also competitive in setting limits, especially in 
the lower mass region. In particular, using electrons and in case of no evidence at 
Tevatron with 10 fb~^, the Z'q_j^ boson can be excluded for values of the coupling 
down to 0.03 (0.04 for muons) for Mz' = 600 GeV. For the LHC, to set the same 
exclusion limit for the same mass, 1 fb^^ of integrated luminosity is required, allowing to 
exclude g'^ > 0.02(0.35) using electron(muons) in the final state. For the same integrated 
luminosity, the LHC has much more scope in excluding a Z'^_j^, for Mz' > 1.0 TeV. 

For a coupling of 0.1, the Z'^_^ boson can be excluded up to 1.40(1.25) TeV at the 
LHC considering electrons(muons) for 1 fb~^, and up to 1.0(0.9) TeV at Tevatron for 
10 fb~^ of data. For g[ = 0.05, the LHC when looking at muons(electrons) will require 
400(100) pb-^ to start improving the current available limits, while with 100 pb ^ it 
can set limits on g'^ = 0.158, out of the reach of Tevatron. It will ultimately be able to 
exclude Z'^_j^ boson up to Mz' = 1-6 TeV for 1 fb^^ (both with electrons and muons). 
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Figure 3.10: (3.10a) Contour levels for 95% C.L. plotted against g'l and Mz' at the LHC 
for selected integrated luminosities and integrated luminosity required for observation at 
3a and 5a vs. Mz> for selected values of g'l (in which only the allowed combination of 
masses and couplings are shown), for (3.10b) the LHC at y/s = 7 TeV and (3.10c) at 
Tevatron (y/s = 1.96 TeV), for both electrons and muons. The shaded areas and the 
allowed {Mz'^g'i) shown are in accordance with eq. (3.2) (LEP bounds, in black) and 
table 3.1 (Tevatron bounds, in red for electrons and in green for muons). 

The 95% C.L. exclusions for the LHC at y/s = 7 TeV and at Tevatron are summarised 
in table 3.4, for selected values of couplings and integrated luminosities. 
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LHC 
L (fb'^) 
0.1 
0.2 
0.5 
1 

Tevatron 



9'i 



0.05 



pp - 
1[ 



e ' e 

: 0.1 



9i 



' - 0.158 
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1075 
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1375 



1100 
1250 
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pp 



e ' e 



pp — /i^/x 



g[ = 0.05 



9'i 



0.1 



g[ = 0.158 



700 
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820 
950 
1100 
1250 

pp — ^ fJ-^fJ-" 



1225 
1425 
1575 



£ (fb 



0.05 



9'i 



0.1 



g[ = 0.158 



g[ = 0.05 g[ = 0.1 g[ = 0.158 



10 
12 



750 
775 
800 
825 



950 
975 
1000 
1020 



680 



860 
875 
900 



Table 3.4: 



Maximum Z'^^ 



boson masses (in GeV) for a 95% C.L. exclusion for selected 



g[ and integrated luminosities in the B 
excluded configurations. 



L model. No numbers are quoted for already 



3.2.2.2 LHC at ^ = 14 TeV 

We consider here the design performance, i.e., ^/s = 14 TeV of CM energy with large 
luminosity, L = 100 fb~^. Figure 3.11 shows the discovery potential for the Z'q_j^ 
boson under these conditions, while figure 3.12 shows the integrated luminosity required 
for 3a evidence as well as for 5a discovery as a function of the ^^-L boson mass for 
selected values of the coupling at ^/s = 14 TeV. We consider the integrated luminosity 
in the range between 10 pb"^ up to 100 fb~^ After some years of data analysis, the 
performances of the detector will be better understood. We therefore use the resolutions 
for both electrons and muons quoted in eqs. (3.11) and (3.12), respectively. 

Prom figure 3.11, we can see that the LHC at y/s = 14 TeV will start probing a 
completely new region of the parameter space for L > 1 fb~^. For £ > 10 fb^^ a Z'q_j^ 
gauge boson can be discovered up to masses of 4 TeV and for couplings as small as 
0.01(0.02) if we are dealing with electrons (muons). At L = 100 fb~^, the coupling can 
be probed down to values of 8 • 10~^ in the electron channel, while couplings smaller 
than 1.6 • 10~^ cannot be accessed with muons. The mass region that can be covered 
extends towards 5 TeV irrespectively of the final state. 

As before, figure 3.12 shows the integrated luminosity required for 3{5)a evidence 
(discovery) of the Z'q_j^ boson as a function of the mass, for selected values of the 
coupling. We explore the range in luminosities, from 10 pb~^ to 100 fb"'^. However, 
just the configuration with g[ =0.1 can be probed with very low luminosity, requiring 
30(100) pb~^ and 50(150) pb^^ at 3o"(5cr) considering electrons and muons in the final 
state, respectively. For values of the coupling such as 0.05 and 0.2, 90(220) pb^^ and 
60(200) pb~^ are the integrated luminosities required to start to be sensitive (at 3{5)a) 
if electrons are considered, while 160(500) and 70(220) pb~^ are the least integrated 
luminosity required, respectively, if instead we look at muons. It is worth to emphasise 
here that the first couplings that will start to be probed at the LHC are those around 
g[ = 0.1. 
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Figure 3.11: Significance contour levels plotted against g[ and Mz' at the LHC for 
y/s = 14 TeV for several integrated luminosities for (3.11a) electrons and (3.11b) 
muons. The shaded areas correspond to the region of parameter space excluded experi- 
mentally, in accordance with eq. (3.2) (LEP hounds, in black) and table 3.1 (Tevatron 
bounds, in red). 




(a) (b) 
Figure 3.12: Integrated luminosity required for observation at 3a and 5a vs. Mz' for 
selected values of g[ at the LHC for y/s = 14 TeV for (3.12a) electrons and (3.12b) 
muons. Only allowed combination of masses and couplings are shown, in accordance 
with eq. (3.2) and table 3.1. 



The better resolution in the case of electrons reflects in a better sensitivity to smaher 
Z'^_i^ masses with respect to muons. For Mz' = 600 GeV, the LHC with = 14 TeV 
requires 1.0 fb~^ to be sensitive at ba to a value of the coupling of 0.025 in the electron 
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channel. If we are considering muons, 3.5 fb^^ is the required luminosity to probe at 5a 
the same value of the coupling. 

The 5(7 discovery potential for the LHC at ^/s = 14 TeV is summarised in table 3.5, 
for selected values of Z'q_^ masses and couplings. 



= 14 TeV 


pp — > e+e 


pp n+i^t 


si 


Mz' = 1 TeV 


Mz' = 2 TeV 


Mz' = 3 TeV 


Mz' = 1 TeV 


Mz' = 2 TeV 


Mz' = 3 TeV 


0.025 


2.5(7.0) 


50(>100) 


>100(>300) 


15(30) 


>100(>100) 


>300(>300) 


0.05 


0.4(1.0) 


9(20) 


80(>100) 


0.8(2.5) 


20(50) 


>100(>100) 


0.1 


0.07(0.2) 


1.5(4.0) 


15(50) 


0.1(0.3) 


2.0(6.0) 


20(65) 


0.2 


-(-) 


0.3(1.0) 


3(10) 


-(-) 


0.4(1.2) 


3(12) 



Table 3.5: Minimum integrated luminosities (in fb^^) for a 3a (5a) discovery for selected 
Z'^_j^ boson masses and g[ couplings for the B — L model. No numbers are quoted for 
already excluded configurations. 

Figures 3.13 and 3.14 show a pictorial representation of the Z' properties and line- 
shapes (widths and cross sections) for selected benchmark points on the 3a and 5a lines 
for 10 fb~^ of data at ^/s = 14 TeV, plotting the di-lepton invariant mass to which 
just the cuts of eq. (3.8) have been applied (without selecting any mass window). The 
binning is equal to the typical resolution in each case. 
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7, Z, Z'^_j^ — )• e'^e ) for several masses and couplings 



(Mz'/TeV, g[, Tz'/GeV): (3.13a) ("0.6, 0.0075, 0.006^, ^1.1, 0.015, O.OSj, (l.%, 0.025, 
0.2i;, (^2.0, 0.04, 0.67;, (2.Q, 0.07, 2.7^ and (3.Q, 0.2, 31;,- (3.13b) (Q.Q, 0.009, O.OOQj, 
(l.l, 0.02, 0.09;, (l.Q, 0.04, 0.53;, C2.1, 0.07, 2.2;, ("2.6, 0.12, 7.9; and (3.A, 0.3, 6i;, 
from the 3a and 5a lines at 10 fb~^ of figure 3.11 (y/s = 14 TeV), respectively, using 
15 GeV binning. Notice that the asymmetry of the peaks is the result of our choice to 
consider the full interference structure. 



The improved resolution for electrons allows a measure of the Z'^ ^^ boson width not 
only at high masses, but also opens the possibility of a measurement even for smaller 
masses. 
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(a) (b) 
Figure 3.14: -^-{jpp — )• ^,Z.,Z' — )• /U^/.f~) for some masses and couplings (Mz' / TeV, 
g'^, Tz'/GeV): (3.14a) (2.6, 0.07, 2.7) and (3.4, 0.2, 29) and (3.14b) (2.6, 0.12, 8) and 
(3.4:, 0.3, 65), from the 3a and 5a lines at 10 fb~^ of figure 3.11 (y/s = 14 TeV), using 
125 GeV binning. Notice that the asymmetry of the peaks is the result of our choice to 
consider the full interference structure. 

3.2.2.2.1 Exclusion power 

If no evidence for a signal is found at the LHC either at its designed energy and lumi- 
nosity configuration, very strong 95% C.L. exclusion limits can be derived. 

Due to the improved resolutions for both electrons and muons, they have very similar 
exclusion powers for couplings g'l > 0.1, therefore setting similar constraints. Depending 
on the amount of data that will be collected, several maximum bounds can be set (see 
figure 3.15a): i.e., with 10 fb~^ of data, the LHC at 14 TeV can exclude at 95% C.L. 
up to a mass of roughly 5 TeV for a value of the coupling g'^^ = 0.5, that we arbitrary 
take as the biggest allowed value for the consistency of the model from a RGE analysis 
of its gauge sector, as in section 3.1.2. For 100 fb~^ and for the same value of the 
coupling, the LHC can exclude at 95% C.L. masses up to roughly 6 TeV. For 10 fb~^ it 
will be possible to exclude a Z'q_^ boson for Mz' = 600 GeV if the coupling is greater 
than 1.8 • 10^^(9 • 10"'^) for muons (electrons), and values of the coupling greater than 
1.5 • 10~^(6.5 • 10"'^) for an integrated luminosity of 100 fb^^. Figure 3.15b shows the 
integrated luminosity that is required to excluded a certain Z'^ j^ boson mass for fixed 
values of the coupling. As previously noticed, electrons and muons set the same limits 
for g'-^ > 0.1. An integrated luminosity of 10 fb^^ is required to exclude a Z'^_^ mass 
up to 3.8 TeV for g[ = 0.2, instead 40 fb~^ reduces this to g[ = 0.1. For an integrated 
luminosity of 10 fb~^ the LHC experiments will be able to exclude the Z'^_^ boson for 
masses up to 3.1(3.1) TeV for g[ = 0.1, 2.5(2.1) TeV for g[ = 0.05 and 1.7(1.1) TeV 
for g[ = 0.025, when considering decays into electrons(muons). With 100 fb~^ of data. 
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Mz(TeV) Mz,(TeV) 

(a) (b) 
Figure 3.15: (3.15a) Contour levels for 95% CL. plotted against g'l and Mz' at the 
LHC for selected integrated luminosities and (3.15b) integrated luminosity required for 
observation at 3a and 5ct vs. Mz' for selected values of g'l (in which only the allowed 
combination of masses and couplings are shown), for ^/s = 14 TeV, for both electrons 
and muons. The shaded areas and the allowed {Mz',g'i) shown are in accordance with 
eq. (3.2) (LEP bounds, in black) and table 3.1 (Tevatron bounds, in red for electrons 
and in green for muons). 



more stringent bounds can be derived: for g'^ = 0.05(0.025) the Z'^_j^ boson can be 
excluded for masses up to 3.6(2.6) TeV in the electron channel, and up to 3.1(1.9) TeV 
in the muon channel. 

The 95% C.L. exclusions for the LHC at y/s = 14 TeV are summarised in table 3.6, 
for selected values of Z' masses and couplings. 
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Table 3.6: Minimum integrated luminosities (in fb^^) for a 95% C.L. exclusion for 
selected Z'^_j^ boson masses and g'l couplings in the B — L model. No numbers are 
quoted for already excluded configurations. 



3.3 Conclusion 

In this chapter, we have presented the results of our investigation of the gauge sector. 

First, we have presented the experimental constraints on the free parameter in this 
sector, the Z' boson mass and gauge coupling g'^. We have shown that the LEP indirect 
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limits and the Tevatron direct ones are complementary, being the former(latter) con- 
straining more for heavy(hght) Z' boson. We have also shown that a constraint on g'^^ 
can be derived from a RGE analysis of the running gauge couplings, whose equations 
decouple from the others. The RGE analysis gives an upper bound on 5^ at the EW 
scale, to avoid Landau poles up to a certain cut-off energy scale. 

We have then moved on, to study the phenomenology of the Z' boson in the pure 
B — L model. We have presented production cross sections, decay widths and BRs for 
2-body decay modes. On the parameter space that we have considered, the Z' boson is 
a rather narrow resonance, broadening up for high masses and big couplings. The decay 
into heavy neutrinos is certainly interesting, with a total BR up to ~ 20%, at most. We 
will study it with great details in the next chapter. 

We have shown that the Z' boson is dominantly coupled to leptons (their total BR 
sums up to 3/4). It only holds vectorial coupling to Dirac fermions, while neutrinos, 
being Majorana particles, couple to the Z' boson with pure axial couplings. 

Finally, we have presented the discovery potential for the Z'^_j^ gauge boson at the 
LHC for CM energies of -y/s = 7 and 14 TeV, using the integrated luminosity expected 
at each stage. This has been done for both the Z'^_^ — )• e+e" and Z'^_^ — fi'^fj.~ decay 
modes, and includes the most up-to-date constraints coming from LEP and Tevatron. 
We have proposed an alternative analysis that has the potential to improve experimental 
sensitivities. We also haved looked in detail at the different resolutions, showing that 
electrons and muons present very similar discovery power for values of the coupling 
bigger than roughly 0.1, for ^/s = 14 TeV. 

We are overall confident that the inclusion of further background, as well as a realistic 
detector simulation, will not have a considerable impact on the results we have presented. 
In fact, as noted in section 3.2.2, all detector effects can be casted in the form of a signal 
acceptance, including also the effect of kinematic and angular acceptance cuts. By 
looking at Refs. [4; 5; 66], we estimate an overall acceptance factor of ~ 70%, which 
we found to be approximately constant over the mass regions considered, to be applied 
to our parton level results [once the cuts of eqs. (3.8) and (3.13) are considered]. This 
acceptance is mainly related to the lepton identification, both at Tevatron and at the 
LHC. On the significance, as in eq. (3.17), the reduction is then of ~ 84%. 

A general feature is that greater sensitivity to the Z'^_j^ resonance is provided by 
the electron channel. At the LHC this has better energy resolution than the muon 
channel. A further consequence of the better resolution of electrons is that an estimate 
of the gauge boson width would eventually be possible for smaller values of the ^^-L 
boson mass than in the muon channel. The simultaneous measure of the Z' boson mass 
and gauge coupling (for instance, through a line-shape fit) will also enable to infer the 
B — L breaking VEV, given the simple formula that relates them. Limits from existing 
data imply that the first couplings that will start to be probed at the LHC are those 
around g'l = 0.1. Increased luminosity will enable both larger and smaller couplings to 
be probed. 
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Our comparison has shown that, for an integrated luminosity of 10 fb~ , Tevatron 
is stih competitive with the LHC in the electron channel and in the small mass region, 
being able to probe the coupling at the level of 5a down to a value of 4.2 • 10^^. The 
LHC will start to be competitive in such a region only for integrated luminosities close 
to 1 fb-i at ^ = 7 TeV. Also, at = 7 TeV, the mass reach will be extended 
from the Tevatron value of Mz' = 850 GeV, with electrons, up to 1.25(1.20) TeV for 
electrons (muons). The muon channel at Tevatron needs more than 10 fb~^ to start 
probing the Z'q_^ at 3a. Hence, it has not been studied. 

When the data from the high energy runs at the LHC becomes available, the disco- 
very reach of Z'q_j^ boson will be extended towards very high masses and small couplings 
in regions of parameter space well beyond the reach of Tevatron and comparable in scope 
with those accessible at a future LC [72]. 

If no evidence is found at any energies, 95% C.L. limits can be derived, and, given 
their better resolution, the bounds from electrons will be more stringent than those from 
muons, especially at smaller masses. 



Chapter 4 

Fermion sector 



In this chapter we study the phenomenology at the LHC of the new particles in the 
fermion sector of the pure B — L model, namely, the heavy neutrinos. 

We start in section 4.1 by presenting the existent constraints on the fermion sector. 
Special mention is for the light neutrinos masses, for which no absolute measure exists, 
with the strongest bounds coming from cosmological observables. 

In section 4.2 the properties of the heavy neutrinos are presented, i.e., their produc- 
tion cross sections, widths, and BRs. The possibility for these to be long-lived particles 
is rather interesting, allowing for signatures with distinctive displaced vertices, in which 
a high energetic and isolated pair of leptons point to a different vertex than the primary 
one. Bringing pieces together, the process — )■ Z' — )• Uhi^h is studied. In particular, 
we show that a suitable and peculiar signature to investigate is the tri-lepton decay 
of the Z' boson via heavy neutrino pairs. This signature provides a powerful insight 
in the fermion sector of the model (and in its interplay with the gauge sector), and it 
allows for the measure of the heavy neutrino masses. It is remarkable that with the 
simultaneous measure of the neutrino lifetime (through the displacement of its decay 
vertex) and mass, light neutrino masses can be inferred. For two benchmark points, 
in section 4.2.3 we present a detailed study of the tri-lepton signature (the main result 
of this Thesis), where signal and backgrounds are fully simulated and analysed (at the 
parton level and at the detector level), showing that it yields a good signal over back- 
ground ratio that makes it observable at the LHC, especially as rather generic cuts for 
background rejection are considered. 

Finally, we draw the conclusions for this chapter in section 4.3. 

The work in this section is published, in part, in [47]. 

4.1 Experimental constraints 

LEP-I data [49] have established the existence of exactly 3 neutrinos that couple to the 
SM Z boson, with masses below Mz/2. Therefore, to avoid this bound, we require the 
heavy neutrinos to be heavier than Mz/2 ~ 46 GeV. 
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For the SM (or light) neutrinos, no absolute mass measure exists (see [73]). Nonethe- 
less, oscillation experiments provide a measure for squared mass differences [8; 9]: 

Sm,2^ ~ 7.67t°:}6 . 10-5 eY\ (4.1) 

(5m23^ ~ 2.39^0.11 ' 10"^ eV^ . (4.2) 

and mixing angles [8; 9]: 

(sin.?i2)' = 0.312t^,f,l, (4.3) 

(sin ^23)' = (4.4) 

{sin'disf = 0.016 ±0.010. (4.5) 

A cosmological upper bound of ''^^m^i^ < 0.58 eV also exists [10]. Ultimately, light 

I 

neutrinos have been taken to be degenerate and with a mass of ml,^ = 10^^ eV. 



4.2 Phenomenology of the heavy neutrinos 

After the diagonalisation of the neutrino mass matrix realising the see-saw mechanism, 
as in section 2.2.3, we obtain three very light neutrinos (f/), which are the SM-like 
neutrinos, and three heavy neutrinos (j^/i). The latter have an extremely small mixing 
with the I'l^s thereby providing very small but non- vanishing couplings to gauge and 
Higgs bosons. 

As already stated, we are interested here in delineating their impact onto the collider 
phenomenology, especially concerning the interaction with the gauge sector. Therefore, 
the analysis of the signatures from heavy neutrinos is mainly described in the light of 
discovery at hadronic colliders, in general, and at the LHC, in particular. In this section, 
we consider as discovery channel the pair production via the Z' boson. 

In this respect, to maximise and highlight the interesting patterns that can be ob- 
served, heavy neutrinos will be taken relatively light (with respect to Mz'), and, for 
simplicity, degenerate in mass. As we pointed out in the introduction, the latter is an 
approximation for illustrative though realistic purposes. 

The distinctive features of the B — L model take place because the heavy neutrinos 
decay predominantly to SM gauge bosons, in association with a lepton (either charged 
or neutral, depending on the electrical nature of the SM gauge boson). Being these 
couplings see-saw suppressed, in a large portion of the parameter space heavy neutrinos 
can be long-lived particles. Also, once heavy neutrinos are pair-produced via the Z' 
boson, they give rise to novel and spectacular multi- lepton and/or multi-jet decay modes 
of the intermediate boson. In the following sections are shown the consequences for 
collider searches. 

It is reasonable to conclude that the heavy neutrinos truly carry the hallmarks of 
the B — L model at colliders. 
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4.2.1 Production cross sections and decay properties 

We focus here in the pair production via Z' gauge boson exchange in the s-channel, 
whose Feynman diagram is in figure 4.1. Figure 4.2 shows the cross sections for the 
process 

J^i^K (4.6) 

i 

at the LHC, for yfs = 7 and 14 TeV CM energies, in the M^' — g'\ plane, for several 
heavy neutrino masses: mj^^ = 100 GeV (straight line), 200 GeV (dashed line), and 400 
GeV (dotted hne). 




Figure 4.1: Feynman diagram for heavy neutrino pair production. 
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Figure 4.2: Heavy neutrino pair production cross sections at the LHC for (4.2a) y/s = 7 
TeV and for (4-2b) y/s = 14 TeV, as a function ofm^y^, for several Z' masses and g'^ cou- 
plings. The shaded areas shown are experimentally excluded in accordance with eq. (3.2) 
(LEP hounds, in black) and table 3.1 (Tevatron bounds, in red). Heavy neutrinos have 
been summed over the generations. 

Complementary to this, figure 4.3 shows the production cross sections for the process 
of eq. (4.6) as a function of the heavy neutrino mass, for a choice of A/^/ and g'l, at the 
LHC, again comparing the two foreseen CM energies. A similar plot, for -y/i = 14 TeV 
only, can be found in Ref. [60]. 
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Figure 4.3: Heavy neutrino pair production cross sections at the LHC (for y/s = 7 and 
14 TeV) for several Z' masses and g'l couplings. Heavy neutrinos have been summed 
over the generations. 

It seems clear that the process of eq. (4.6) can hardly be tested at the LHC in its 
early stage (i.e., for y/s = 7 TeV CM energy), where the total integrated luminosity 
will be not above 1 fb~^. Even supposing the maximum value for the cross section, of 
~ 50 fb for Mz' ~ 900 GeV, g[ ~ 0.13 and m,^^ ~ 100 GeV, only around 50 events are 
produced, and considerably wiped out once detector geometrical and kinematical cuts 
are included, even before selecting any particular decay mode. 

On the contrary, this mechanism is suitable for testing at the LHC for its design 
performances (i.e., for y/s = 14 TeV CM energy), where, depending on the Z' boson 
mass and the value of the g[ coupling, there exist configurations in the parameter space 
allowing for thousands (or tens of thousands) of heavy neutrino events. In section 4.2.3, 
a detailed analysis of a particular decay mode of the neutrino pair will be shown, namely, 
the tri-lepton mode, for Mz' = 1.5 TeV and g[ = 0.2 (green solid curve of figure 4.3). 
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4.2.1.1 Decay properties 

As clear from eqs. (2.37) and (2.40), heavy neutrinos are predominantly RH states, 
with a tiny [O(10~^)] component of LH helicity, providing very small but non- vanishing 
couplings to gauge and Higgs bosons. In turn, it is this tiny component that enables the 
following Uh decays, when kinematically allowed: 

Vh ^ t^W^, (4.7) 

Vh ^ 1^1 (4.8) 

Vh nhi, (4.9) 

T^h vih2, (4.10) 

Uh ^ viZ'. (4.11) 
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Figure 4.4: Heavy neutrino BRs versus its mass for the fixed Mz' = 1.5 TeV, rrih^ = 150 
GeV, rrih^ = 450 GeV, a ~ 10~^ rads and g[ = 0.2. Here, W means the sum over 
and W~ . 



Figure 4.4 presents the corresponding BRs versus the heavy neutrino mass for the 
values of the other relevant B — L parameters given in the caption. One can see that 
the BR(i^/i — )• I'^W^) is dominant and reaches the 1/2 level in the my^^ » M^y, Mz,mh^ 
limit, while BIi{i'h i^lZ) and BR(z^/i — )• i^ihi) both reach the 1/4 level in this regime 
(for a — )• 0). Schematically, as in Ref. [34], 

BRii^h e^W) w BR{uh rW+) w BR{uh uiZ) , (4.12) 
w BR{uh^nhi)+BR{uh^vih2). (4.13) 
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In contrast, the Uh — )• I'lZ' decay channel (when kinematically open) is well below 
the percent level and is negligible for our study. To simplify the discussion in this 
section, we assume that the heavy neutrino masses are smaller than both Higgs boson 
masses. Under this assumption z^/, — >• uihi [i = 1,2) is not kinematically possible and 
BR(i/fe ^'^W^)[B'R{vh viZ)] reaches the 2/3[l/3] level in the m^^ > Mw, Mz limit. 

Nonetheless, the possibility of producing the light Higgs boson from a heavy neutrino 
decay is quite interesting and peculiar, as in once it involves new particles from all the 
sectors of the B — L model. In section 5.2, we will show the cross section for this process, 
highlighting its importance for a future LC. 

We recall that the heavy neutrino couplings to the weak gauge bosons are propor- 
tional to the (squared root of the) ratio of light to heavy neutrino masses (see the 
Appendix A), which is extremely small [O(10~^) for rriyf^ ~ 100 GeV]. Hence, the decay 
width of the heavy neutrino is correspondingly small and its lifetime large. The heavy 
neutrino can therefore be a long-lived particle and, over a large portion of parameter 
space, its lifetime can be comparable to or exceed that of the b quark. (In fact, for 
m^, = 10~^ eV and = 200 GeV they are equal). 

In figure 4.5 we present the heavy neutrino lifetime (top panel) in picoseconds and 
the proper decay length (or mean path) (bottom panel) in micrometers as a function of 
the light neutrino mass for three different choices of heavy neutrino masses (m^^ = 100, 
150 and 200 GeV). The proper decay length is defined as ctq, where tq is the lifetime 
of the heavy neutrino. The purple (horizontal) band presents the proper decay length 
of the b quark [73], while the blue band indicates the range of a typical microvertex 
detector [63] . The red band shows the region of light neutrino masses excluded by direct 
measurements of neutrino oscillations [see eq. (4.1)], by taking the lightest neutrino to 
be massless (thus the other neutrinos cannot populate this region, only the lightest one 
if massive). One should also note that the lifetime and the proper decay length of the 
heavy neutrinos in the laboratory frame will actually be equal to those given in figure 4.5 
times the Lorentz factor equal to py^/m^^ defined by the ratio Mz'/my^ which can be 
as large as about a factor of 10. We can then see that there exists a region where the 
heavy neutrino lifetime is of the same order as that of the b quark (shown as a purple 
band). The mean path and the respective lifetime of heavy neutrinos can therefore be 
measured from a displaced vertex inside the detector. The heavy neutrino can however 
be distinguished from a b hadron through the observation of vertices consisting of only 
two isolated leptons. (Such a SM B meson decay, while possible, would have a very 
small BR < 10"^ [73].) 

An experimentally resolvable non-zero lifetime along with a mass determination for 
the heavy neutrino also enables a determination of the light neutrino mass. The lifetime 
measurement allows the small heavy-light neutrino mixing to be determined and, as one 
can see from eqs. (2.34)-(2.37), this, along with the heavy neutrino mass, gives the light 
neutrino mass. Considering only one generation for simplicity, this is expanded upon 
below. 
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Figure 4.5: Heavy neutrino lifetime (top panel) and proper decay length (or mean path) 
CTo (bottom panel) as a function of the light neutrino mass for three different choices of 
heavy neutrino masses. The purple (horizontal) band presents the proper decay length 
of the b quark while the blue band indicates the range of a typical microvertex detector. 
The red band shows the region excluded by neutrino oscillation direct measurements. 




Neutrino mass eigenstates are related to gauge eigenstates by eq. (2.37), hence the 
eigenvalues are given by solving the equation 

which yields 

m^j = sin 2au ruD + sin^ M , (4-14) 
ruiyf^ = — sin 2a,y m£) + cos^ M , (4-15) 

with given by eq. (2.36). We have then three parameters (m/j, M, and a^) and a 
constraint [given by eq. (2.36)], that can be used to eliminate one parameter from the 
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above equations. 

The Feynman rules given in section A. 1.2 demonstrate that heavy neutrino interac- 
tions are determined entirely by the (sine of the) mixing angle a,^, as is the total width 
(and therefore the mean decay length). Hence, it is convenient to keep m^i^ and as 
independent model parameters eliminating niD from eq. (2.36), 

ruD = mD{au,muJ . (4.16) 

By measuring the heavy neutrino mass we can also invert eq. (4.15) 

M = M{au,m^J , (4.17) 

to finally get a fully known expression for the SM light neutrino mass as a function of 
our input parameters rrii,^ and (that can be measured independently), by inserting 
eqs. (4.16)-(4.17) into eq. (4.14), 

m^^{mD,M) = mui{ai,,muj . (4.18) 

This simple picture shows that within the B — L model we have an indirect way of 
accessing the SM light neutrino mass by measuring the mass of the heavy neutrino 
and the kinematic features of its displaced vertex. If the whole structure of mixing 
is taken into account, including inter-generational mixing in the neutrino sector (as in 
Ref. [34]), the task of determining the light neutrino masses in this way would become 
more complicated but the qualitative features and the overall strategy would remain the 
same, thereby providing one with a unique link between very large and very small scale 
physics. This in turn will also be a direct test of the specific realisation of the see-saw 
mechanism within the gauged B — L framework, as described in section 2.2.3. 

In the next subsection we will combine the heavy neutrino decay patterns with the 
pair production mechanism via the Z' boson. The overall process can be seen as non- 
standard decay patterns of the Z' boson. We will analyse their BRs by distinguishing 
the number of charged leptons in the final state, individuating the most suitable one for 
a detailed study, aimed to demonstrating that there exists a signature that enables the 
heavy neutrino mass to be measured at the LHC. 

4.2.2 Phenomenology of the neutrino pair production 

The possibility of decays of the Z' gauge boson into pairs of heavy neutrinos is one 
of the most significant results of this work, since, in addition to the clean SM-like di- 
lepton signature, it provides multi-lepton signatures where backgrounds can be strongly 
suppressed. In order to address this quantitatively, we first determine the relevant BRs. 

As in eqs. (4.7)-(4.8), a single heavy neutrino decay will produce a signature of 0, 1, 
or 2 charged leptons, depending on whether the heavy neutrino decays via a charged or 
neutral current and on the subsequent decays of the SM and Z gauge bosons. We 
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can have both chains 



and 



Vh^i^W^ ^£^ + \ (4.19) 

hadrons 



Vh^viZ^vi + \ . (4.20) 

I z^; z^;/ hadrons 

The pattern in eq. (4.19) provides 1 or 2 charged leptons, while that in eq. (4.20), or 
2. Hence, multi-lepton signatures may arise when the Z' gauge boson decays into a pair 
of heavy neutrinos, producing up to 4 charged leptons in the final state. There exist in 
total 5 different topologies, distinguished by the number of charged leptons in the final 
state, that are: 



lepton mode: Z' 



6ui via 2Z bosons 

4z// 2j via 2Z bosons (4-21) 
2z/; Aj via 2Z bosons 



f £^ Zvi 2j YiaW + Z bosons 
1 lepton mode: < , (4.22 



2 lepton mode: Z' 



^ ui 4j via W + Z bosons 

£^ Aui via 2Z oiW + Z bosons 

£^ £=F 2vi 2j via 2Z or W + Z bosons (4.23) 
£± £=F 4j via 2W bosons 

3 lepton mode: Z' 21^ 2j vi via 2W oi W + Z bosons (4.24) 

4 lepton mode: Z' 2£^ 2t^ 2ui via 2W , 2Z oi W + Z bosons (4.25) 

Eq. (4.21) shows that the Z' boson could decay completely invisibly, or in a multi-jet 
final state with large missing energy. As well, a mono-leptonic decay is possible, as in 
eq. (4.22), with multi-jets and, again, possibly large missing energy. Although rather 
interesting (especially the completely invisible decay), they are usually overwhelmed 
by the backgrounds. Therefore, we consider only multi-lepton decay modes, i.e., with 
2 or more charged leptons in the final state. Figure 4.6 shows the BRs of a Z' boson 
decaying into 2 (top-left panel) and 3 or 4 (top-right panel) leptons (plus possibly missing 
transverse momentum and/or jets, as appropriate) as a function of ^j^^, where a lepton 
can be either an electron or a muon and these contributions are summed. While the 
former are clearly dominant, the latter are not at all negligible. For M^y± < m^j^^ < Mz, 
the Vh — ^ £^ decay is the only one kinematically possible, whereas for m,y^ < M^y± 
the heavy neutrino can decay only via an off shell W/Z and can therefore be very long- 
lived (depending on rrii^^). For a very massive Z' boson (2 TeV < Mz' < 5 TeV) the 
multi-leptonic BRs are roughly 2.5% in the case of Z' — ?• 3£ and 0.5% in the case of 
Z' —7- 4£, for a wide range of heavy neutrino masses. 

Finally, as in eq. (4.23), also non-standard di-lepton decays are possible, with the Z' 
boson that decays into 2 leptons plus a large amount of missing transverse momentum 
and/or highly energetic jets, see figure 4.6 (bottom- left panel). Particularly interesting 
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Figure 4.6: Z' boson BRs, as a function ofm^i^ , into: 2 leptons (both e and fi, top-left 
panel); 3 and 4 leptons + X (both e and fi, top-right panel); 2 leptons -h X jets (both 
e and /i, bottom-left panel); zoom of the previous plot with same legend (bottom-right 
panel). 



is the decay into 2 leptons and 4 jets (only via W bosons), since here there is no missing 
transverse momentum at all and its BR is rather large with respect to the other non-SM 
signatures, as we can see in figure 4.6 (bottom-left panel). 

4.2.3 Tri-lepton signature analysis (parton level vs. detector level) 

The non-standard decay modes of the Z' boson (via the decay into heavy neutrino pairs) 
have been presented in section 4.2.2. We have already discussed that Z' decay patterns 
with less than 2 charged leptons have no scope at the LHC, being overwhelmed by the 
backgrounds. On the other side, non-standard 2 leptons decay modes are meagre with 
respect to the prompt di-lepton decay modes, and their analysis would be a more refined 
study than the one of section 3.2.2. We think that only the 2£ + 4j deserves a closer 
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look at, as backgrounds can be effectively suppressed and neutrino masses reconstructed, 
although the very high jet multiplicity might complicate the analysis. Ref. [74] showed 
that the task could be feasible, although in the slightly different context of a 4*^^ family 
extension of the SM, where the heavy neutrino pair is produced via SM Z and Higgs 
bosons, therefore in a completely different kinematic region. 

The 4 leptons decay mode of eq. (4.25) has the advantage of being rather background 
free. However, its corresponding BR is small [0(0.5%)]. In the context of the B — L 
model it has been studied in Ref. [60], showing that it can be observable for £ = 300 
fb~^, and the intermediate heavy neutrino mass poorly measurable. 

As this last reference was in progress, we looked in details at the signature with 3 
charged leptons, in eq. (4.24), as first reported in Ref. [47]. For the price of considering 
bigger backgrounds, the tri-lepton signature has higher cross section, thereby a lower 
integrated luminosity is sufficient for its discovery, if compared with the 4 leptons mode. 
It is crucial to note that the 2 jets come always from a W boson, and this is fundamental 
for reducing the backgrounds, as we will show. A pictorial representation of the tri-lepton 
decay mode of the Z' boson is in figure 4.7. 




Figure 4.7: Feynman diagram for the tri-lepton decay mode of the Z' boson (i = e, ^i). 

The tri-lepton decay mode also offers a nice framework to study the heavy neutrino 
properties. It is remarkable that in this signature the heavy neutrino flavours can be 
univocally identifled, as the charged lepton that comes along the jets can be unambigu- 
ously determined, and so its flavour ^. Just by counting, the relative BR and the possible 
mixing in the heavy neutrino sector can also be potentially measured. Nonetheless, we 
leave these features for future investigation and focus here in delineating the strategy for 
the discovery, in the approximation of degenerate and non-mixed heavy neutrinos. For 
simplicity, we also limit ourselves to the case without leptonically decaying r's, leaving 
also this case for future investigation. 

The chosen benchmark points, at the LHC for = 14 TeV and L = 100 fb~^, 
are for Mz' = 1.5 TeV, g[ = 0.2 [for which, from flgure 3.2b, a{pp Z') = 0.272 
pb], and two different heavy neutrino masses: m^^^ = 200 GeV, for which the heavy 
neutrino pair production cross section via Z' boson is 46.7 fb [also, from flgure 3.4, 
second panel, BR(Z' — )• Yl^h^'h) ~ 17.1%], and = 500 GeV, for which, again, 
a{pp ^Z' ^Y. ^hyh) = 23.4 fb [and BR(Z' ^ ^ v^Vh) ~ 8.7%]. 

^Notice that the same is true also for the 2t + 4j mode, but in general it is not for tlie other 
multi-lepton decay modes. 
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These two benchmark points provide two very kinematically different examples. In 
the first case, the heavy neutrinos are much lighter than the Z' boson, producing highly 
boosted events. In the second case, their mass is comparable to Mz'/2, hence close to 
their production threshold, resulting in minimal boost. From a merely kinematic point 
of view, all other cases will be somewhere between these two. 

When the heavy neutrino decays via the £^W^ mode, with a subsequent leptonic 
decay of the W^, the charged pair of leptons can carry an invariant mass equal to or 
lower than the heavy neutrino mass, with the maximum invariant mass configuration 
occurring when the light neutrino is produced at rest, so that the edge in this distribution 
corresponds to the I'h mass. A peak in such a distribution corresponding to the SM-like 
Z boson, coming from the I'Z decay mode for the heavy neutrino, will also be present 
in this distribution (see, e.g., plots in Ref. [47]). 

While the invariant mass distribution can provide some insights into the mass of 
the intermediate objects, this is not the best observable in the case of the tri-lepton 
signature, because the final state light neutrino escapes detection. A more suitable 
distribution to look at is the transverse mass defined in Ref. [75], i.e., 

mf. = (^^AP{vis) + P^{vis) + - (^Pxivis) +fTf > (4-26) 

where {vis) means the sum over the visible particles. For the final state considered 
here, if we sum over the 3 leptons and the 2 jets, this distribution will peak at the Z' 
mass. We can also see evidence for the presence of a heavy neutrino by just considering 
the 2 closest (in Ai?, where Ai? = Ar/^ + A(/>2) leptons and the missing transverse 
momentum, since this is the topology relevant to a Vh decay. The results show that this 
transverse mass peak for the heavy neutrino is likely to be the best way to measure its 
mass [47]. The striking signature of this model is that both of the above peaks occur 
simultaneously. 

We stress that by choosing event by event the 2 closest (in A.R) leptons, the peak 
corresponding to the heavy neutrino is well reconstructed, as clear from figure 4.12. 
As we pointed out previously, this prescription also allows for the almost unambiguous 
identification of the single lepton coming from the heavy neutrino that decays semi- 
hadronically. Its flavour can therefore be identified. 

In the evaluation of the background we considered three sources (including generation 
cuts, to improve efficiency): 

- WZjj associated production {a^i = 246.7 Po, £ = e,iJ,,T; A.Rjj > 0.5, Pj^^ > 
40 GeV, \7]j,^^\ < 3); 

- ti pair production {a2e = 29.6 pb, i = e, fi {b quark not decayed); QCD scale 
= Mt/2 to emulate the NLO cross section; no cuts applied); 

- tiiu associated production {a^i = 8.6 ib, i = e, ^, r; QCD scale = Pj' > 20 
GeV). 
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In the case of WZjj associated production, 3 leptons come from the subsequent 
leptonic decays of the two gauge bosons. This is the main source of background. From 
ti pair production one can obtain 2 isolated leptons via the decay of the VF^'s pro- 
duced from (anti)top decays and produce a further lepton from a semi-leptonic B meson 
(emerging from one of the 6's) decay. This lepton though will not be isolated, because 
of the large boost of the b quark, and this property is used to suppress this potential 
background. (In practice, the requirement of isolated leptons, made on the angular sepa- 
ration of a lepton from the nearest jet, strongly suppresses ti.) Finally, iiii' will produce 
3 isolated leptons, 2 of which from the VF^'s arising from the (anti)top decays (as above), 
resulting in a significant background despite the small production cross section. 

In the following subsections we present the result of the analysis at the parton level 
and at the detector level, respectively, of the tri-lepton signature and of the backgrounds. 
We define the cuts and show their effect in enhancing the signal-to-background ratio. 

4.2.3.1 Parton level analysis 

This subsection presents the analysis of signal and backgrounds for the tri-lepton signa- 
ture at the parton level. 

The first set of cuts we use is designed to impose generic detector angular acceptances, 
lepton and jet transverse momentum minimal thresholds and to provide isolation for 
leptons and jets: 

Selection 1 





< 


2.5, 




< 


3; 




> 


15 GeV, 


p'^ 


> 


10 GeV, 


Jl,2 


> 


40 GeV; 




> 


0.5 V£ = 1 . . . 3, j 




> 


0.2 V£,/ = 1...3, 


^Rkj 


> 


0.5. 



Special care should be devoted to the treatment of the tt background, given its large 
production rate. However, as previously mentioned, we expect that this noise can be 
eliminated efficiently by enforcing a suitable lepton-jet separation. The impact of the 
first set of cuts on the signals and ti background is illustrated in table 4.1. The ARij 
requirement is indeed extremely effective and reduces this background by a factor of 
2 • 10~^. The loss of signal due to this cut is instead minimal. Also note that the signal 
events with the smaller boost have a higher efficiency for passing the angular isolation 
cuts. 
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Figure 4.8: Signal (rriy^ = 200 GeV) and background distributions after the Selection 1 
cuts. (Events per .if = 100 fb^^.) 



Figures 4.8-4.9 show the distributions in Mjj, M^^gji M,^ after Selection 1 

cuts, for the signal with the two heavy neutrino masses, 200 and 500 GeV, that we are 
considering, and for the backgrounds. 

In signal events both jets come from the therefore we apply the following con- 
straint: 

Selection 2 

\Mjj - Mw\ < 20 GeV. (4.28) 

After the application of this cut the other distributions considered are shown in 
figures 4.10 and 4.11 for the 200 and 500 GeV heavy neutrino masses, respectively (here, 
we now also show the difference between the and M^jj distributions). From these 



57 




100 200 300 400 

M.. (GeV) 




200 400 600 800 

M,.j (GeV) 




500 1000 1500 2000 

M\,,j(GeV) 




500 1000 1500 

M\,(GeV) 



Figure 4.9: Signal (my^ = 500 GeV) and background distributions after the Selection 1 
cuts. (Events per = 100 fb'^.) 



plots it is clear that transverse mass Mg^gj provides good discrimination between signal 
and background. The following cut is then used to further suppress the background: 



Selection 3 



\M^i,2j -Mz'\< 250 GeV. (4.29) 

After this set of cuts we end up with a very clean signal for both a 200 and 500 
GeV I'h mass in the di-lepton transverse mass distribution, in fact practically free from 
background, as shown in figure 4.12. Notice that this variable was formed from the 
two closest (in ARu) leptons since they are likely to originate from the same boosted 

In order to establish the signal, we finally select events around the visible peak, 
by requiring: 
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Cut 


m^^ = 200 GeV 
# of events Efficiency % 


m^^ = 500 GeV 
# of events Efficiency % 


tt 

# of events Efficiency % 


No cuts 
77 cuts 
AR+Pt cuts 


482.0 100 
346.0 71.8 
68.0 14.1 


239.0 100 
171.0 71.4 
73.7 30.8 


1.28-10** 100 
5.1 • 10^ 43.9 
99.7 0.014 



Table 4.1: Efficiencies of the Selection 1 cuts for the two benchmark signals and the tt 
background, for events with three leptons and with two or more jets in the final state for 
= 100 fb~^ . In the case of ARjj < 0.5 partons were merged into one 'jet' at the very 
beginning of the selection. 



Selection 4 

0<M^g< 250 GeV or 400 GeV < M^g < 550 GeV, (4.30) 

depending on the benchmark signal under consideration. A fit of the signal in these 
regions will finally give the heavy neutrino mass. The efficiencies of the Selection 1-4 
cuts, are given in table 4.2. 



m., = 200 GeV 



Cuts 


Signal 


Efficiency 


WZjj 


Efficiency 


tt 


Efficiency 


tt£u 


Efficiency 






events 


(%) 


events 


(%) 


events 


(%) 


events 


(%) 




1 


68.0 


100 


5875 


100 


99.6 


100 


89.1 


100 


0.87 


2 


68.0 


100 


498.8 


8.5 


5.38 


5.4 


19.3 


21.8 


2.97 


3 


58.8 


86.5 


10.58 


12.7 





0.8 


0.0667 


2.2 


18.0 


4 


56.0 


94.1 


4.48 


67.6 





56.4 


0.0305 


64.8 


26.3 










m^^ = 


500 GeV 








Cuts 


Signal 


Efficiency 


WZjj 


Efficiency 


tt 


Efficiency 


ttiv 


Efficiency 


s/^/b~ 




events 


(%) 


events 


(%) 


events 


(%) 


events 


(%) 




1 


73.6 


100 


5875 


100 


99.7 


100 


89.1 


100 


0.95 


2 


73.6 


100 


498.8 


8.5 


5.38 


5.4 


19.3 


21.8 


3.22 


3 


68.8 


93.4 


10.58 


12.7 





0.8 


0.0667 


2.2 


21.1 


4 


46.3 


66.0 


2.879 


27.1 





8.7 


0.00952 


10.1 


27.6 



Table 4.2: Signal (m^j^ = 200 GeV at the top and m^^^ = 500 GeV at the bottom) and 
background events per = 100 fb^^ and efficiencies following the sequential application 
of Selection 1-4 cuts. 

We should stress that during all the step of our analysis we have included detector 
effects and applied these to the parton level events simulated using calchep. We have 
used Gaussian smearing of leptons and quarks energies to simulate a typical electromag- 
netic energy resolution given by 0.15/\/!E and a typical hadronic energy resolution given 
by 0.5/\/^, to mimic the ATLAS [64] and CMS [63] detector performances. The reso- 
lution in missing transverse momentum was derived in turn from the above mentioned 
smearing of leptons and quarks energies. The final figure 4.12 and table 4.2 present 
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Figure 4.10: Signal (ttIi,^ = 200 GeV) and background distributions after the Selection 1 
and 2 cuts. (Events per ^ = 100 fb~^.) 



results taking into account these effects. One should also mention that, since we deal 
with a multi-lepton signature in the final state, detector resolution effects are typically 
smaller or of the order of the chosen 20 GeV bin width for the final transverse mass 
distributions. Finally, we have verified that all results presented here are actually very 
stable against the implementation of electromagnetic and hadronic calorimeter energy 
resolution effects. 

In the next section we will present the study of the tri-lepton signature in a more 
realistic framework, with a full detector simulation in which also jet hadronisation and 
showering are included. The main difference with respect to the analysis presented 
here concerns the jet multiplicity, that at the parton level we cannot emulate. Also, 
the lepton isolation requirement will be treated in a more comprehensive way than the 
simple approximation used here. In the analysis, we will adopt the same strategy and 
we will use cuts emulating those presented here. 

It is interesting to note here that, in the strategy described in this section, we did 
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Figure 4.11: Signal (m^,^ = 500 GeV) and background distributions after the Selection 1 
and 2 cuts. (Events per ££ = 100 fb~^.) 



not impose any requirements on the sign of the leptons in order to reconstruct the heavy 
neutrinos, nor on the overall balance of their total sign. Also, flavour violating processes, 
clear indication of new physics, are possible. We decided to not include any of these 
features in our analysis to be the most model-independent possible, highlighting that 
simple cuts are enough for heavy neutrino discovery. 

4.2.3.2 Detector level analysis 

This subsection presents the analysis of signal and backgrounds for the tri-lepton sig- 
nature at the detector level. Its aim is to validate the analysis strategy outlined in the 
previous subsection within a more realistic framework. The parton level events gener- 
ated for the previous analysis have been here hadronised and showered with pythia, 
version 6.2.40 [76]; the CMS detector has been emulated with the fast detector simulator 
DELPHES, version 1.9 [77]; the jets have been reconstructed with the built-in SIScone 
algorithm, version 1.3.3 [78], with cone size 0.5. 
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Figure 4.12: Signal (m^^ = 200 GeV, top panel, and m^^ = 500 GeV, bottom panel) 
and background distributions after the Selection 1, 2 and 3 cuts. (Events per S£ = 100 



Most of the Selection 1 cuts of eq. (4.27) is built-in in the object definition. The 
minimum p^ for leptons is 10 GeV, while the detector angular acceptance is \r](\ < 
2.5(2.4) for electrons(muons). Regarding the jets, minimum px is 20 GeV and \rjj\ < 3. 
Our choice here of a cone size of 0.5 matches the last statement of eq. (4.27). The 
tri-lepton signature is finally defined by the presence of exactly 3 leptons in the final 
state, both electrons and muons, and we require at least 2 jets. The 2 most energetic 
jets are requested having pT > 40 GeV. A summary of the Selection 1 cuts here used is 
in eq. (4.31). 
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Selection 1 





- 


3, 


#J 


> 


2; 




< 


2.5(2.4) 


\m\ 


< 


3; 


■^-^1,2,3 


> 


10 GeV 


Jl.2 


> 


40 GeV 


jet cone radius 




0.5, 


isolated leptons 




true. 



The DELPHES package allows for the lepton isolation to be more carefully defined. 
To match the CMS definitions, we define a lepton to be 'isolated' if no other tracks 
in a cone of radius 0.3 within the tracker has > 2 GeV. The latter appears more 
restrictive than the corresponding ARw requirement we used at the parton level. 

In the previous section we highlighted that the requirement of isolated leptons was 
crucial to suppress the tt background, in which a further lepton stems from a b quark 
decay. Table 4.3 shows the effect of the Selection 1 cuts for both signal points and the tt 
background. We see that, in the detector framework, the requirement of isolated leptons 
reduces the background by a further factor 2.4%, an order of magnitude less efficient 
than at the parton level. Regarding the signal, it still holds that the isolation requirement 
suppresses more the low rrij/^ point, where final state objects are more boosted. 



Cut 


m^^ = 200 GeV 
# of events Efficiency % 


m^^ = 500 GeV 
# of events Efficiency % 


tt 

# of events Efficiency % 


Basic cuts 
Isolation cut 


429.4 100 
178.3 41.5 


214.4 100 
117.9 55.0 


4.31 • 10'' 100 
1.04 • 10* 2.4 



Table 4.3: Efficiencies of the isolation cut for the two benchmark signals and the tt 
background, for events with three leptons and with two or more jets in the final state, 
for 5£ = 100 fb~^. 'Basic cuts' refer to the Selection 1 cuts of eq. (4-31) without the 
isolation requirement. 



At the parton level there is no ambiguity concerning the jets, since, for the signal, 
two and only two are present. At the detector level, instead, the average jet multiplicity 
is around 4, thus we need to identify the two jets to include in our distributions. As 
pointed out, for the signal only, the topology of the tri-lepton signature imposes the jets 
to come exclusively from a W boson. Hence, on a event by event basis, we can choose 
the pair of jets that better reconstructs the W boson (meaning, whose invariant mass 
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is closest to Mw = 80.4 GeV). Selection 2 cut is hence defined as a cut in the invariant 
mass distribution of this jet pair, jets that will subsequently enter in the definition of 
our distributions, such as Mj^2j- To avoid to spoil considerably the signal, while still 
rejecting most of the background, the cut of eq. (4.28) has to be relaxed to 

Selection 2 



\Mjj - Mh^I < 30 GeV. (4.32) 

It shall be noted that this is not the most efficient way to identify the pair of jets to 
be used. We comment on this at the end of this section. 

Once the pair of jets has been identified, we can use it to plot the Mj^2j distribution, 
shown in figure 4.13 (summing up all backgrounds and signal, and for the signal only, 
in the top panels and in the bottom panels, respectively). The events for the signal are 
visible as an excess over the background, peaking at the Z' boson mass, as clear from 
figures 4.13c and 4.13d (in figures 4.13a and 4.13b the peak is not evident only due to 
the log scale). 
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Figure 4.13: Distribution of Mj^2j f(^f signal (4-13a) m^^ = 200 GeV and (4-13b) rriy^ = 
500 GeV, summed with the backgrounds, and for signal only, in (4.13c) and (4.13d), 
respectively, after the Selection 2 cuts. (Events per ^ = 100 fb~^). 
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As in the previous section, we define the following cut to further suppress the back- 
ground: 



Selection 3 



|M|^2i - Mz'\< 250 GeV, (4.33) 

in which we consider the Z' mass as known, measured from Drell-Yan processes (see 
section 3.2.2). The distribution after this cut is shown in figure 4.14. 
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Figure 4.14: Distribution o/M|^ for signal (4-14C') ''t^uh ~ ^0*-* GeV and (4-Hh) 
500 GeV, superimposed to the backgrounds after the Selection 3 cuts. (Events per J£ 
100 fb~^}. 



In order to establish the signal, we finally select events around the visible peak, 
by requiring: 

Selection 4 



100 < mI(^ < 250 GeV or 300 GeV < M^^ < 550 GeV, (4.34) 

depending on the signal benchmark under consideration. Again, a fit of the signal in 
these regions will finally give the heavy neutrino mass. The efficiencies of the Selection 
1-4 cuts are given in table 4.4. Notice that the number of events for the signal, especially 
at the beginning, is bigger than at the parton level. This is because all the other multi- 
lepton signatures can impinge in the tri-lepton one, but they do not survive the whole 
chain of cuts. This summary clearly confirms the feasibility of the extraction of the 
signal even with less than 100 fb~^ of integrated luminosity. 

As we specified earlier, table 4.4 shows that the procedure we used to identify the 
pair of jets to include in our distributions is not the most efficient. In fact, the Selection 
2 cut has a lower efficiency for the signal (~ 40%, irrespectively of the benchmark point) 
than for the background (~ 70% 90%). The reason for this is that we are shaping 
the background as the signal, forcing to select those jets that mimic it. Nonetheless, 
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m., = 200 GeV 



Cuts 


Signal Efficiency 
events (%) 


WZjj Efficiency 
events (%) 


tt Efficiency 
events (%) 


ttiu Efficiency 
events (%) 


s/Vb 


1 
2 
3 
4 


178.3 100 

74.2 41.6 

34.3 46.3 
31.9 92.8 


3154 100 
2189 69.4 
27.96 1.28 
11.51 41.18 


1.04- 10-* 100 
9.36 • 10^ 89.8 
<1 <0.1 



112.5 100 
84.4 75.0 
0.286 0.34 
0.086 30.0 


1.52 
0.69 
6.46 
9.36 


ruy^ = 500 GeV 


Cuts 


Signal Efficiency 
events (%) 


WZjj Efficiency 
events (%) 


tt Efficiency 
events (%) 


ttiu Efficiency 
events (%) 


S/VB 


1 
2 
3 
4 


117.9 100 
48.9 41.5 
33.5 68.5 
22.1 66.0 


3154 100 
2189 69.4 
27.96 1.28 
9.05 32.4 


1.04-10* 100 
9.36 • 10=* 89.8 
<1 <0.1 



112.5 100 
84.4 75.0 
0.286 0.34 
0.114 40.0 


1.01 
0.45 
6.31 
6.50 



Table 4.4: Signal (m^^ = 200 GeV at the top and m^,^ = 500 GeV at the bottom) and 
background events per = 100 fb^^ and efficiencies following the sequential application 
of Selection 1-4 cuts. 

mimicking the W boson mass is not sufficient to imitate all the signal properties. In 
fact, the subsequent Selection 3 cut strongly suppresses only the background. 

4.3 Conclusion 

In this chapter, we have presented the results of our investigation of the fermion sector. 

First, we have summarised the existing experimental constraints on the neutrino 
masses. We have also discussed the approximation used in the phenomenological study, 
i.e., to take heavy neutrinos relatively light (with respect to the Z' boson) and degener- 
ate. 

Then, we have presented a detailed study of the heavy neutrino properties (cross 
sections, decay widths and BRs). We have shown that the production cross sections 
via the Z' boson, the channel of interest, are up to 50 fb for the LHC in its early stage 
{^/s = 7 TeV), value that induced us to consider this channel as viable only for the LHC 
in its designed performances {^/s = 14 TeV and full luminosity). 

Regarding the decay properties of the heavy neutrinos, we have shown that the 
decay into a SM gauge boson (together with a lepton) is preferred. This partial width 
is very small, and, hence, the heavy neutrino has a very small intrinsic width. The 
decay patterns of the heavy neutrino pair production via Z' boson have been analysed 
in details. They can be thought of as multi-lepton (and multi-jet) decays of the Z' 
boson. Up to 4 charged leptons in the final state can be produced. 

We have concluded that for a large portion of the parameter space, the heavy neu- 
trinos are rather long-lived particles, so that they produce displaced vertices in the LHC 
detectors, that can be distinguished from those induced by h quarks. In addition, from 
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the simultaneous measurement of both the heavy neutrino mass and decay length one 
can estimate the absolute mass of the parent light neutrino, for which at present, only 
limits exist. 

To address the neutrino mass measurement, we have chosen to study the so-called 
tri-lepton decay mode, in which the Z' boson decays into 3 charged leptons (either 
electrons or muons) and 2 jets, always stemming from a W boson. 

For benchmark scenarios of the B — L model, we have chosen two that should be 
accessible at the LHC, having a Z' boson mass and fermion couplings not far beyond 
the ultimate reach of Tevatron and LEP and displaying two extreme relative conditions 
between the Z' boson and heavy neutrinos, that is, one with the latter produced at rest 
and the other highly boosted in the direction of the Z' boson. 

A detailed parton level simulation of the aforementioned signal benchmarks has been 
performed, including all the relevant backgrounds. We have shown that rather generic 
cuts are already suitable for an efficient background rejection. This analysis strategy 
has, then, been validated with a full detector simulation of both signal and background. 
To do so, a strategy for the jet identification has been proposed, based on the signal 
topology at the parton level (in which, as mentioned, the jets come exclusively from the 
W boson). 

The transverse mass has been found to be the suitable distribution for the tri-lepton 
signature analysis. By carefully choosing what particles to consider for its definition, 
peaks for both the Z' boson and for the heavy neutrino emerge. Their simultaneous 
observation is a striking signature of the B — L model. In particular, the sharpness of 
the peak corresponding to the heavy neutrino, as well as its definition well above the 
background, as shown, makes the heavy neutrino mass well measurable. 

As a result, the tri-lepton signature holds a very good signal over background ratio 
that makes it observable at the LHC, for yfs = 14 TeV and £ < 100 fb~^. 



Chapter 5 

Scalar sector 



In this chapter we study the scalar sector of the B — L model. The new states here are 
two CP-even Higgs bosons, whose masses are free parameters. Another free parameter 
is the angle a that controls their mixing. 

Section 5.1 presents the constraints on the free parameters in this sector. We will 
first review in section 5.1.1 the existent experimental limits, coming from both direct and 
indirect searches at LEP. In section 5.1.2 is then presented the analysis of the theoretical 
bounds on the scalar sector, i.e., the so-called unitarity bound and the triviality and 
vacuum stability bounds, the latter two coming from the study of the RGEs. 

In section 5.2, the properties of the scalar sector are delineated (i.e., production cross 
sections, intrinsic widths, BRs), and the capabilities for Higgs discovery at the LHC are 
summarised. The main focus is on the impact of the gauge and fermion sectors. In 
fact, although no new consistent production mechanisms arise, the decay patterns in 
this model are rather peculiar, such as the decay into pairs of Z' bosons and into heavy 
neutrinos (besides the decay of the heavy Higgs boson into pairs of the light scalar 
bosons). Cross sections for some full processes are presented in section 5.2.3. 

Finally, we draw the conclusions for this chapter in section 5.3. 

The study of the unitarity bound is published in [3]. The RGE study (triviality and 
vacuum stability) is published in [40] . The analysis of the phenomenology of the Higgs 
bosons is published in [79]. 

5.1 Constraints 

In this section is delineated the viable parameter space for the new independent param- 
eters in the scalar sector, i.e., m/j^, m/jj, and the mixing angle a. 

First, we summarise the experimental constraints from LEP, from direct searches 
and from the analysis of the precision tests. Regarding the former, as discussed in the 
introduction, we did not pursue a complete beyond the tree-level analysis and we will 
refer to the literature. 

Finally, we present our analysis of the theoretical constraints, i.e., the unitarity 
bound and the constraints coming from the study of the RGEs. 
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5.1.1 Experimental constraints 

Past and current experiments have set limits on the parameters of the scalar sector 
in the SM as well as in various extensions of it: see, for example Ref. [1] for LEP and 
Ref. [80] for Tevatron. For the model discussed here, the relevant analysis is summarised 
in figure 5.1 (figure 10a in the LEP combined analysis of Ref. [1]), in which a generic 
overall factor ^ has been introduced. Such parameter is defined as the coupling(s) to 
the Z boson of the Higgs particle(s) in the considered extension normalised to the SM: 

(5.1) 

Uhzz 

hence it parametrises the deviations of the new model with respect to the SM. 




120 

mjj(GeV/c^) 

Figure 5.1: The 95% C.L. upper bound on ^ = gHZZ /gfjzz l-^l- -^^ B — L model, 
= cosa(sina) for H = /ii(/i2). The solid line represents observed values; the dotted 
line represents values expected for the background. 



In the minimal U{1) extension of the SM, the argument of this Thesis, two real scalar 
degrees of freedom exist: the one coming from the Higgs singlet, required to break the 
extra U{1)b-l gauge factor (and therefore giving the Z' gauge boson a mass), and the 
one coming from the Higgs doublet, required to break the EW gauge symmetry to give 
masses to the W and Z bosons. In all generality, these two scalars will mix (as, for 
example, in any other extension of the SM with one scalar singlet in addition to the 
Higgs doublet). With reference to eq. (2.28), hi is the lightest eigenstate, that couples 
to the Z boson proportionally to cos a, /i2 is the heaviest one, that couples to the Z 
boson proportionally to sin a. Hence, the LEP lower bounds on the scalar masses of 
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the U{1)b-l extension here considered are read straightforwardly from figure 5.1 by 
considering: 



i.e., the hmit for /ii(/i2) is extracted by considering ^ as the cosine(sine) of the mixing 
angle in the scalar sector [see eq. (2.28)]. 

Figure 5.1 shows the lower bound on the Higgs boson mass as a function of The 
SM Higgs boson is recovered by the condition = 1. We see that we can have significant 
deviations from the SM Higgs boson mass limit, nih > 114.4 GeV, only for values of the 
angle a > 11/4:, for the lightest state hi. For example, for a = it/3, the LEP limit on 
the lightest Higgs state reads as rrih-^ > 100 GeV. That is, in this model, a light Higgs 
boson with mass smaller than the SM limit can exist only if it is highly mixed, i.e., the 
light Higgs boson is mostly the singlet state. For the same value of the angle, the limit 
for m/j2 is more stringent than the condition nifi^ > rrij-i^, in fact for a = vr/S, m^^ > 114 
GeV must be fulfilled. 

Direct searches have put limits on scalar masses in the B — L model below the 
LEP limit of 114.4 GeV (for the SM). The existence of scalar bosons heavier than this 
limit can also be restricted. One way is by analysing their impact on the EW precision 
measurements [49], encoded in the 5, T and U parameters [81]. In the case of a SM 
scalar sector augmented by just one scalar singlet, the constraints on the scalar masses 
are analysed in Ref. [2], and summarised in figure 5.2. 

Generally speaking, the scalar sector of the B — L model is equivalent to the one of 
generic scalar singlet extensions of the SM, so that we can straightforwardly interpret 
the results of Ref. [2] as valid also here. It ought to be noticed, though, that this 
statement relies on the simplification that the impact of the extra matter content of 
the B — L model, i.e., Z' boson and heavy neutrinos, on the precision parameters is 
negligible. This is indeed the case for the direct contribution of the Z' boson, as it has 
vanishing mixing with the SM Z boson at the tree level [51] (and the running on g, 
turning this mixing on at the NLO, is a two-loops effect). Nonetheless, the Z' boson 
and the heavy neutrinos can still alter the scalar boson total widths entering in the 
precision parameters' calculation. An increment in the total width will relax the results 
of Ref. [2], although marginally as the total width changes are mild for values of the 
mixing angle bounded therein. 

With reference to [2], the mixing matrix V corresponds to the matrix in eq. (2.28), 
so that I Vol I = sin a. Important consequences of this analysis are that there is an upper 
bound on the light scalar, ?n/jj < mf^_^'^^ = 165 GeV, and that the complete inversion 
in the scalar sector, i.e., for a = vr/2, is forbidden by precision data for m/^j > m^^^^ . 

5.1.2 Theoretical constraints 




for H = hi 
for H = h2 



(5.2) 



In this section we present the analysis of the theoretical constraints on the new param- 
eters in the scalar sector, namely the unitarity bound and those coming from the study 
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1 Singlet 




m„ (GeV) 

Figure 5.2: Allowed region (at 95% confidence level) in a model with one additional 
singlet in addition to the usual SU{2)l doublet. The lightest (heavier) scalar mass is 
iTT-o (n^i) o,nd the mixing matrix is defined in eq. (A) of Ref. [2]. In our notation, 
\Vqi\ = sina. The region below the curves is allowed by fits to S, T and U. 

of the RGEs. The unitarity assumption [3] rehes on demanding that the cross sections 
for the scattering of all the particles in the model are unitary, i.e., with an occurring 
probability less than one. The RGE constraints [40] come from the model to be well- 
defined at one loop, assuming that the evolved parameters do not hit any Landau pole 
or destabilise the vacuum of the model. The latter request is known as the vacuum 
stability condition, while the former is the triviality condition. 

5.1.2.1 Unitarity bound 

It is generally not possible within the SM framework (or any of its non-supersymmetric 
extensions encompassing the Higgs mechanism) to predict the mass of the Higgs boson. 
Hence, several theoretical methods have been developed to constrain its value (see [82], 
[83], [84]). For example, to stay with the SM, the pioneeristic work of Ref. [83] showed 
that, when the Higgs boson mass is greater than a critical value of around 1 TeV (known 
as unitarity bound), the spherical partial wave describing the elastic scattering of the 
longitudinally polarised vector bosons at very high energy {^/s — )■ oo) violates unitarity 
at the tree level and the theory ceases to be valid from a perturbative point of view. 

In the high energy limit (s » ^vF±zz')' amplitude involving the (physical) 
longitudinal polarisation of gauge bosons approaches the one involving the (unphysical) 
scalar one (equivalence theorem, see [85]), and, following the BRS invariance [86], the 
amplitude for emission or absorption of a 'scalarly' polarised gauge boson becomes equal 
to the amplitude for emission or absorption of the related Goldstone boson. Since it is 
the elastic scattering of longitudinally polarised vector bosons that gives rise to unitarity 
violation, the analysis of the perturbative unitarity of two-to-two particle scatterings in 
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the gauge sector can be performed, in the high energy hmit, by exploiting the Goldstone 
sector. 

Hence, in the high energy hmit, we can substitute the vector boson and Higgs bo- 
son sectors with the related (would-be) Goldstone and Higgs boson sectors. We will 
therefore focus on the scalar interacting Lagrangian of the Higgs and would-be Gold- 
stone sectors (in the Feynman gauge), i.e., the scalar Lagrangian of eqs. (2.2), (2.3) and 
(2.13), neglecting here the gauge couplings in the covariant derivative of eq. (2.7), and 
we will calculate tree-level amplitudes for all two-to-two processes involving the full set 
of possible (pseudo) scalar fields ^. The relevant Feynman rules are in section A. 3, where 
Vp is the Goldstone counterpart of the vector V. 

5.1.2.1.1 Evaluation of the unitarity bound 

As already intimated, the equivalence theorem allows one to compute the amplitude of 
any process with longitudinal vector bosons Vl {V = W^, Z, Z'), in the limit my <^ s, 
by substituting each one of them with the related Goldstone bosons v = w^,z,z', and 
its general validity is proven [85]; schematically, if we consider a process with four 
longitudinal vector bosons: M.{VlVl — >• VlVl) = 'M{vv — )• vv) + 0{m'y/s). 

Given a tree-level scattering amplitude between two spin-0 particles, M(s,^), where 
6 is the scattering (polar) angle, the partial wave amplitude with angular momentum J 
is given by 



where Pj are Legendre polynomials. It has been proven in Ref. [87] that, in order to 
preserve unitarity, each partial wave must be bounded by the condition 



It turns out that only J = (corresponding to the spherical partial wave contribution) 
leads to some bound, so we will not discuss the higher partial waves any further. 

We have verified that, in the high energy limit, only the four-point vertices (related 
to the four-point functions of the interacting potential) contribute to the J = partial 
wave amplitudes, and this is consistent with many other works that exploit the same 
methodology (for example, see [88; 89; 90]). The main contributions come from the 
zz — )• zz and z' z' — )• z' z' channels. To a lesser extent, also hihi — )• hihi and h2h2 — >■ /i2^2 
play a relevant role, being equal to the most constraining one among the main channels 
in some regions of the parameter space (for a — and a — 7r/2, respectively). Due to 
this, we will not discuss them any further (details of this can be found in Ref. [3]). 

^Moreover, while evaluating scalar bosons' scattering amplitudes, it has been explicitly verified that, 
in the search for the Higgs boson mass limits, the contribution that arises from intermediate vector 
boson exchange is not relevant. 




(5.3) 




(5.4) 
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Moving to the results, the unitarity bound is, in general, a function of all the param- 
eters that take part to the four-point vertices, i.e., the scalar masses, the mixing angle 
a and the singlet VEV x. It is convenient to define a 'high-mixing' and a 'low-mixing' 
domains in the mixing angle. The former is defined by noticing that for a value of the 
angle in the range 



the allowed parameter space [in the {nih^-mh^) plane] is completely defined by the 
zz —7- zz channel. The 'low-mixing' domain is the complementary region. For instance, 
since x > 3.5 TeV as shown in section 3.1.1, due to the LEP bound [51], if we choose 
exactly x = 3.5 TeV, the high-mixing domain, in this case, is the one for 0.07 < a < tx 12 
(and, conversely, the low- mixing one is the interval < q < 0.07). 

Figure 5.3 shows the allowed regions in the (m/i^-m/ij) plane, for a value of the 
angle in the low-mixing domain (figure 5.3a), and for values in the high-mixing one 
(figures 5.3b, 5.3c and 5.3d), for several values of the singlet VEV x (affecting only the 
limit coming from the z' z' — )• z' z' channel). 

We see that in both cases, as expected, the upper bound on the light Higgs boson 
mass does not exceed the SM one (which is around 700 GeV [87]), and it approaches 
the experimental lower limit from LEP (according to [1]) as the heavy Higgs boson mass 
increases. This is because the two Higgs bosons 'cooperate' in the unitarisation of the 
scattering channels, so that, if one mass tends to grow, the other one must become 
lighter and lighter to keep the scattering matrix elements unitarised. 

In the high- mixing domain, as shown in figures 5.3b, 5.3c and 5.3d (for a = 0.1, 
a = 7r/4, a = 0.97r/2, respectively ^), the allowed region is completely determined by 
the zz — >• zz scattering, irrespectively of the value of the singlet VEV x. The maximum 
allowed value for the heavy Higgs boson mass only depends on the mixing angle 



Only in the low-mixing domain, as in figure 5.3a (for a = 0.01), we are able to 
appreciate some interplay between the two scattering processes in setting the unitarity 



a heavy Higgs boson with mass of more than 10 TeV, the z' z' — t- z' z' channel sharply 
limits the allowed mass region, with a 'cut-off' on the maximum m/jj almost insensible 
to the light Higgs boson mass (and to the value of the mixing angle), that scales linearly 
with the singlet VEV x 




(5.5) 




(5.6) 



bound. In this case, while the zz 



zz scattering channel would allow the existence of 




(5.7) 



^For the last of these values of the mixing angle, the lower limit from LEP experiments on the light 
Higgs boson mass is > 40 GeV [1], while for the other ones it is almost equal to the SM lower limit 
(jnh-i^ > 114.4 GeV), as illustrated in figure 5.3. 
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(c) (d) 
Figure 5.3: Higgs bosons mass limits in the B — L model coming from the unitarity 
condition \Re{ao)\ < ^ applied to the zz — )• zz and z'z' — )• z'z' scatterings for several 
values of x: (5.3a) for a = 0.01, (5.3b) for a = 0.1, (5.3c) for a = tt/A, and (5.3d) 
for a = 0.9 tt/2. The (blue) horizontal shadowed region is forbidden by our convention 
m/i2 > TTi/jj. The (red) vertical shadowed region is excluded by the LEP experiments. 

This is pictorially shown in figure 5.3a: the sohd green area represents the allowed 
portion of the m/^^-m/ij space at x = 3.5 GeV, that, at x = 7 TeV, increases until the 
(green) crossed shadowed region, to relax even further to the (green) single line shadowed 
region for x = 35 TeV. 

To summarise, given a value of the scalar mixing angle a, the upper bound on the 
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light Higgs boson mass varies between the SM hmit and the experimental lower limit from 
LEP as the upper bound for the heavy Higgs boson mass increases. Moreover, when 
a assumes values in the high-mixing domain, the strongest bound comes exclusively 
from the z-boson scattering, independently from the chosen singlet VEV x, while, in the 
low-mixing domain, the z' boson scattering can also be important, imposing a cut-off, 
linear in x, on the heavy Higgs boson mass (when more constraining than the zz — )■ zz 
scattering). 

In the following subsection we investigate the constraints coming from the analysis 
of the RGEs, collected in appendix B. 

5.1.2.2 Renormalisation group equations 

The RGE evolution can constrain the parameter space of the scalar sector in two com- 
plementary ways. From one side, the couplings must stay perturbative. This condition 
reads: 

0<Ai,2,|3|(Q')<l VQ'<g, (5.8) 

and it is usually referred to as the 'triviality' condition. Notice that A|3| = IA3I. On 
the other side, the vacuum of the theory must be well-defined at any scale, that is, to 
guarantee the validity of eqs. (2.11) and (2.12) at any scale Q' < Q: 

< Ai,2(gO and 4Ai(Q')A2(g') - Ai(Q') > V Q' < Q . (5.9) 

Eq. (5.9) is usually referred to as the 'vacuum stability' condition. In contrast to the 
SM, in which it is sufficient that the Higgs self-coupling A be positive, in the case of this 
model the vacuum stability condition [and especially the second part of eq. (5.9)] can 
be violated even for positive Ai^2,3- 

One should notice that our conventional choice ruh^ < m/jj, as noted previously, 
allows us to consider a and —a as two independent solutions, although the theory is 
manifestly invariant under the symmetry a — >• —a. These two solutions are complemen- 
tary, meaning that the region excluded by the choice rrih^ < nih^ at a certain value of the 
angle a is precisely the allowed one for the complementary angle 7r/2 — a. The special 
case Q = 7r/4 is symmetric, and corresponds to maximal mixing between the scalars. 
a = corresponds to a SM scalar sector totally decoupled from the extended one, and 
hi is the usual SM Higgs boson, q = 7r/2 is the specular case, in which /12 plays the role 
of the SM Higgs boson. 

Notice also that, again in contrast to the SM in which the gauge couplings have a 
marginal effect, in our case the RH neutrinos play for the extra scalar singlet the role 
of the top quark for the SM Higgs in the vacuum stability condition ^. Their RGEs 
are then controlled by the Yukawa couplings with a negative contribution coming from 

^Also notice that we have three RH neutrinos, as we have three colours for the top quark. However, 
they are Majorana particles rather than Dirac ones, so they carry half the (independent) degrees of 
freedom of the top quark. 
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g'^ [see eq. (B.9)]. Therefore, in some regions of the parameter space, the impact of 
the gauge sector is not marginal and can effectively stabilise the otherwise divergent 
evolution of the Majorana Yukawa couplings for the RH neutrinos. 

A final remark is in order about eq. (B.34), the evolution of A3, the mixing parameter 
of the scalar potential [see eq. (2.3)]. This RGE is almost proportional to A3 itself, so a 
vanishing boundary condition is almost stable ^. Non-proportional terms arise from the 
new gauge couplings (g and g[), i.e., deviations from the vanishing boundary conditions 
are of the order of the gauge couplings, hence quite small. They are particularly negligi- 
ble in the pure B — L model, as also g has a vanishing boundary condition, with a weak 
departure from it due to the mixing in the gauge coupling sector [29]. Nonetheless, other 
benchmark models in our general parameterisation could show different behaviours. 

The results we are going to present are obtained by analysing the RGEs in ap- 
pendix B. For their numerical study, we put boundary conditions at the EW scale 
on the physical observables: m/j^, m^j, a, v, Mz' , g'i,g,ml,'^'^ , that we trade for m, fj,, 
Ai, A2, A3, x,yj^23 using, for the relevant parameters therein, eq. (2.33). Notice that 
a\g denotes the mixing angle\coupling in the scalar\gauge sector. Where stated in the 
text, we impose boundary conditions on some parameters of the Lagrangian rather than 
on the physical observables. This is done for consistency of those studies. 

For the pure B — L model, object of the numerical analysis in this Thesis, the 

definition g = holds, and as a consequence, we also have that the B — L breaking VEV 

Mz' 

X can be easily related to the new Z' boson mass hy x = ^ ^ . Here we fix 5^ = 0.1. 

Regarding the heavy neutrinos, for simplicity we consider them degenerate and we fix 
their masses to mi'i^'^ = m^^ = 200 GeV (whenever not specified otherwise), a value 
that can lead to some interesting phenomenology (see section 4.2). The free parameters 
in this study are then nih^, m/ij, a and x. The general philosophy is to fix in turn some 
of the free parameters and scan over the other ones, individuating the allowed regions 
fulfilling the conditions of eqs. (5.8) and (5.9). 

5.1.2.2.1 Triviality and vacuum stability bounds 

Given the simplicity of the scalar sector in the SM, the triviality and vacuum stability 
conditions can be studied independently and they both constrain the Higgs boson mass, 
providing an upper bound and a lower bound, respectively. In more complicated models 
as the one considered here, it might be more convenient to study the overall effect of 
eqs. (5.8)-(5.9), since there are regions of the parameter space in which the constraints 
are evaded simultaneously. This is the strategy we decided to follow. 

Figure 5.4 shows the allowed region in the (m/j^-m/ij) parameter space, for increasing 
values of the mixing angle a, for fixed VEV x = 7.5 TeV and heavy neutrino masses 
TTij^^ = 200 GeV, corresponding to Yukawa couplings whose effect on the RGE running 
can be considered negligible. For a = 0, the allowed values for m/j^ are the SM ones 
and the extended scalar sector is completely decoupled. The allowed space is therefore 

^From the last line of eq. (2.33), setting A3 = corresponds to a = 0, but not vice versa. 
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the simple direct product of the two, as we can see in figure 5.4a. When there is no 
mixing, the bounds we get for the new heavy scalar are quite loose, allowing a several 
TeV range for m/jj, depending on the scale of validity of the theory. We observe no 
significant lower bounds (i.e., rrih^ > 0.5 GeV), as the impact of the RH Majorana 
neutrino Yukawa couplings is negligible. 




50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 



m,,(GeV) m^,(GeV) 

(c) (d) 

Figure 5.4: Allowed values in the rrih^ vs. space in the B — L model by eqs. (5.8) and 
(5.9), for (5.4a) a = 0, (5.4b) a = 0.1, (5.4c) a = it/A and (5.4d) a = 7r/3. Colours 
refer to different values of Q / GeV: blue (1^^), red (^lO'^j, green ("10^°;, purple (10^^) and 
cyan (^10^^ J. The shaded black region is forbidden by our convention m/j^ > f^hu while 
the shaded red region refers to the values of the scalar masses forbidden by LEP. Here: 
X = 7.5 TeV, = 200 GeV. 

As we increase the value for the angle, the allowed space deforms towards smaller 
values of m^.^. If for very small scales Q of validity of the theory such masses have 
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already been excluded by LEP, for big enough values of Q, at a small angle as a = 0.1, 
the presence of a heavier boson allows the model to survive up to higher scales for smaller 
hi masses if compared to the SM (in which just h\ would exist). Correspondingly, the 
constraints on rrih^ become tighter. Moving to bigger values of the angle, the mixing 
between hi and /12 grows up to its maximum, at a = 7r/4, where hi and /i2 both contain 
an equal amount of doublet and singlet scalars. The situation is therefore perfectly 
symmetric, as one can see from figure 5.4c. Finally, in figure 5.4d, we see that the 
bounds on m/jj are getting tighter, approaching the SM ones, and those for rrih^ are 
relaxing. That is, for values of the angle 7r/4 < a < vr/2, the situation is qualitatively 
not changed, but now /12 is the SM-like Higgs boson. Visually, one can get the allowed 
regions at a given angle 7r/2 — a by simply taking the transposed about the m/^^ = m^^ 
line of the plot for the given angle a. 

Per each value of the angle, we can then fix the lighter Higgs mass nih-^ to some 
benchmark values (allowed by LEP for the SM Higgs) and plot the allowed mass for the 
heavier Higgs as a function of the scale Q. This is done in figure 5.5, where the allowed 
masses are those contained between the same colour lines. Notice that here the VEV 
X is fixed to a different value, x = 3.5 TeV. The effects of changing the VEV x will be 
described in section 5.1.2.2.3. 

As previously noticed, the allowed range in 111^2 gets smaller as we increase the angle. 
Apart from the case a = where there is no dependency at all from m/j^, there is a 
strong effect of m/^^ on the bounds on nifi^- Not all the allowed regions (for m/jj) at a 
fixed hi mass are inside the regions that are allowed for a smaller ruh^ . This is true only 
for m/ij > 160 GeV. For smaller m/^^'s, the distortion in the allowed region constraints 
tightly m/j2 for the survival of the model to big scales Q. This is because such distortion 
is just towards smaller hi masses, see figure 5.4. 

Complementary to the previous study, we can now fix the light Higgs mass at specific, 
experimentally interesting ^, values, i.e., m/j^ = 100, 120, 160 and 180 GeV, and show 
the allowed region in the m^^ vs. a plane. This is done in figure 5.6. 

From this figures it is clear the transition of ^2 from the new extra scalar to the 
SM-like Higgs boson as we scan on the angle. As we increase nih^ (up to nifi^ = 160 
GeV), a bigger region in is allowed for the model to be valid up to the Plank scale 
(i.e., Q = 10^^ GeV, the most inner regions, in cyan). Nonetheless, such a region exists 
also for a value of the light Higgs mass excluded by LEP for the SM, m/,^ = 100 GeV, 
but only for big values of the mixing angle. No new regions (with respect to the SM) 
in which the model can survive up to the Plank scale open for rrih-^ > 160 GeV, as the 
allowed space deforms towards smaller values of ruh^ . 



^The chosen values maximise the probability for the decays hi — >■ bb, hi — >■ 77, hi — >■ W^W and 
hi — >■ ZZ, respectively (see figure 5.15a). 
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Figure 5.5: Allowed values (that are those between the same colour lines) for irih^ as a 
function of the scale Q in the B — L model by eqs. (5.8) and (5.9), for several values 
of rrih^ and (5.5a) a = 0, (5.5b) a = 0.1, (5.5c) a = vr/S and (5.5d) a = 7r/4. Also, 



3.5 TeV and m, 



200 Ge V. Only the allowed values by our convention m^^ > m^^ 



are shown. 



5.1.2.2.2 Heavy neutrino mass influence 

As already intimated, the RH neutrinos play for the extra scalar singlet the role of the 
top quark for the SM Higgs. This is particularly true for the vacuum stability condition, 
as the fermions in general provide the negative term that can drive the scalar couplings 
towards negative values. Figure 5.7 shows how the allowed regions in the (m/i^-m/jj) 
plane change for a RH Majorana neutrino Yukawa coupling y*^ = 0.2 (that for x = 3.5 
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Figure 5.6: Allowed values in the m/ij vs. a space in the B — L model by eqs. (5.8) 
and (5.9), for (5.6a) ruh, = 100 GeV, (5.6b) ruh, = 120 GeV, (5.6c) nif,^ = 160 GeV 
and (5.6d) ruh^ = 180 GeV. Colours refer to different values of Qj GeV: blue (^^), 
red (\^^ ), green ("10^°;, purple flO^^; and cyan (\^'^). The plots already encode our 
convention m/^j > and the shaded red region refers to the values of a forbidden by 
LEP. Here: x = 3.5 TeV, m^^ = 200 GeV. 

TeV correspond to m,y^ = 1 TeV), not negligible anymore. For y^^ = 0.4, the changes 
are even more drastic, shrinking the allowed region even further. 

The effect of having non negligible y^^ couplings is evident if we compare figure 5.7 
to figure 5.4. Notice that also the VEV x is changed (from 7.5 TeV to 3.5 TeV), but this 
is only responsible for the smaller upper bounds of in figures 5.7a and 5.7b. For 
small values of a it is evident our analogy between the top quark and the RH neutrinos, 
as now m/i2 has a sensible lower bound too. The analogy holds also for bigger values 
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m,,(GeV) m^,(GeV) 



Figure 5.7: Allowed values in the nih^ vs. nih,^ space by eqs. (5.8) and (5.9), for (5.7a) 
a = and (5.7b) a = 0.1, (5.7c) a = it/A and (5.7d) a = vr/S, for = 1 TeV 
and X = 3.5 TeV. Colours refer to different values of Q/GeV: blue (W"^), red (W'^), 
green (10^^), purple (10^^) and cyan (10^^). The shaded black region is forbidden by 
our convention m/ij > mh^ , while the shaded red region refers to the values of the scalar 
masses forbidden by LEP. 

of the angle, as the allowed region of masses is shrunk from below as we increase the 
RH Majorana neutrino Yukawa coupling, while the upper bound stays unaffected. The 
effect is even more evident for big values of the scale Q, with the Plank scale (i.e., 
Q = 10^^ GeV) precluded now for whatever Higgs boson masses at a = 7r/4 and tightly 
constraining the allowed ones at a = vr/S. 

Moving to the (nifi^-a) scan at fixed ruh^ values, figure 5.8 shows the effect of the 
heavy neutrinos in this case, to be compared to figure 5.6. It is evident that this model 
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Figure 5.8: Allowed values in the ruh. 



(d) 



VS. 



a space in the B — L model by eqs. (5.8) 



and (5.9), for (5.8a) ruh^ = 100 GeV, (5.8b) ruh^ = 120 GeV, (5.8c) rrih^ = 160 GeV 
and (5.8d) ruh^ = 180 GeV. Colours refer to different values of Q/ GeV: blue (^10"^ 
red flO^j, green ("10^°;, purple (10^^) and cyan (10^'^ ). The plots already encode our 
convention > rnh^ and the shaded red region refers to the values of a forbidden by 
LEP. Here: x = 3.5 TeV, m^^ = 1 TeV. 



can survive until very large scales Q with massive heavy neutrinos (for which, y^^ > 0.2) 
only for the light Higgs boson masses allowed in the case of the SM, that is, m^^ ~ 160 
GeV. The mixing angle must also be small, q < vr/5, providing a tight constraint on 
m/i2. For smaller h\ masses, the effect of a large y^'^ is to preclude scales Q > lO"^ GeV 
almost completely, with for example just a tiny strip for rrih-^ = 120 GeV for which there 



exists a combination of m/^j and a such that the model is consistent up to Q 



10 



10 
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GeV. Finally, figure 5.8d is not visibly different from figure 5.6d just because we are 
showing only the rrih^ > nih^ region, the shrunk region being below. 

5.1.2.2.3 VEV effect 

The last effect to evaluate comes from changing the values for the B — L breaking VEV 
X. Figure 5.9a shows the allowed regions in the m/tj vs. a plane for fixed ruh^ = 160 
GeV and y*^ = 0.2 (that is, a particular case that shows all the interesting features at 
once). As expected, since A2 is a function of rrih^/x [see for instance eq. (2.32)], at a = 
the bound on rrih^ simply scales linearly with the VEV. Regarding the upper bound, 
increasing the VEV x naively increases the allowed region for the heavy Higgs mass, 
but it is remarkable that the effects are present only for small angles, a < 0.1 radians, 
being the bigger angles unaffected. Concerning the lower bound, or the vacuum stability 
of the model, at fixed y^^ , increasing the VEV x requires to increase rrih^ to keep A2 
constant at the EW scale. This explains why, with non negligible y^^ , the allowed heavy 
Higgs masses are shrinking from below when we increase the VEV x, as one can see in 
figure 5.9a and comparing figure 5.9b with figure 5.7c, both for q = 7r/4 and y^^ = 0.2, 
but for X = 3.5 and x = 7.5 TeV, respectively. 

In general, for the model to survive up to very large scales Q ~ Mpianck, it is preferred 
the heavy neutrinos to be light with respect to the VEV x, in such a way that their 
Yukawa couplings are negligible in the RGE evolution of the scalar sector. 




Figure 5.9: Allowed values by eqs. (5.8) and (5.9) (5.9a) in the rrih^ vs. a space for 
nih^ = 160 GeV and y^' = 0.2, for Q = 10^ GeV (straight line) and Q = 10^^ GeV 
(dashed line) for several B — L breaking VEV values (x = 3.5, 7.5 and 35 TeV, giving 
m^^ = 1, 2 and 10 TeV, respectively), and (5.9b) in the m^-^ vs. m^^ space, for a = '7r/4, 
X = 7.5 TeV and y^^ = 0.2, where colours refer to different values of Q / GeV: blue (10^), 
red (lO"^ ). The plots already encode our convention nih^ > fn^-^ and the shaded red region 
refers to the values of a forbidden by LEP. 
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5.2 Phenomenology of the Higgs sector 

In the previous section we have defined the regions of the scalar sector parameter space 
that are compatible with experimental and theoretical constraints. In this section we 
summarise some of the phenomenology of the Higgs bosons in the B — L model at the 
LHC, first presented in Ref. [79]. 

The scalar Lagrangian of the B — L model is part of a rather general family of 
extensions of the SM, in which the scalar content is augmented with a scalar singlet 
(see, e.g., Refs. [19; 20; 21; 22; 23; 24; 25]). The interactions of the Higgs bosons with 
SM particles are, therefore, the same as in the traditional literature. However, a richer 
phenomenology arises in the B — L model because of the interplay with the gauge and 
the fermion sectors, where the new particles, the Z' boson and the heavy neutrinos, 
interact with the scalar bosons. The interplay of the sectors can be strong, and, in some 
region of the parameter space, both the Z' boson and the heavy neutrinos can be lighter 
than the Higgs bosons, allowing for new decays, beside those into SM particles. 

Concerning the strength of Higgs interactions, some of the salient phenomenological 
features can be summarised as follows: 

- SM-like interactions scale with cosQ;(sinQ;) for /ii(/i2); 

- those involving the other new B — L fields, like Z' boson and heavy neutrinos, 
scale with the complementary angle, i.e., with sinQ;(cosa) for hi{h2); 

- triple (and quadruple) Higgs couplings are possible and can induce resonant be- 
haviours, so that, e.g., the /12 — )• hi hi decay can become dominant if 171^2 > 2m/,^. 

We first present production cross sections for the two Higgs bosons at the LHC, 
for CM energies of ^/s = 7 and 14 TeV, as well as their BRs, for some fixed values of 
the scalar mixing angle a. Its values have been chosen in each plot to highlight some 
relevant phenomenological aspects, such as the decay into the new B — L particles. 

It turns out that the most efficient production mechanisms at the LHC are still the 
SM ones. However, new signatures arise in the decay processes. Section 5.2.3 will bring 
pieces together and present event rates for some phenomenologically viable signatures, 
as, among those relevant here, four lepton decays of a heavy Higgs boson via pairs of Z' 
gauge bosons (which, e.g., in the SM also occurs via ZZ but in very different kinematic 
regions) and heavy neutrino pair production via a light Higgs boson (yielding, e.g., the 
same multi-lepton signatures discussed in section 4.2.2 for the Z' boson). 

5.2.1 Production cross sections and decay properties 

In figure 5.11 we present the cross sections for the most relevant production mechanisms, 
i.e., the usual SM processes such as gluon-gluon fusion, vector-boson fusion, ti associated 
production, and Higgs-strahlung (whose Feynman diagrams are in figure 5.10). For 
reference, we show in dashed lines the SM case (only for hi), that corresponds to a = 0. 
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Comparing figure 5.11c to figure 5.11a, there is a factor 3-^4 enhancement passing 
from ^/s = 7 TeV to ^/s = U TeV CM energy at the LHC. 

The cross sections are a smooth function of the mixing angle a, so as expected every 
subchannel has a cross section that scales with cosa(sinQ), respectively, for /ii(/i2). As 
a general rule, the cross section for hi at an angle a is equal to that one of /12 for 7r/2 — a. 
In particular, the maximum cross section for /i2 (i.e., when a = vr/2) coincides with the 
one of hi for a = 0. 

We notice that these results are in agreement with the ones that have been discussed 
in [20; 23; 25] in the context of a scalar singlet extension of the SM, having the latter 
the same Higgs production phenomenology. Moreover, as already shown in Ref. [20], 
also in the minimal B — L context a high value of the mixing angle could lead to 
important consequences for Higgs boson discovery at the LHC: a sort of rudimentary 
see-saw mechanism could suppress hi production below an observable rate at ^/s = 7 
TeV and favour just heavy Higgs boson production, with peculiar final states clearly 
beyond the SM, or even hide the production of both (if no more than 1 fb~^ of data is 
accumulated). Instead, at ^/s = 14 TeV we expect that at least one Higgs boson will be 
observed, either the light one or the heavy one, or indeed both, thus shedding light on 
the scalar sector of the B — L extension of the SM. The region of the parameter space 
that would allow the scalar sector to be completely hidden, for example for a ~ 7r/2 
and m/i2 heavy enough to not be produced, whatever the value of m/^^ , is experimentally 
excluded by precision analyses at LEP [2]. 




Figure 5.10: Feynman diagrams for the dominant SM Higgs boson production mech- 
anisms in hadronic collisions. From top-left, clockwise: Higgs-strahlung, vector boson 
fusion, it associated production, and gluon-gluon fusion (V = W,Z;H = hi,h2)- 
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Figure 5.11: Cross sections in the B — L model for hi at the LHC (5.11a) at y/s = 7 
TeV and (5.11c) at ^ = U TeV, and for /12 (5.11b) at yfs = 1 TeV and (5. lid) at 
y/s = 14 TeV. The dashed lines in figs. (5.11a) and (5.11c) refer to a = 0. The dotted 
part of the lines in figs. (5.11b) and (5. lid) refer to /i2 masses excluded by unitarity (see 
Ref [3]). 



5.2.1.1 New production mechanisms 

All the new particles in the B — L model interact with the scalar sector, so novel 
production mechanisms can arise considering the exchange of new intermediate states. 
Among the new production mechanisms, the associated production of the scalar boson 
with the Z' boson and the decay of a heavy neutrino into a Higgs boson are certainly 
the most promising, depending on the specific masses. The Feynman diagrams for these 
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processes are in figure 5.12. Figures 5.13 and 5.14 show the cross sections for these new 
production mechanisms, for ^/s = 14 TeV and several values of a. 




Figure 5.12: Feynman diagrams for the new Higgs boson production mechanisms in 
hadronic collisions in the B — L model. From left to right: the associated production of 
a Higgs boson and the Z' boson and Higgs production via heavy neutrino (H = hi,h2)- 
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Figure 5.13: Cross sections in the B — L model for the associated production with the 
Z'^_i^ boson (5.13a) of hi at a = tt/4 and (5.13b) of h2 at a = 0, at ^/s = 14 TeV. 

Figures 5.13a and 5.13b show the cross sections for the associated production with 
the Z' boson of hi and of /i2, respectively, for several combinations of Z' boson masses 
and g'l couplings. The process is 

qq^Z'*^Z'hn2), (5.10) 

and it is dominated by the Z' boson's production cross sections (see section 3.2.1). 
Although never dominant (always below 1 fb), this channel is the only viable mechanism 
to produce /i2 in the decoupling scenario, i.e., a = 0. 

In figure 5.14 we plot the cross sections for the other new production mechanism (the 
Higgs production via heavy neutrino) against the light Higgs mass, for several choices of 
parameters (as explicitly indicated in the labels). We superimposed the red-shadowed 
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Figure 5.14: Cross sections in the B — L model for the associated production of hi with 
one heavy and one light neutrino at y/s = 14 Te V. The red shading is the region excluded 
by LEP constraints (see section 5.1.1). 

region in order to avoid any mass-angle combination that has been aheady excluded by 
LEP constraints, as discussed in section 5.1.1. The whole process chain is 

qq ^ Z' ^ Vhh'h ^ i'hnhi{2) , (5-11) 

and it requires to pair produce heavy neutrinos, again via the Z' boson (see section 4.2.1). 
Although rather involved, this mechanism has the advantage that the whole decay chain 
can be of on-shell particles, beside the peculiar final state of a Higgs boson and a heavy 
neutrino. For a choice of the parameters that roughly maximises this mechanism (M^/ = 
900 GeV, g[ = 0.13 and m^^ = 200 GeV, from figure 4.2b), figure 5.14 shows that the 
cross sections for the production of the light Higgs boson (when only one generation of 
heavy neutrinos is considered) are above 10 fb for m/i^ < 130 GeV (and small values of 
a), dropping steeply when the light Higgs boson mass approaches the kinematical limit 
for the heavy neutrino to decay into it. Assuming the transformation a — ?• it/2 — a., the 
production of the heavy Higgs boson via this mechanism shows analogous features. 

These new production mechanisms have rather small cross sections at the LHC. 
Nonetheless, a future LC could be the suitable framework to probe them, as shown in 
Ref. [91]. We showed therein that these mechanisms can also be the leading ones for 
accessing the scalar sector of the B — L model. 
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5.2.2 Branching ratios and total widths 

Moving to the Higgs boson decays, figure 5.15 shows the BRs for both the Higgs bosons, 
hi and /12, respectively. Only the 2-body decay channels are shown here. For a descrip- 
tion of the 3-body decays, see Ref. [79]. 
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Figure 5.15: (5.15a) BRs for hi for a = 27r/5 and m^i^ = 50 GeV and (5.15c) hi total 
width for a choice of mixing angles and (5.15b) BRs for h2 for a = 37r/20 and nifi^ =120 
GeV, Mz' = 210 GeV and m^^ = 200 GeV and (5.15d) /12 total width for a choice of 
mixing angles. 



Regarding the light Higgs boson, the only new particle it can decay into is the heavy 
neutrino (we consider a very light Z' boson unlikely and unnatural), if the channel is 
kinematically open. In figure 5.15a we show this case, for a small heavy neutrino mass. 
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i.e., niu^ = 50 GeV, and we see that the relative BR of this channel can be rather 
important, as the decay into b quark pairs or into W boson pairs, in the range of masses 
110 GeV < rrihi — GeV. Such a range happens to be critical in the SM since here the 
SM Higgs boson passes from decaying dominantly into b quark pairs to a region in masses 
in which the decay into W boson pairs is the prevailing one. These two decay channels 
have completely different signatures and discovery methods/powers. The fact that the 
signal of the Higgs boson decaying into b quark pairs is many orders of magnitude below 
the natural QCD background, spoils its sensitivity. In the case of the B — L model, the 
decay into heavy neutrino pairs is therefore phenomenologically very important, besides 
being an interesting feature of the B — L model if rriy^ < Mw, as it allows for multi- 
lepton signatures of the light Higgs boson. Among them, there is the decay of the Higgs 
boson into 3£, 2j and (studied for the Z' case in section 4.2.3), into 4£ and (as, 
again, already studied for the Z' case in Ref. [60]) or into M and 2j (as already studied, 
when £ = /X, in the 4*^ family extension of the SM [74]). All these peculiar signatures 
allow the Higgs boson signal to be studied in channels much cleaner than the decay into 
b quark pairs. 

In the case of the heavy Higgs boson, further decay channels are possible in the 
B — L model, if kinematically open. The heavy Higgs boson can decay in pairs of the 
light Higgs boson (/12 — ?• hi hi) or even in triplets (/i2 — )• hi hi hi), in pairs of heavy 
neutrinos and Z' bosons. Even for a small value of the angle, figure 5.15b shows that 
the decay of a heavy Higgs boson into pairs of the light one can be quite sizeable, at the 
level of the decay into SM Z bosons for rrih^ = 120 GeV. 

The BRs of the heavy Higgs boson decaying into Z' boson pairs and heavy neutrino 
pairs decrease as the mixing angle increases, getting to their maxima (comparable to 
the W and Z ones) for a vanishing a, for which the production cross section is however 
negligible. As usual, and also clear from figure 5.15b, the decay of the heavy Higgs boson 
into gauge bosons (the Z' boson) is always bigger than the decay into pairs of fermions 
(the heavy neutrinos, even when summed over the generations, as plotted), when they 
have comparable masses (here, Mz' = 210 GeV and m^^ = 200 GeV). It is important to 
note that all these new channels do not have a simple dependence on the mixing angle a 
(see figure 5.16 for Z' bosons and neutrinos final states and Ref. [79] for the light Higgs 
boson case). 

The other standard decays of both the light and the heavy Higgs bosons are not 
modified substantially in the B — L model (i.e., the Higgs boson to W boson pairs is 
always dominant when kinematically open, otherwise the decay into b quarks is the pre- 
vailing one; further, radiative decays, such as Higgs boson decays into pairs of photons, 
peak at around 120 GeV, etc.). Only when other new channels open, the standard de- 
cay channels alter accordingly. This rather common picture could be altered when the 
mixing angle a approaches '7r/2, but such situation is phenomenologically not viable [2]. 

Figures 5.15c and 5.15d show the total widths for hi and /i2, respectively. In the first 
case, few thresholds are clearly recognisable as the heavy neutrino one at 100 GeV (for 
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angles very close to 7r/2 only), the W and the Z ones. Over the mass range considered 
(90 GeV < m/ij < 250 GeV), the particle's width is very small until the W threshold, 
less than 1 — 10 MeV, rising steeply to few GeV for higher hi masses and small angles 
(i.e., for a SM-like light Higgs boson). As we increase the mixing angle, the couplings of 
the light Higgs boson to SM particles is reduced, like its total width. 

On the contrary, as we increase a, the /12 total width increases, as clear from fig- 
ure 5.15d. Also in this case, few thresholds are recognisable, as the usual W and Z gauge 
boson ones, the light Higgs boson one (at 240 GeV) and the t quark one (only for big 
angles, i.e., when /12 is the SM-like Higgs boson). When the mixing angle is small, the 
/i2 total width stays below 1 GeV all the way up to m/jj ~ 300 500 GeV, rising as the 
mass increases towards values for which T^^ ~ nifi^ ~ 1 TeV and /12 loses the meaning 
of bound state, only for angles very close to 7r/2. Instead, if the angle is small, i.e., less 
than 7r/10, the ratio of width over mass is less than 10% and the heavy Higgs boson is 
a well-defined particle. In the decoupling regime, i.e., when a = 0, the only particles /12 
couples to are the Z' boson and the heavy neutrinos. The width is therefore dominated 
by the decay into pairs of them and it is tiny, as is clear from figure 5.15d. 
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Figure 5.16: Dependence on the mixing angle a of (5.16a) BR(h2 — )• I'h^h) <ind of 
(5.16b) BR{h2 Z'Z'). 

As already mentioned, figure 5.16 shows the dependence on the mixing angle a 
of the BRs of /i2 into pairs of Z' bosns and heavy neutrinos, /12 — s- I'h and /12 — )■ 
Z' Z' , respectively (not influenced by m^J. As discussed, the interaction of the heavy 
Higgs boson with SM (or non-SM) particles has an overall sin a (or cos a, respectively) 
dependence. Nonetheless, the BRs in figure 5.16 depend also on the total width, that 
for a > 7r/4 is dominated by the /i2 — )• W^W~ decay. Hence, when the angle assumes 
big values, the angle dependence of the /12 BRs into heavy neutrino pairs and into Z' 
boson pairs follows a cot a behaviour. For a study of BR(/i2 — )■ /ii /ii), see Ref. [79]. 
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5.2.3 Cross sections for processes involving Z' bosons and neutrinos 

In this subsection we combine the results from the Higgs boson cross sections and those 
from the BR analysis in order to perform a detailed study of typical event rates for some 
Higgs signatures that are specific to the B — L model. 

As in section 3.2, we identify two different experimental scenarios related to the LHC: 
the 'early discovery scenario' (i.e., with CM energy of ^/s = 7 TeV and an integrated 
luminosity of L = 1 fb~^), and the 'full luminosity scenario' (i.e., with CM energy of 
^/s = 14 TeV and an integrated luminosity of L = 300 fb~^). 

Combining the production cross sections and the decay BRs presented in the previous 
subsections, the two different scenarios open different possibilities for the detection of 
peculiar signatures of the model, involving heavy neutrinos and the Z' boson. In the 
early discovery scenario there is a clear possibility to detect a light Higgs state yielding 
heavy neutrino pairs, while the full luminosity scenario affords the possibility of looking 
for the heavy Higgs state decaying into Z' boson pairs, as well as into heavy neutrino 
pairs. In addition to them, the decays of the heavy Higgs state into light Higgs boson 
pairs are possible, as discussed in Ref [79]. 

The Feynman diagrams for these processes are shown in figure 5.17, with the gluon- 
gluon fusion being the dominant production mechanism. 



Figure 5.17: Feynman diagrams for the new Higgs boson discovery mechanisms in 
hadronic collisions in the B — L model. From left to right: the decay of a Higgs bo- 
son into heavy neutrino pairs and into Z' boson pairs. 

First, we focus on the early discovery scenario: in this experimental configuration, 
the most important B — L distinctive process is the heavy neutrino pair production via 
the light Higgs boson: 



In figure 5.18 we show the cross sections, in the {vrih^-a) plane, for this process at the 
LHC for = 7 TeV, for two different heavy neutrino masses {m^^^ = 50, 60 GeV). 
These plots are the result of the combination of the light Higgs boson production cross 
sections via gluon-gluon fusion only (that represents the main contribution) and of the 
BR of the light Higgs boson to heavy neutrino pairs. 

At this energy configuration, the expected integrated luminosity is not above 1 fb~^, 
yielding the production of up to around 200 (or more) heavy neutrino pairs via the light 
Higgs boson, for m^^^ = 50 GeV, 110 GeV< nih^ < 150 GeV and tt/A < a < O.AItt. 
For m,^^ = 60 GeV, up to 100 heavy neutrino pairs can be produced, for 130 GeV< 



9 



9 





(5.12) 
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nihi < ISO GeV, and tt/4 < a < 0.457r. Notice that the maximum cross section [of 
about 250(100) fb, for ttIj,^ = 50(60) GeV] occurs for a mixing angle of around vr/S, due 
to the interplay between the production mechanism (maximised for a — )■ 0) and the BR 
into heavy neutrino pairs (maximised for a — ?■ 7r/2), with the latter being the dominant 
part. 
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(a) (b) 
Figure 5.18: Cross section times BR contour plot for the process pp ^ hi ^ ^^h^h o,t 
the LHC for ^/s = 7 TeV, plotted in the mh^-a plane, with (5.18a) m^^ =50 GeV and 
(5.18b) ruy^ = 60 GeV. The red-shadowed region is excluded by the LEP experiments 

[!]■ 



In the full luminosity scenario, several important distinctive signatures appear. We 
focus here on the Z' boson and heavy neutrino pair production via /i2, 

pp^ /i2 ^Uh^h, (5.13) 
pp^ h2 -^Z'Z'. (5.14) 

In analogy with the previous case, figure 5.19 shows the cross sections for the process 
of eq. (5.13) (via the heavy Higgs boson), at the LHC for y/s = 14 TeV, for ?n,y^ = 150 
GeV and ruy^ = 200 GeV. In this case, as expected, the maximum cross section is for the 
complementary angle, i.e., for about tt/6. The expected high integrated luminosity for 
this energy stage of the LHC is such that a good amount of events can be produced even 
for cross sections of around 1 fb, so that for both choices of the heavy neutrino mass, 
there exist values of a and m^^ for which the event rate could lead to the observation 
of this signature. In particular, for tti^^,^ = 150 GeV, we find a cross section times BR of 
0.85 fb (that corresponds to around 250 events, at most) for 320 GeV< m/jj < 520 GeV 
and 0.037r < a < 0.337r. For m,y^ = 200 GeV, the same value of the cross section is for 
450 GeV< lUh^ < 550 GeV and 0.037r < a < 0.2l7r. 
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(a) (b) 
Figure 5.19: Cross section times BR contour plot for the process — )■ /12 — )• I'h'^h o-i the 
LHC for y/s = 14 TeV, plotted in the (m/ij-aj plane, with (5.19a) m^^ = 150 GeV and 
(5.19b) m,^^ = 200 GeV. The red-shadowed region is excluded by unitarity constraints 
(see section 5.1.2.1). 
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Figure 5.20: Cross section times BR contour plot for the process — t- /12 — Z' Z' at the 
LHC for y/s = 14 TeV, plotted in the (rrih^-a) plane, (5.20a) with mz' = 210 GeV, and 
(5.20b) mz' = 280 GeV. The red-shadowed region is excluded by unitarity constraints 
(see section 5.1.2.1). 



In figure 5.20 we show the cross sections, in the (01^.^-0) plane, for the process of 
eq. (5.14) (i.e., the Z' boson pair production via the heavy Higgs boson) at the LHC for 
y/s = 14 TeV, for mz' = 210 GeV and mz' = 280 GeV. This process is also favoured for 
small mixing angles, with a maximum cross section of 0.85 fb for 420 GeV< < 650 
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GeV and O.OSvr < a < 0.257r, for the former Z' mass. For a bigger value of the Z' boson 
mass, such as Mz' = 280 GeV, a noticeable parameter space is still potentially accessible, 
with a maximum cross section of 0.3 fb (corresponding to around 100 events, at most) 
for 560 GeV< m/jj < 800 GeV and O.OSvr < a < 0.197r. The Z' bosons can subsequently 
decay into pairs of leptons, both electrons and muons, with BR(Z' — t- ~ 30% 

£=e,fi 

(see figure 3.4). Thus, in total, up to around 25 events of heavy Higgs bosons decaying 
into 4 leptons could be observed, in a completely diff^erent kinematical region from the 
decay into SM Z bosons. 

5.3 Conclusion 

In this chapter, we have presented the results of our investigation of the scalar sector, 
with focus on the interaction with the gauge and fermion sectors. 

First, we have presented a summary of the existing experimental constraints and the 
full analysis of the high energy theoretical constraints on the scalar sector. 

We have briefly summarised the unitarity bound, whose scope was to clarify the 
role of the two Higgs bosons in the unitarisation of vector and scalar bosons scattering 
amplitudes, that we know must hold at any energy scale. 

Using the equivalence theorem, we have evaluated the spherical partial wave ampli- 
tude of all possible two-to-two scatterings in the scalar Lagrangian at an infinite energy, 
identifying the zz — )• zz and z'z' — )• z'z' processes as the most relevant scattering chan- 
nels for this analysis {z^ ^ is the would-be Goldstone boson of the Z^ ^ vector boson). 

Hence, we have shown that these two channels impose an upper bound on the two 
Higgs masses: the light one cannot exceed the SM bound while the limit on the heavy 
one is established by the singlet Higgs VEV, whose value is presently constrained by 
LEP and may shortly be extracted by experiments following the discovery of the Z'q_j^ 
boson. 

We have then investigated the triviality and vacuum stability conditions of the pure 
B — L model with a particular view to define the phenomenologically viable regions of 
the parameter space of the scalar sector, by computing the complete set of RGEs (gauge, 
scalar and fermion) at the one-loop level in presence of all available experimental con- 
straints. The RGE dependence on the Higgs masses and couplings (including mixings) 
has been studied in detail for selected heavy neutrino masses and couplings as well as 
discrete choices of the singlet Higgs VEV. 

Altogether, we have found that there exist configurations of the model for which 
its validity is guaranteed up to energy scales well beyond those reachable at the LHC 
while at the same time enabling the CERN hadron collider to probe its scalar sector in 
Higgs mass and coupling regions completely different from those accessible to the SM. 
Furthermore, we have shown that investigations of the Higgs sector of this extended 
scenario may also lead to constraints on other areas, such as the (heavy) neutrino and 
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Z' sectors (the latter indirectly, through the VEV of the singlet Higgs state directly 
intervening in the scalar RGEs). 

Finally, we have investigated the production and decay phenomenology of both Higgs 
states of the minimal B — L model at the LHC, at both the foreseen energy stages of 
the LHC (i.e., y/s = 7 and 14 TeV) and with the corresponding integrated luminosities. 
While virtually all relevant production and decay processes of the two Higgs states of 
the model have been investigated, we have eventually paid particular attention to those 
that are peculiar to the described B — L scenario. The phenomenological analysis has 
been carried out in the presence of all available theoretical and experimental constraints 
and by exploiting numerical programs at the parton level. While many Higgs signatures 
already existing in the SM could be replicated in the case of \ts B — L version, in either 
of the two Higgs states of the latter (depending on their mixing), it is more important 
to notice that several novel Higgs processes could act as hallmarks of the minimal B — L 
model. These include Higgs production via gluon-gluon fusion, of either the light or 
heavy Higgs state, the former produced at the lower energy stage of the CERN collider 
and decaying in two heavy neutrinos, and the latter produced at the higher energy stage 
of such a machine and decaying not only in heavy neutrino pairs but also in Z' boson 
pairs (and eventually in light Higgs boson pairs). For each of these signatures we have 
in fact found parameter space regions where the event rates are sizable and potentially 
amenable to discovery. 

While, clearly, detailed signal-to-background analyses will have to either confirm or 
disprove the possibility of the latter, our results have laid the basis for the phenomeno- 
logical exploitation of the Higgs sector of the minimal B — L model at the LHC. 



Chapter 6 



Discussion 



This Thesis studies the phenomenology of a minimal extension of the SM, namely, the 
B — L model. This well-motivated framework naturally implements the see-saw mech- 
anism for neutrino mass generation, being the RH neutrinos required for anomaly can- 
cellation. It also explains the accidental and non-anomalous U{\)b-l global symmetry 
in the SM, being this, here, the remnant after spontaneous symmetry breaking of its 
local counterpart. It finally introduces scope in the Higgs boson sector, that lacks any 
experimental observation. 

We think that the goal of a proficient interaction with experimentalists has been 
achieved. A joint effort has been fundamental in the studies concerning the signatures 
of the model and the strategy for their observation. Also, some of the gaps in the 
preparation towards real data have been filled in this work. The model under discussion 
has been studied in many of its aspects, identifying novel and exciting signatures. 

The phenomenology of the new particles of the model has been analysed, possibly, 
in the most comprehensive way. The peculiarities with respect to the SM, and to some 
of its extensions, have been described and highlighted, focusing on the most interesting 
signatures one may observe. 

In the attempt of a complete study of the model, our work has defined the allowed 
regions of its parameter space and presented a detailed study of some of its signatures, 
in the gauge, fermion and scalar sectors, at the LHC, waiting for the data to confirm or 
disprove this TeV scale realisation of the minimal B — L model. 

6.1 Summary of the results 

We summarise here the most important results we achieved in this Thesis. For a more de- 
tailed description, see the partial conclusions for each chapter, in section 3.3, section 4.3, 
and section 5.3, for the gauge, fermion, and scalar sector, respectively. 

The non-observation of any Z' boson nor of Higgs boson at past and recent colliders 
results in lower bounds on the parameter space of the model. We have presented the most 
up-to-date constraints for the former at Tevatron in section 3.1.1.1. Complementary to 
it, theoretical arguments can mainly give upper bounds on the free parameters. The 
topics we studied in this Thesis are the unitarity bound of the scalar sector (that we 
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summarised in section 5.1.2.1) and the constraints coming from the analysis of the 
RGEs of the model at one loop (that we calculated and collected in appendix B). The 
requirement that the evolution of the parameters does not hit any Landau pole, up to 
a specific cut-off scale, is the so-called triviality bound, that we studied for the gauge 
couplings in section 3.1.2, and for the scalar parameters in section 5.1.2.2. In the latter, 
this analysis has been combined with the study of the vacuum stability condition, leading 
also to lower bounds on the Higgs boson masses. It is remarkable that there exist 
configurations of the model for which its validity is guaranteed up to very high energy 
scales, well beyond those reachable at the LHC, that are completely different form the 
SM ones. In particular, a hi boson lighter than the SM limit of 114.4 GeV is allowed 
for a > tt/3 and a suitable choice of the other parameters. 

We have then investigated the collider signatures of this model at the LHC, focusing 
on the gauge and fermion sectors. In chapter 3 we have shown that the Z'^_^ boson 
is dominantly coupled to leptons, whose total BR sums up to 75%. Also, it can decay 
into heavy neutrino pairs, with a BR up to around 20%. We finally presented (in 
section 3.2.2) the parton level discovery potential of the Z'^^i^ boson at Tevatron and the 
LHC in Drell-Yan processes. The LHC at y/s = 7 TeV will be able to give similar results 
to Tevatron, in the electron channel and in the low mass region, and is able to extend 
its reach for heavier Z' boson, up to masses of 1.25(1.20) TeV for electrons(muons). 
When the data from the high energy runs at the LHC become available, the discovery 
reach of the Z'^_j^ boson will be extended towards very high masses and small couplings 
in regions of parameter space well beyond the reach of Tevatron. We also have shown 
that greater sensitivity to the Z'q_^ resonance is provided by the electron channel, that 
at the LHC has better energy resolution than the muon channel. This in turn could 
enable the estimate of the gauge boson width. Finally, we have commented that the 
inclusion of further background, as well as a realistic detector simulation, will not have 
a considerable impact on the results that we have presented. 

The phenomenology of heavy neutrinos has been presented in chapter 4. Due to the 
see-saw mechanism, that allows for a tiny LH component, the main decay modes for the 
heavy neutrinos (i.e., the decay into SM gauge bosons) have very small partial widths. 
Hence, in a large portion of the parameter space, the heavy neutrino is a long-lived 
particle, with a proper lifetime that can be comparable to (or exceed) the b quark's 
one (for example, they are equal for mi,^ = 200 GeV and m^^ = 10~^ eV). We have 
emphasised that an independent measurement of the heavy neutrino lifetime (through 
the displacement of its decay vertex) and its mass may lead to infer the absolute mass 
of the SM-like neutrinos. Therefore, we have proposed a way to measure the heavy 
neutrino mass. We have described that heavy neutrino pairs, produced by the Z' boson, 
give rise to multi-lepton decays of the latter and we have identified the tri-lepton one 
as the most interesting one for a detailed study. We have proposed cuts and a search 
strategy, identifying the suitable distributions one should look at. A parton level study, 
in section 4.2.3.1, has shown that simple cuts can effectively reduce the backgrounds, 
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that the peak corresponding to the heavy neutrino enables a rather precise measure of 
its mass and that a striking signature for this model is the simultaneous presence of 
peaks for the intermediate Z' boson and for the heavy neutrinos. Finally, a full detector 
simulation, in section 4.2.3.2, has validated the results of the parton level study in a 
more realistic framework, showing that a discovery could be possible with less than 100 
fb-i. 

Finally, in chapter 5, we have laid the basis for the investigation of the Higgs sector. 
The most appealing result is that multi-lepton decays of a Higgs boson, via heavy neu- 
trino pairs, are at the same time rather peculiar and they might allow for the discovery 
of the intermediate scalar in the early stage of the LHC. 

We can reasonably conclude that the heavy neutrinos truly carry the hallmarks of 
the B — L model at colliders. 

6.2 Outlook on future work 

In this Thesis, we have mentioned several topics that were left for future investigations, 
primarily for time limitations. We further discuss here some of them as a possible 
continuation of the work undertaken. Also, further developments and new studies are 
addressed. 

When discussing the phenomenology of the Z' boson, in chapter 3, we have discussed 
that the asymmetries of the decay products stemming from the Z'^_^ are trivial at the 
peak. Hence, we have decided to not study them. However, as already pointed out, 
asymmetries become important in the interference region, especially just before the 
Z' boson peak, where the Z — Z' interference will effectively provide an asymmetric 
distribution somewhat milder than the case in which there is no Z' boson. A first 
attempt to study them is in Ref. [92], though further improvements may be possible, for 
example by studying the asymmetries of the heavy neutrinos, which, as we have shown, 
have just axial couplings to the gauge boson. 

Another important aspect that needs to be quantitatively addressed is the impact on 
the model of the electroweak precision observables. Although the impact on the gauge 
sector has been found to be negligible in the pure B — L model, at one loop (see, e.g., 
[51]), one is left with the scalar and fermion sectors. The former has been studied in 
Ref. [2], though without considering heavy neutrinos, that, as we highlighted, might 
contribute. 

In chapter 4, we have highlighted that the tri-lepton signature via heavy neutrino 
pairs is one of the most distinctive signatures of the pure B — L model at colliders. 
We have successfully studied it with a full detector simulation, concluding that a rather 
precise measurement of the heavy neutrino mass is a realistic task. It might be appealing 
to improve the undertaken analysis: the implementation of the full light neutrino masses 
and mixing angles, as well as considering the heavy neutrinos as not degenerate, would 
be a further step, whose aim is to show the feasibility of extracting the mass of each 
light neutrino in a realistic case. Further, the precision of the heavy neutrino mass 
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measurement and the impact of fully decayed tau leptons ought to be quantitatively 
addressed. 

Likewise, in chapter 5, it has been pointed out that the neutrino pair production could 
be important for the light Higgs boson discovery in the pure B — L model, in a range of 
masses that happens to be critical in the SM, since here the SM Higgs boson passes from 
decaying dominantly into h quark pairs to a region in masses in which the decay into W 
boson pairs is the prevailing one. It is also important to note that the former signature 
provides poor sensibility, being it many orders of magnitude below the natural QCD 
background. Therefore, the decay into heavy neutrino pairs is phenomenologically very 
important, especially at the LHC in its early stage (i.e., for ^/s = 7 TeV). Beside being 
an interesting feature of the B — L model if rriy^ < Mw, as it allows for multi-lepton 
signatures of the light Higgs boson (some of them being much cleaner than the decay into 
b quark pairs), these decays could provide the first hints of a fundamental scalar particle. 
The analysis of the tri-lepton signature for the Higgs boson is certainly a very interesting 
task to pursue at the LHC at both foreseen CM energies. Furthermore, a full analysis of 
the kinematic distribution of the decay products is also crucial for the identification of 
the spin of the intermediate resonance. Although in a very different kinematic region, 
the Z' boson gives rise to the same decay products, hence the intermediate resonance 
has to be distinguished, also for model discrimination purposes. 

Finally, the embedding of the model in GUTs could be pursued. In this respect, 
we note that the supersymmetric version of the pure B — L model has already been 
considered: see, for example, [93; 94; 95; 96; 97; 98] and references therein. This 
implementation combines the virtues of the pure B — L model (namely, the natural and 
dynamic implementation of neutrino masses) with the elegant solution of the hierarchy 
problem of the SM, as well as naturally providing a dark matter candidate, that are 
typical for a supersymmetric framework. Because of the enlarged spectrum, it might 
be interesting to apply the collider studies of this Thesis to the supersymmetric B — L 
model, analysing the discovery potential at the LHC of the new particles and the new 
signatures arising in this enlarged framework. 



Appendix A 

Model implementation 



In this appendix is presented the implementation of the pure B — L model into the 
numerical tools that have been used. Also, the Feynman rules are collected in section A. 3. 

A.l Implementation of the pure B — L model 

Most of the numerical analysis of this Thesis ^ has been carried out using the calchep 
package [41], an automated tool for tree level calculation of physical observables. 

We implemented the pure B — L model, for which the covariant derivative given by 
eq. (2.7) reads 

D^ = d^ + igsT'^G^^ + zffT'^P^/ + ig^YB^ + ig[YB-LB'^ . (A.l) 

For its straightforward implementation in CALCHEP, we have used the lanhep module 
[42]. This package, beside providing the suitable output for calchep, also derives the 
Feynman rules of the model, that we collected in section A. 3. 

The availability of the model implementation into calchep in both the unitary 
and t'Hooft-Feynman gauges allowed us to perform powerful cross-checks to test the 
consistency of the model itself. 

The implementation of the gauge sector is quite straightforward. Since there is no 
mixing between the SM Z and the Z'g_j^ bosons, one just needs to define a new heavy 
neutral gauge boson together with the simplified covariant derivative given by eq. (A.l), 
and the charge assignments in table 2.1. For the scalar sector, the mixing between 
mass and gauge eigenstates of the two Higgs bosons, as well as the reformulation of the 
Lagrangian parameters in terms of the physical quantities, has been done accordingly 
to eqs. (2.28)^(2.33). 

Finally, regarding the fermion sector, the neutrinos ought to be carefully handled, 
as explicitly described in the following subsection. 

part from the full detector simulation of the tri-lepton signature (see section 4.2.3.2), the study 
of the RGEs (see section 5.1.2.2), done with mathematica [43], and the study of the unitarity bound 
(see section 5.L2.1). 
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A. 1.1 Neutrino eigenstates 

The implementation of the neutrino sector is somewhat more compHcated. Majorana- 
hke Yukawa terms are present in eq. (2.10) for the RH neutrinos, therefore one must 
implement this sector such that the gauge invariance of the model is explicitly preserved. 
This can be done as follows. As a first step we rewrite Dirac neutrino fields in terms of 
Majorana ones using the following general substitution: 

where is a Dirac field and i^l{R) its left(right) Majorana components. If we 
perform the substitution of eq. (A. 2) in the neutrino sector of the SM, we will have 
an equivalent theory formulated in terms of Majorana neutrinos consistent with all 
experimental constraints. 

The derivation of the mass eigenstates is the next step, and it has been explained in 
detail in section 2.2.3. 

The last subtle point is the way the Lagrangian has to be written to meet the 
requirements of lanhep. In particular, this regards the Majorana-like Yukawa terms 
for the RH neutrinos [the last term in eq. (2.10)]. In order to explicitly preserve gauge 
invariance, this term has to be written, in two-component notation, as: 

V^'^^^^^X + h.c, (A.3) 

where ly is the Dirac field of eq. (A. 2), whose Majorana components ui^ji mix as in 
eq. (2.37). 

Altogether, this implementation is an explicit gauge-invariant formulation of the 
neutrino sector suitable for the lanhep tool. The specific interactions pertaining to the 
heavy neutrinos are collected in the following subsection. 

A. 1.2 Feynman rules involving heavy neutrinos 

Given their importance, and to help the reader, we list here the Feynman rules involving 
the heavy neutrinos in the pure B — L model. The intervening quantities are defined in 
the main text. 
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2x V / 



U- Ui 

(/^sin2ajy 5^ cos 2ai/ 

where y'^ = ^- — ^— , sin2a,y = —2 — ^ , cos2ai, = — 

A. 2 One-loop vertices 

For the correct analysis of the scalar sector of chapter 5, one- loop vertices need to be 
considered. They couple the Higgs fields to massless particles, such as photon and gluon, 
for which no tree level interaction exists. In fact, such couplings are mediated by a loop of 
either (heavy) quarks or massive gauge bosons (or both). Despite the mass suppression, 
due to the integration of the particle in the loop, these interactions can still be sizable. 
In particular, the coupling to the gluons is responsible for the most effective production 
mechanism for the Higgs boson at the LHC, the so-called gluon-gluon fusion. 
A detailed description of the implementation is as follows. 

• The one-loop vertices g — g — hi{h2), 7 — 7 — hi{h2) and 7 — Z{Z') — hxih^) via W 
gauge bosons and heavy quarks (top, bottom, and charm) have been implemented, 
adapting the formulas in Ref [99]. 

• Running masses for top, bottom, and charm quarks have been considered, evalu- 
ated at the Higgs boson mass: Q = ray^^imh^ (depending on which scalar boson 
is involved in the interaction). 



• Running of the QCD coupling constant has been included, at two loops with 5 
active flavours. 
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Finally, the NLO QCD fc-factor for the gluon-gluon fusion process [100; 101; 102] ^ has 
been used. Regarding the other processes, we decided to not implement their fe-factors 
since they are much smaller in comparison. 

It shall be noticed that the implementation of the formulas for the running quark 
masses and for the running QCD coupling constant, as well as the fc-factor for the gluon- 
gluon fusion, have been done for the consistency of the NLO evaluation and to get in 
agreement with the cross sections in the literature (see, e.g., Ref. [102]). 

In the following subsection we list the complete set of tree level Feynman rules 
for the pure B — L model, as given by lanhep. Where possible, the couplings have 
been expanded upon the mass eigenstate basis. The one-loop vertices are not included. 
We remand the reader to the formulas in Ref. [99] for the SM, of which ours are a 
straightforward extension. 

Concerning the notation we use, the following applies. 

• A is the photon, G is the gluon. 

• Vp is the Goldstone field of the vector V . 

• is the ghost field related to the vector V , and is the corresponding conju- 
gated field. 

• Pn is the four-momentum of the n!'^ vector field (as it appears in the left column). 

• sw = sinvK is the sine of the EW angle (and cw the cosine). 

• Sa = sin a is the sine of the scalar mixing angle a (and Cq the cosine). 

• = sina^ is the sine of the neutrino mixing angle Uy (and Ca,, the cosine). 

^Notice that in Ref. [101] (Ref. [102]), mt = 174(178) GeV, while we used rut = 172.5 GeV as top 
quark pole mass value. 
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1 e-Mc-Sa X X , 

2 Mw-Su, "pi ""^ 
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A. MODEL IMPLEMENTATION 



Fields in 


the vertex 


Variational derivative of Lagranj 


^ian by fields 










l^)cb - ■ (1 + -f^)a) 


_ 


Cbq 


Zp 


1 i-e-Mc x 5 
2MwSn, PI '^ab 




- 


Cha 


z ^ 


1 / r p r 

~ "^Ql^pqlac ■ Ocb 




_ 


dbq 


TT7-I- 








dbq 


TTT-I- 


li-e-M^-V2-V,d /I ^.5\ r 
4 Mw-s^u ^ ' )abOpq 




C(ip 


Sba 

uq 


TT7'+ 








Sba 

uq 




l^l^^^li^s ■ (1 + l')ab - M, . (1 - 7^),,) 




Cbq 


w ^ 






-J 


^bq 




1 i-e-Mc-V2-V^d (-1 1 -,5^ r 
4 Mw'Sw ^ ' ^ 




^ap 


^oq 


\ 


-^G5pq^ac ' ^cb 






dhri 
oq 




9s ' InhKa 
-"^ 'do p<j 




dap 


dbq 




-T27^^9<c{2sJ-il + j%b 


- (3 - 2s J) ■ (1 - j^),b) 




dbn 
uq 


Z p. 


-la'l^pqlac ■ 5cb 




dfnp 


tba 




4 Syj ^ 1 J <^u pq /(^L- 






tba 
uq 


Wp 






"'ap 


111. 








Ea 


Cb 




e^ac ■ Scb 




P ^ 


Pu 
^b 


1 

hi 


1 Ca-e-Me X 

2 Mw-Sn, "''^ 




Ca 


Pb 


h2 


1 e-Ms-Sa X , 




Ea 


Pb 


Z^, 


i^7a^c((l-2..2).(l-7^U 


-2sJ-{l + ^%b) 




Pu 
^b 


Zp 


1 i-e-Me „5 

2 Mtjv-Sw lab 




ea 






g'llac ■ Scb 




ea 


Vlb 








ea 


l^lb 






2s^Sc,^Mwy^,il + l^)ab) 


ea 


^hlb 


w-p 
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Fields in 


the vertex 


Variational derivative of Lagrangian by fields 


ea 


^hlb 




i M^^KeMe V2(l - 7^)a6 + 2s^Ca^MwyS,{l + 7%b) 


G^p 


Guq G pf 




C^p 


G'^q Gf^r 


gs ' P2fpqr 




hi 


hi 


-l-^:^{2cJ'Sy,CaMh2'^Mwx - 2cc^s^C2aCaMh^'^Mwx 

+2Ca SwCaMh/ MwX + 2Ca SwC2aCaMhj^^ MwX 

-sJeMh^^v'^ + sJc2aeMh^^v'^ - sJeMh^^v^ - sJc2aeMh^^v'^ 

-Cq^ SaeMh^'^ S2aVX + 2Su,Sa^ CaM^^'^ Mw S2aV 
+Ca^SaeMh^'^S2aVX - 2SwSa^ CaMh^"^ Mw S2av) 


hi 


hi 


h2 


-\7^{^Ca^SwSaMh^'^MwX - QCa^SwSaC2aMh^'^MwX 
+QCq^ SyjSaMh^^^ My/X + QcJ^ SwSaC2aMh^^ My/X 

„9 71^99 .-.9 -iij'99 n Ti^99 

+ZSa CaeMh^'^V^ - 3Sa C2aCaeMh/v'' + 3Sa CaeMh^ V 
+3Sa^ C2aCaeMh2^v'^ + (1 - 'iSa^)CaeMh2^ S2aVX 
-2(2 - 2,Sa^)SyjSaMh^'^ MwS2aV - (1 - 2,Sa^)CaeMh^^ S2aVX 
+2(2 - ?,sJ)s^So,Mh^'^MwS2aV) 


hi 


h2 


h2 


-\7^{^Su,Sa^CaMh^'^MwX - QSy,Sc? C2aCaMh^^ My/X 
+QSu,Sa^CaMh^^MY/X + Qs^Sa^ C2aCaMh^'^ My/X 

o Ti4'99 t-i9 Ti/r99 f-i9 Ti^'99 

— SCc^So-eM/jj t)^ + SaC2aeMh^'^V^ — 3Ca SaeMh2 V 

o9 7i^'99 /r\ o9\ 7i/r9 
— OCq, SQ,C2QeiW/j2 U + (2 — 6Sa )SaeMfi2 S2aVX 

+2(1 - 3Sa'^)s^CaMh2'^MwS2aV - (2 - 3Sa2)s„eM,j^ ^^^aVX 

-2(1 - 3s„^)s^CaM,j/Mvi.S2a't;) 


hi 


W+ 




Ca-e-Mw uu 

Sw ^ 


hi 


W+ 






hi 


w+ 






hi 


w+ 




+2SwC2aCaMh^'^ Mw - SaeMh^"^ S2aV + SaSUhj^^ S2av) 
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A. MODEL IMPLEMENTATION 



Fields in 


the vertex 


Variational derivative of Lagrangian by fields 




Z, 










Zf 


2c^-s^\P3 Pi) 


' ''1 




Zf 


-l^{2s^,CaMh^^Mw - 2s^C2aCaMh^'^Mw + 2s^CaMh,'^Mw 








+2SyjC2aCaMhi'^ Mw - SaeMh2^S2aV + SaeMh^'^ S2av) 


' ''1 


z', 




-8sa -g'l^-x- g^" 




z', 


Z'f 


-2t-sa-gi [Pi-Ps) 


' ''1 


z' 


Z'f 


11/ 7ir9 71^9 ■7iyr9 

2Firs('5a^^^/ii ^ ~ SaC2aeMh/v + S^eMh^ V 








+SaC2aeMh^'^V - 2SwCaMh^'^ Mw S2a + 2SwCaMh^'^ Mw S2a) 


ho 


ho 


h2 


-^-^::^{'^Su,Sa^Mh2^Mwx - 2swsJc2aMh,^'^Mwx 








O Q q Tir9Ti«- 9 Ti*'99 

+2Su,Sa Mh^ My^x + 2Su,Sa: C2aMh^'^MwX + Co, CaeMh^^V^ 








9 Ti*'99 9 Ti*'99 9 7i^99 

-Ca^C2aCaeMh^'^V^ + Ca^CceMfe^ ^t;^ + Ca^C2aCaeMh2 V 








9 7i^'9 r»9 71^'971/r 








-Sa CaeMh/S2aVX - 2Ca SwSaMh/MwS2aV) 


11-2 






p. \/I-iiT - ci „ 111/ 


ho 




, Wf 




h2 


wt 




1 L W ^^2) 


ho 


wt 

F 




-IT^i^SwSaMh/Mw - 2SwSaC2aMh2^Mw + 2s^^SaMh/Mw 








+2Su,SaC2aMh^'^Mw + CaeMh^"^ S2aV - CaSMh-,'^ S2av) 


ho 


z„ 






ho 


z„ 


Zf 


2c„-s„l^'3 P'lj 


I, 

h2 


Zf 


Zf 


-l^i'^S^SaMh^^Mw - 2s^SaC2aMh2'^Mw + 2s^SaMh,'^Mw 








+2SyjSaC2aMh^'^ Mw + CaeA42^S2QU - CaeM/i^2s2a^;) 


h2 




z'. 


8c„ • . ^ . 


h2 




Z'f 


2i • Ca • g[ {p'^ - p'^) 


h2 


Z'f 


Z'f 
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Fields in 


the vertex 


Variational derivative of Lagrangian by fields 








+C2aCaeMh2 V + 2SwSaMh2 MwS2a " 2Su)SaMhi MwS2a) 


Ta 




\ 


eiac ■ 5cb 


' a 


' 


hi 


1 - p- A/Jt- r 

J. L-Q/ <Z- 1V±'J- A 


fa 


Tb 


h2 


1 e-Mr'Sa X 

2 Mw-Sw 


fa 


Tfj 




4 ^J,^,7S^=((l 2s^^) • (1 7&)cb 2s^^ ■ (1 + 7&)cfe) 


fa 


Tb 


Zp 


1 i-e-AU 5 

2 Mvi^-Sti, 'aft 


fa 


Tb 




ffl7ac • deb 


fn 

' a 


U'^h 
oo 


TJT" — 

w ^ 


1 Ca^-e-\/2 , 5\ , M 
4 s„ 7^)cfe7ac 


fa 


V2,b 


Wp 


-iM^^(c«.eM,v2(l - 75)^,, - 2s^Sa,^MwyZil + 7^)a6) 


fa 






IP - .S^ •a/2 /I f^x // 


1 a 


i^hib 




1 i / i\ JT /7T/-I 5n ir» Ti^" J~) / -1 , 

4 Af^^s^ (sa.eM^V2(l 7 )ab + 2s^c„,MH/y^3(l + r)ab) 


fia 


/J, J, 


\ 


e^ac ■ Scb 


H'a 




hi 


2 Mvvs„ ""-b 


fia 


Uh 


h2 


1 A^i, -.s^ [- 

2 Mws„ "ab 


fia 






4cjs„7ac((l 2S»^)-(1 7^)cb 2S^^ • (1 +7^)cfe) 


fia 


fJ-b 




1 i-e-Mfj, ^,5 

2 Mvy-Sm ^afi 


fia 


f^b 


z ^ 


Qllac ■ Ocb 


fia 


V2b 




1 Crv,,-e-\/2 /-I u 

4 s„ • (1 7^)cfe7^c 


fia 


V2b 


Wp 


4Mv^s„('^"-^^mV2(1 Y)ab 2s^Sa,Mwyn2[^+T)ab) 


j2a 


^h2b 


T IT" — 

w ^ 


1 e-Sci,,-\/2 /-I 5\ n 

4 •(! r)cb7Sc 


fia 


^h2b 


TIT" — 

Wp 


1 if i\ /r /o /1 in Ti/T D / 1 \ 

4 Mv^s„ (sa,eM^V2(l 7 )a6 + 2s^c„^ M,4/y^2 (I + 7^)a6) 


Sap 


Cbq 


w-^ 


-i^-(l-7%«c 


Sap 


Cbq 


Wp 




Sap 


Sbg 


A, 


■^aSpq^ac ■ (^cft 


Sap 


Sbq 




9s • 7afe V 
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A. MODEL IMPLEMENTATION 



Fields in 


the vertex 


Variational derivative of Lagrangian by fields 






hi 


1 Cg-e-M, X . A , 

2 Mw-Sw^PI 






Q 

'-'dp 


^oq 


h2 


1 e-Ms-Sa r X 

2 Mw-Sw "PQ 






e 

'^dp 




^/^ 




cfe (3 2su; ) • 


(1 - l%b) 


'^ap 


^bq 


Zp 


2MwSiu PI ^ab 








Sbq 




Pq i UiL, "cu 






_ 


tbq 










S(ip 


tba 


Wp 


-¥'M^f'"^Pl{Ms-{l-l^)ab 


-Mf{l + -f%b) 








1 e-V2-V^s . (1 ^5) A 






^ap 


Ubq 




1 i-e-Ms-V2-Vus fi ^.5\ ^ X 






tap 


ha 










tap 


hq 




4 ik/xv-Sif; /^y \ " ^ ' / 


Mf{l- J^)ab) 






dbq 


TIT' -I- 


4 ^ i /cu pq iu,t. 






- 

^ap 


^oq 




4 Mw^-s^ ^ ' lab^pq 








Sbq 




4 Su, ^ ' Jco^pq lau 








Sbq 




4 Mw-s™ PyV 5^1/ jao 


Mi • (1 - j%b) 






tbq 




3 ^^pq jac ^cb 








tbq 










f 

^ap 


•^bq 


hi 


1 C„-e-A/t r r 

2 A'/m.-s„ "P? "afe 






tap 


^oq 


h2 


1 e-MfSc : X 
2MwSu, PI "-^ 








tbq 




12 Cuj-Sn> PI '"•^\^ f 1 \ 


- -f^)cb - 4S^„^ • 


(1 + 7%) 


f 

^ap 


^bq 


Zp 


2MwSi,"P<} lab 






tap 


tbq 




-la'l^pqlac ■ Scb 






'^ap 


hq 










Uap 


hq 


w+ 


1 i-e-Mi-V2-Vub (-1 1 „5^ I 
4 Mws„ l-L -t- 7 ;a60p<j 






Uap 


dbq 











Ill 



Fields in the vertex 


Variational derivative of Laj 


2;rangian by fields 






Shn W+„ 


\ sj r)cb5pqlL 






U(Xp 


Shn Wt 


1 Mo-a/^-I^.o /-I , '^\ r 
LLC ivis V ^ yus I 1 «/0 \ A 






^ap 


7/L,„ A , 


2 r Mr 






Uap 


Uhn Gnr- 

^bq ^ fir 


9s ■ lab\q 






Ti 

"'ap 


"•bq ^ fjL 




(1 - j^)^b - 4.sJ 


(l + 7^)c6) 


Uap 


Ubq ^ fi 


-\9'l^pqlac ■ Scb 






T/f/+ 


W~ 7 
yy V 


-<^{p1gl^P - pPgl^- - p^g-P + pPgt^- + p^g-P 




T/f/+ 


W~ 7 

yyp 








w+ 
yy n 


T/f/~ 7^ 

yyp 


-h^{P2-P3) 






r<w+ 


o yy p, 








r<w+ 


yy p 


—i-e- My/ 






r<w+ 


fiW- A 


-e • p7 








V_ 1 L\ 


1 Ca-e-Mw 

2 Su; 










1 e-Mx4^-Sa 

2 s-tij 










Pi 








C^- Zp 


1 j-e-Miv 

2 
















nW+ 


VV p 


1 i-{l~2su,'^)-e-Mw 

2 ■ s 






VV p 


W~ Z 
yy fi 


i-e-Mw-Sw „uv 
Cw 9 






yy p 


yy fj, 


-htiP3-Pi) 








yyp z^ 


\'-^^{p'2-p'l) 




















i-e - Mw 


















hi 


1 Ca-e-Mw 

2 S7() 
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A. MODEL IMPLEMENTATION 



Fields in the vertex 


Variational derivative of Lagrangian by fields 






1 e-Mw-Sa 

2 Sj]j 








z, 










Zf 


1 i-e-Mw 

2 Syj 










s^, -Pi 








wt 

w * * p 


2 CwS^ 








(JW+ yi^- 


S-u, -^1 










1 i-e-Mw 




















1 i-e-Mw 

2 C yj " S uj 








hi 


1 Ca-e-Mw 

— 9 ^ — 










1 e-My/'Soi 








C^' hi 


2M^/ • Sa ■ g'l 








C^' h2 


-2ca ■ Mz' ^ ■ g'l 










ic e.v/2( 5) M 

4 Syj ^ ' ' '-''^ ' 










i (c«.eMeV2(l + 75)^;, - 




l%b) 


pC 

^ a 




i'^7-^(l + 7'U7^c 






^ a 




-ZTT-^ (sa.eMe \/2(l + 75)„fe 




- l%b) 


' a 


vzb 


i^-;-^(i+7')c67^c 






' a 




i (c".eM,V2(l + 75)„fe - 




l%b) 


t'^ 
' a 


J^hsb 








' a 


hob * * p 




+ 2s,nC„ Mu^v^ ( 1 


I Jab J 




V2b W^f, 










V2b W+ 




2SwSa,Mwyn2{'^ - 


l')ab) 




J^h2b W+fj, 










J^h2b 


iM,:s,.,KeM^^^(i+7^u 


+ 2s^Ca^MwyZi'^ 


- 7%b) 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 


i ct 




hi 


— ISiR. ( r . Crv ■ \Pl ■ X ■ ■ 5nh — ■ Sn ' 7Tli,h 

rjr. \^Oi ^Oij^ V ^ Jn\ ^LLO -'OL ^ Oil, Vtl 


'-'CtO / 






Z^l A> 


ho 


_£a^ 1 _ _ /2 . rj. . yD _ r _|_ g ■ Crv ■ m,,h 

y^'-'ci Oily V " ^(10 i ^cii/ ''"tz/i 


* 5n7>) 






Z/i A 




2 c„-s^ '"c 'cb 








Z/i 

i 




i • Cn • Sn, • a/2 • • 






Via 




Z'u 


i/l /tit, lcb\ OLu ^Oi) J 






VA n 




z'p 

F 








7/1 . 


^hlb 




2 X V tit/ ^cif V ^•^yni^o.b ^a^au V yni^o,b ^ ^"^a 




771 , u 1 


771 ^ 


lJu 1 u 


ho 


lif e Q ^\/27:?/^ r^„f, — ,s„ r„ ^\/27;7/^ (^^^ — 2,s„ 






ia 


^ tlLO 


z„ 


1 Ca,y-e-SQ^ ^ /^ __5 






7/1 


IJU 1 I, 


Zt? 


2 ' V ^ Un-i K'^a-v lab 'aby 








l^h 1 h 
ilLO 


Z'u 


2Cn • Sn ■ q\ Inr ■ 7^ 

^Lti/ /tit. / 






i ct 


ILLO 


z'„ 


X 'ab 








Z 




rj. (^ti ^ctjy "s/^ *^ yTi2 ^oh ^OL ^ Oil/ ^^uh 


* bnh 1 








h2 


— ^^(Sn • Cn- • \/2 • X • • Snh + Sn ' Cn • m,,h 


• ) 




za 


Voh 


z„ 








^ AO. 


Voh 


Zp 


i ■ ■ Sn ■ a/2 • ■ 

t- "-a^ •'o^ V ^ t/n2 /ab 






1^9 n 


VOh 
^ Lu 


Z'a 


i/l /tic /cov t^i^ ^ti^ y 






zct 


Z 


Z'p 

F 


i-mi,h-Sa„'^ „,5 
X 'afe 






7/0 ^ 


Uh.nl. 




i i fs^ ^c^ \/2xv^ 5 h - c„c^ "^xflxv^ 5 h + 2s^ 




TTi u/) I, 1 

"^vh^ab ) 


za 


nzo 


ho 


2 X V t^ t^*^ ^ ■^*^i/n2 tio '^cif'-'Cifiy V tfn2 Ob ^^ot 






^2a 


^h2b 


Zp. 


1 Oaij-e-Sau „,5 

2 c„-s„ >c • 7cb 






J^2a 




Zp 


-ii • \/2 • 2/;?^ . ^5^ _ ^^_^2 . ^5^) 






J^2a 


Vh2b 










J^2a 


^h2b 


Z'f 


X 'ah 






^3a 


Z^36 


hi 
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A. MODEL IMPLEMENTATION 



Fields in the vertex 


Variational derivative of Lag 


rangian by fields 




^3 a 






(Sa • Ca, • \/2 • X • ?/^3 • 5ab + Sq, • Cq • ITlyh ' ^ab) 




^ od 




z„ 


1 Ca^^-e u 5 
9 ^ „ Jac Irh 






^ od 




Zp 

^ r 








oct 






-ffi7ac7cb(ca.^ - Sa.^) 






^3 a 


00 


z'p 

F 








Uo ^ 


^h3b 






'^V2xy^.^5ab + 2SaSa^Ca^ 


rriuhOab) 




nob 


ho 






muhOab) 




^h3b 


Z. 


1 Cc^-e-Sc^ ^ __5 

-2 Cm-Sw ' 'CO 






^6a 




Zt? 








^3 a 


^h3b 


Z'a 








OCl 


llOO 


z'„ 

F 








Vh 1 n 
flLCL 


Uh 1 h 
flLO 






+ Sa ■ Ca^ ■ rriyh ■ Kb) 










^{Sa- Sa,-V^-X- ■ dab 


Cq, • Ca^ ■ TTli/h ' '^aft) 




Uh 1 

IlLCL 




z„ 


1 e-So,^'^ n 5 

Z C-w'Sw ''-^'^ 'CO 






'^rila 


^hlb 


Zp 


-i ■ Ca, • • \/2 • y^^ ■ 








Uhlb 


r- 








IlLCL 




z'„ 

F 








flA CI 








+ Sa • Ca^ ■ m^h ■ Kb) 




^h2a 


^h2b 




^{Sa-Sa,-V2-X-y^^-Sab 


Cq, • ■ 'TIT'i/h ' Kb) 




il A (I 


nzo 


z„ 


2 c„-s„ icb 






^h2a 


Vh2b 


ZjF 


-i ■ Ca, • • \/2 • yg ■ 7^ J 






^h2a 


^h2b 










^h2a 


Vh2b 


Z'f 


X 'aft 






^h3a 




hi 




+ Sa ■ ■ rriyh ■ Kb) 




^h3a 


l^hSb 


h2 




Cq • Cqc^ ■ 'Tflyh ' Kb) 
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Fields in tlie vertex 


Variational derivative of Lagrangian by fields 


nod 


noo 


1 (i-^ai, 5 


nod 


floO 




nod 




/ U ^ / 9 9\ 

-9llacYcb[Sa. - Ca/) 


'^nod 




9 

X 'db 




Ay W+o 


-e^{2gf-'''gP'' - g^^g'''' - g^^'g'"^) 




Ay Wt 


2e^ • g^^ 


A., 




2 s„ y 


An 


h^ Wt W~u 


2 s u, y 


A 




1 f - - o 111/ 

2 s„ y 


An 






A, 


w+y w-p 


Cin-C^ /o lltT 1J D UU DfT iin 1J(t\ 

'-'ID ^ tg^J.U gUIJ g^J^iy glJU g ^J, /J g 1^ U j 


An 




2 s„ y 


A, 


W-y Zf 


2 s„ y 


\ 


W+ Zy 


(l-2s„2).g2 




G„r, Gnr- G„, 
' uq ^ pT ^ (TS 


9s''{9'''9'"'fpqtfrst - 9'"'9'''fpqtfrst + 9^"" 9'" fprtfqst 
-g^'^g'^^fprtfqst + 9'"'9'"'fpstfqrt " 9^'' 9"'" fpstfqrt) 


hi 


hi hi hi 


+sa\Mh,'' + Mh/y + sjc2a{My - M^y 

+2cJsa^{My - My)s2aVx) 


hi 


hi hi /i2 


-I'^ica'iMy + My)x^ - ca^c2a{My - My)x^ 
-sa^My + uyy - sa^c2a{My - uyy 

-(1 - 2Sa^)MyS2aVX + (1 - 2So:^)My S2aVx) 


hi 


hi /i2 h2 


-^Cct'^Sci^C2a{My - + ?,Ca^Sci^{My + M/,,2)i;2 
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A. MODEL IMPLEMENTATION 



Fields in 


the vertex 


Variational derivative of Lagrangian by fields 








/ ^9 9\ Ti r 9 /I ,-1 9 9\ Ti 9 \ 

+ (1 - QCa'^Sa^)Mh^^S2aVX - (1 - 6Cq"= • Sa^)iV4/S2ara) 




hi 




1 Cr^ ^ -e^ tfv 

9^2 y 




hi 




11/9ti*'9 9 Ti*'9 9ti^9 








+Ca^C2aMh^^X + S^^Mfej^Sac?^ " S„^M/j/s2aV) 


' ''1 


hi 




Or 2 y 

^ <--l(J 


"i 


hi 
III 


Zp Zp 


11/9ti4'9 9 Ti^9 9Tir9 

-ilF^\^o^^h2 X - Ca^C2aMh/x + C„^iV4/x 








+C„^C2aM/,/x + Sa M,,2^S2q?^ " Sa^Mh/s2aV) 


hi 


hi 


Z' Z' y 


9 / 9 /II/ 

Ssa^ -g'l^- g'"' 






Z'p Z'p 










9 Ti*-9 9ti4'9 9Tir9 \ 

+Sa''c2aMh2^V + Ca^M/^^^SaaX - Ca''Mh/s2aX) 


hi 


ho 












911^-99, 9 71J-99 971^99 9 tij'99 








+ (1 - 2Sa^)Mh2^S2aVX - (1 - 2s«2)^^^ 2^^^^^^ 


hi 


ho 




2 s ^ 5^ 


hi 


ho 


Wp 


-\^{^h2 X - C2aMh^^X + Mh.^X + C2aMh^^X 
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Variational derivative of Lagrangian by fields 
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Appendix B 

Renormalisation group equations 



In this appendix we present the complete set of one-loop RGEs for the minimal U {1)b-l 
extension of the SM. For some parameters, the equations will be equal to those of the 
SM, as no extra contribution arises at one loop. 

B.l Gauge RGEs 

The RGEs for the SU{3)c and SU{2)l gauge couplings gs and g are [103] 
d 



dt 
d 
dt 
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-11 + -Ua 
3 ^ 
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22 4 1 
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3 3 ^ 6 



9 
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where = 3 is the number of generations. 

Following standard techniques, we obtain for the Abelian couplings [30; 32]: 



-dt'' 



d 

di 
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[A^^gf + 2A^^g'r9 + A^^g[f] , 

[A^'^-g^ + 2gl) + 2A^^ g[& + gf) + A^^gfg] , 



with 



/ 



(B.l) 
(B.2) 



(B.3) 
(B.4) 
(B.5) 

(B.6) 



where the first sum is over the left-handed two-component fermions and the second one 
is over the complex scalars. For the model we are discussing (Y is the SM weak hyper- 
charge, X = B — L is the B — L number), the coefficients of eq. (B.6) are, respectively. 



A 



YY 



41/6, 



A^^ = 12 , 



A 



YX 
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B.2 Fermion RGEs 

From straightforward calculations we obtain: 

d /9 2 . 2 9 2 17 2 17^2 2 ,2 5^ A 

diy' = 16^ U^* " " 4^ " ~ - 3^^^ ) ■ ^^-'^ 

For the right-handed neutrinos, it is not restrictive to consider the basis in which the 
Majorana matrix of couplings is real, diagonal and positive: y^'^ = diag (yf^, yf^, 2/3'^). 
Then we get [104; 105] ^ 

= Tt2 + [(y^)^] - 6gA , (z = 1 . . . 3) . (B.9) 



B.3 Scalar RGEs 

A very straightforward way to find the one-loop RGEs for the parameters of the scalar 
potential is to compute the one-loop effective potential and to impose its independence 
from the renormalisation scale. To one-loop level, the scalar potential V reads 

y = 1/(0) + Ay(^) , (B.IO) 

where V^^^ is the tree- level potential and Ay^^^ indicates the one-loop correction to it. 
To compute the latter it is useful to rewrite the tree-level potential 

V^^\h,x) = m^H^H + fi^ I X 1^ +Xi{H^Hf + X2 \ x 1^ +hH^H \ x P (B.ll) 

in terms of the real scalar fields 

The only combinations of fields that are involved are (j)'^ = cpf + (j)2 + (p'^ + 4>\ and 
rf = + so that eq. (B.ll) becomes 

= imV^ + + ^Ai0^ + \X,^^ + \Xs<l>W ■ (B.13) 



The one-loop correction to the tree- level potential (B.13) is, in the Landau gauge. 



AV^'H^,V) = j;(-l)2^^(2s. + m^i^^r^') 



In- 



(B.14) 



where Cj are constants that depend on the renormalisation scheme (for example, in the 
MS scheme, it is Cj = 3/2 for scalars and fermions, Cj = 5/6 for vectors). Expand- 
ing eq. (B.14) and keeping the contributions of the scalar fields (Higgs and Goldstone 



^Notice the we get a difference of a factor 3 in the third term in the RHS of the last expression in 
eq. (14) contained in Ref. [105]. The authors of Ref. [105] acknowledged the difference and will correct 
their paper. 



121 



bosons), of the top quark, of the gauge bosons, and of the RH neutrinos only, we obtain 
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where the field-dependent squared masses are, in a self-explanatory notation. 
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(B.15) 
(B.16) 

(B.17) 

(B.18) 
(B.19) 

(B.20) 

(B.21) 



As usual, we define the beta functions Pi {i = 1 ... 3) for the quartic couplings, the 
gamma functions 7m, ^ for the scalar masses and the scalar anomalous dimensions ^(i)^r) 
as follows {t = InQ): 

dt 
dm? 



dt 
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27r,??' • 


(B.26) 



Now we can extract the RGEs for the parameters of the scalar potential just by 
requiring that the first derivative of the effective potential with respect to the scale t 
vanishes, 

^y{i)^^(y(o)+AyW) = 0, (B.27) 
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keeping only the one-loop terms. Reorganising it in a more convenient way, we see that 
eq. (B.27) implies the following equations: 



m (j) 
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Imposing that each term between square brackets vanishes, we can obtain the RGEs 
for the parameters of the scalar potential after inserting the explicit expressions of the 
scalar anomalous dimensions 7^ and 7,^. The latter are easily computed and read [105; 
106; 107] 
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Inserting eqs. (B.28) and (B.29) into the RGEs, we finally obtain the RGEs for the five 
parameters in the scalar potential: 
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